APPLICATION OF CHARGE SIMULATION METHOD FOR ESTIMATION 

OF HIGH VOLTAGE FIELDS 


A t-heais Siibmllled 

In Part-1 al Fulfilment- of the Requirements 
for- the Degr-oe of 


MASTER OF TECHNOLOGY 


By 

A JAY ARORA 


to the 


DEPARTMENT OF ELECTRICAL ENGINEERING 
INDIAN INSTITUTE OF TECHNOLOGY JCANPUR 


MAY, 1 000 



CERTIFICATE 


This Is to certify that the work presented in this 
entitled "APPLICATION OF CHARGE SIMULATION METHOD 

ESTIMATION OF HIGH VOLTAGE FIELDS" submitted by A jay 
has been done under my supervision and it has not 
submitted elsewhere for a degree or diploma. 


May, 1990 


// 



Dr. R. Arora 
Assistant Professor 
Department of Elect leal Engine 
Indian Institute of Technology, 

KANPUR 


thesis 

FOR 

Arora 

been 


er ing. 
Kanpur 



E-B- /v^ 

-i;-::: :39f 

C.EiVT-'? ^ f. 

-1”'" .Vo. A)J0f}^7"‘ 


H-r 

A RO - f \ pp 


O '■■ 



ACKNOWLEDGEMENT 


I am greatly indebted to Dr. R. Arora for his perspective 
and motivation throughout the course of this thesis work. 

Thanks to Dr. P.K. Kalra for allowing me to use his PC 

Lab . 


Thariks are also due to all my friends and all those 
without whom the completion of this thesis would not have been 
poas i bl e . 


MAY, IQdO 


AJAY ARORA 



ABSTRACT 


The subject for eatiomtion o£ HV flelda gained 
importance ever since the beginninA of Electrical Engineering 
by the advent of digital computers. The conventional method 
for the estimation like field sketching by hand and by 
electrolytic tanks gave up their way. Numerical metiioda like 
FDM Si FEM found their breakthrough in this field, however, a 
revolutionary change was brought up by Steinbegler in the late 
1960 's when he Introduced CSN, a method most appropriate and 
suitable for the estimation of electric field between 
compl 1 cat ed unsymmet r leal electrode con figurations. A1 though 
this method is being used in practice in advaticed countries it 
has just made a breakthrough in India. 

In this thesis computer programmes using CSM have been 
developed for simple electrode configuration having rotational 
symmetry and involving one and two dielectrics. Also, basic 
mathematical background such as plotting of equlpot ent lal 
surfaces, for finding vertical and horizontal components of 
electric fields on arbitrarily inclined Interfaces and 
arbitrary location of simulation charges etc. has been 
developed for further work in tills field. Certain missing 
links in the existing theory have been found out. 

Thus this work can provide a good vStartlng for tackling 
complicated electrode configuration and involving 
mul t idi el ectr Ic cases. 
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INTRODUCTION 


The optimum design of insulation in high voltage apparatus 
between phase and earth is based on the knowledge of electric 
field distribution and the dielectric properties of the 
combination of insulating materials used in the system. The 
principal aim is that the insulation should withstand the 
stresses with adequate reliability and at the same time the 
insulation should not be overdimens ioned . 

The withstand voltage of the insulation of apparatus 
design with non-self-restoring insulation is determined by the 
field intensity developed at any point at the sparking gap on 
the electrodes. Corona discharge can be eliminated by 
employing properly designed high voltage shielding electrodes. 
Corona threshold voltage of shielding electrodes is a 
parameter which governs the radio and TV interferences and 
also in most cases the breakdown of the insulation between the 
electrode arrangement. Therefore, it requires a comprehensive 
study of the characteristics of corona discharges of the high 
voltage electrode system. 

The pr e- est imat i on of the electric stresses makes it 
possible to design apparatus with a high withstand voltage, 
free from rorona discharge and generating low electric 
AtrAPtflAfl ar even resulting electric field low enough to 
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provide sufficient safety inar£:n against insulation failure. 
Corona discharge magnitude and surface flux density of an 
existing electrode system can be determined experimentally 
with moderate accuracy but the pre- estimation of the electric 
stress magnitude at operating voltage can be calculated 
accurately with the help of computation before finalization of 
shield design. 

Several methods are used today to calculate electric field 
either analytically or numerically. 

For estimating electric fields, analytical solution of 
Laplace’s equation can be obtained for relatively simple 
conductor configurations. However, field distributions of some 
of the geometries that are frequently used in high voltage 
apparatus, cannot be expressed in simple analytical terms. 

As an alternative to purely analytical techniques, 
numerical methods are often used to solve such problems. These 
include the Finite Difference Method (FDM), the Finite 
Element Method (FEM), and the Charge Simulation Method 
(CSM) . 

FDM & FEM are directly based on the differential form of 
the Maxwell's equation. Solving differential equations either 
analytically or numerically involve difficulties inherent in 
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the formulation of inaccuracies arising within numerical 
procedures. These methods are also quite time consuming. On 
the other hand CSM is relatively simple as field strengths can 
be calculated analytically with reasonable accuracy because 
HV apparatus curved surfaces are generally preferred over- 
sharp edges. It is also less time consuming for many 
geometries in HV technology. 

This thesis work describes the application of CSM for some 
simple rotation symmetry fields as a numerical extension of 
image charge method (used for calculation of electric fields 
analytically for some simple configurations), thus, 
eliminating the need for formal solution of Laplace’s and 
Poisson’s equations in differential form. 



CHAPTER 1 


THEORY & DISCUSSION 

1 . 1 Image charges and their applications : 

Directly related to the application of Gauss's Lav is the 
method of images (or image charges), which could be used to 
compute analytically some important problems by means of 
ready-made solutions, thus eliminating the need for formal 
solutions of Laplace ' s or Poisson's equations in differential 
£ orm . 

The image charge theory and the CSH are based on the 
uniqueness theorem which states that a solution V(x,y,z) to a 
problem consisting of a charge distribution p(x,y,z) and 
grounded conducting boundaries through , is the unique 
solution if it satisfies the differential equation 

= -pCx, y, z)/ic 

and the boundary conditions 

VCx,y,z) =0 on [5] 

1 Z n 
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The principle of uniqueness can be applied to replace the 
effect of conducting boundaries with image charges* Some of 
the important cases are presented in the following. These 
cases are : 

1. point charge and a grounded conducting plate (or say 
Earth) of semi~inf init e dimensions. 

2. point charge and a grounded conducting sphere. 

3. point charge and a sphere at voltage V. 

4. a line charge and a conducting cylinder at voltage V. 

5. Two conducting cylinders at voltage +V and -V 
respect i vely . 

6. A conducting cylinder at voltage V and a grounded 
conducting plate (or earth surface). 

7. Two conducting spheres at voltages +V and ~V r espect i vel y . 

8. A conducting sphere at voltage +V and a grounded 
conducting plate (or earth surface); and similar few other 
cases . 

Now the thesis discusses each of the above cases one by 
one. Significance of this lies in the fact that analytical 
methods are used to determine the location as well as the 
value of charges accurately. In CSM location is generally 
assumed (surely with certain optimization criteria) and the 
charge values are determined numerically. 
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Case 1 : Point charge and grounded conductlriig plate Cor say 


Ear tlO of Semi -inf 1 nit e dimerislon CFlg. ID. 


Here the plate forms an equipot ent ial surface at 
Voltage = 0 . Thus, this plate can be replaced by a charge 
equal in magnitude and opposite in polarity placed 
symmetrically on the other side of the plate as shown in 
Fig. 1 . 

Thus voltage at any point P(r,2) would simply be 


VCr,z) = KQ 


y/r^+CZ-a)^ 


yr^ + (Z + a)^ 


and Electric field components would be 


E 

r 


^r 


and 


■ 2 KQr 


/r^+CZ-a)^ 


y^r^+CZ + a)^ 


E 

2 


-dv 

^2 


- 2 KQr 


(z-a) 

2 2 
_r +(2-a)^ 


(2 + a) 

2 2 
r +(2+a) 


Case 2 : Point charge and a grounded conducting sphere CFig.2I> 

Let the image charge of q be placed at and let it’s 
magnitude be q'. Now potential at any point (x,y,2) on the 
sphere is 
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y(x,y,z) = 0 = 


Kq 


Kq ' 




2^ 2 
) +y +z 




2 2 
) + y +Z 


2 ,, ^2 2 2 ^ , 2 ,, ^-2 2 2 ^ 
or q ((x-x^) +y +z ) = q’ ((x-b) +y +z ) 


i . e . 


2 ^ n 2 ,2 . 2 ^ ^ 2 , 2 . 2 ^. 2 , 2,2 _ 

q (x-x ) -q’ (x-b) +(q -q ’ )y +(q -q ’ )z = 0 

o 


i - e . 


2 ,2,, 2 2 2,, ^ 2 OU.2 ^2,2 22 

(q -q )(x ■+y +z ) = 2xx^q - 2xbq + b q *- ^^q 


2 2 2 2 

Comparing with x '+y +z = a we get 


2 K .2 

X q - bq ' =0 

o 


“ 2 
q 


Ci) 


and 


^2 ,2 22 

b q’ - x^q 2 

2 72 "" 
q - q ’ 


1 . e . 


,.2 ,2 2 2 

b q’ - x^qq 


2 2 2,2 

a q - a q 


(ii) 


,,2^ 2, ,2 , 2^ 2- 2 

(b +a ) q’ = (a +x^) q 
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Therefore , 


,2 , 2 
b q 


bq’ 


q 


2 2 2.2 

a q -a q 


2 2 b^q’'’ 

b q - x 


2 2 2 ,2 

= a q ~ a q 




22 ., . 2 , 2 - 

=aqCl-q’/q) 


b 


2 

q 


2 


2 2 
a q 


^ q’ 






2 

a 

F“ 


Thus the grounded conducting sphere can be replaced by a image 
charge q', placed at a distance from the centre. 


The potential at any point (x,y,z) outside the sphere 
would be 


V(x,y,z) = Kq 


and E = -V V 


1 

- Kq’ 

1 



[/Cx-b)^ + y^ + z^ 


/f ^2^ 2^ 2 

yf (x-x^) +y +z 


K ^ = 4n£ 

O 

which can also be calculated. 
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Case III : Point cliarge and a sphere at voltage V 


In this case apart from the image chage q' placed at 
one more charge = 4njc^V is placed which gives a voltage V 

on the sphere surface. 

Thus voltage and electric 
sphere can be calculated as 

V(x,y,z) = — — — " 

/ 2 2 2 
y Cx-b) +y +z 

E = -V V 

Case IV> Case V & Case VI [61 : 

A line charge and a cylinder at voltage V (Fig. 3) or two 
parallel cylinders at voltage +ve and -ve, respectively. 

Consider a line charge running parallel to a conducting 
cylinder of radius R and carrying a charge = -p^^ per unit 
length. The cross-section of the system is shown in Fig. 4. If 
the electric field of charges induced on the surface of the 
cylinder can be reduced to the field of a certain image line 


charge , by 

symmetry 

the image 

line charge p^ 

must 

be 

somewhere in 

the 

plane 

of the line 

charge pj^ and the 

axis 

0 f 

the cylinder. 

as 

indicated in Fig. 

4. The magnitude 

0 f 

the 


field at any point outside the 


^ 


r 72 2 2 2 2 

y (x-x^) +y +z y .x +y +z 
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charge and its position x with respect to the cylinder 
axis are unknown. If possible, we roust determine these so that 
the surface of the cylinder is equipot ent ial . 


Now, according to Gauss" law, the total charge enclosed by 
a cylinder coinciding with the conductor surface equals the 
flux of the electric field intensity through the cylinder , 
multiplied by £ . If the field is assumed to be identical in 
the real case and that, when the cylinder is removed and the 
image line charge introduced instead, the flux of the electric 
field intensity is also the same, it follows that p^ , that is 
to -p^ must be equal to , that is to indicated in 

Fig . 4 . 

The potential at a point M" at a distance r from a line 
charge p^ , with respect to a reference point at a distance 
from the line charge, is 


r" '’i H '^1 

2ni r 2n* 

Jo o 

r 


The distance r^^ cannot be taken to be infinite, since in this 

case there are charges at infinity. If we denote by r" and r" 

K, 

the corr esponding from the line charge p^ = "“Pj » then the 
potential at the same point due to p^ is 
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So the total potential at M' is 



(ii) 

From (ii) it is clear that the equation of any 
equ ipot ent i al surface due to the two equal line charges of 
opposite signs is r'/r = constant. Hence, if the distance x of 
the image line charge from the axis of the cylinder can be 
determined so that r’/r = constant for all points on the 
surface of the cylinder, the image line charge is completely 
determined. 

From Fig. 4, if triangles OH ' P and OP'H' are similar (if x 
is chosen appropriately, this can always be the case for a 
fixed point M’) then 

r’/r = R/b = x/R 

[b stands for (x+d), that is, for the distance from the line 
charge p^ to the cylinder axis.] Thus, equation r* /v = 
constant can be satisfied for any point on the cylinder 
surface is 

X = R^ /b . 

The Program 1. calculates for various points along the line 
joining the centres of the cylinder. The ratio of electric 
field at that point to the max field using image charge 
method . 
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Case VII & Case VIII : Two conducting spherea/A conducting 


sphere infront of a grounded conducting plate (See Fig. 4) [5], 

Some fields may be determined by the method of images 
through successive approximations. As an illustration the case 
of a charged sphere near a grounded plane has been considered. 
Both the spheres and the plane have been replaced by a set of 
point charges which will maintain these surfaces as 
equipotentials . 

First a charge is put at the centre of the sphere, as 
in Fig. 6. This makes the sphere, but not the plane, an 
equipot ent ial . Next we put the image of to the right of 
the plane. This makes the plane an equipot ent ial but destroys 
the spherical equipot ential , so we put the image of 

inside the sphere. This makes the sphere again an 
equ ipot ent i al but upsets the plane. Ue continue the process, 
which converges rapidly, until we have the required precision. 

The maximum field strength at the points H and N’ is given 
by 


max 




00 y 00 y 

y EL- + y £ 

n-o (a-x n-o C2D-a-x 

n n 
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where 


^n-1 ^ 5 D-xn- 1 ^ 


2 

and with n = 1,2,3... 


and X = 0 . 
o 

Above equations are used to calculate the field 
intensities between two oppositely charged metal spheres 
along a field line of highest field strength, i.e. between the 
shortest distance M-fl’, 


Numerical calculation of image charges has been used for 
above case and the error in voltage at the tip obtained by 
numerical technique is extremely small for a fairly large 
number of charges. 

In case of sphere and a plate, another hypothetical sphere 
is assumed to be placed symmetrically on the other side of the 
plate for the simulation purpose. 

1 . 2 CSM Theor y C 1 1 ; 

CSH as mentioned above is the numerical counterpart of the 
image charge method which is an analytical method. CSM for 
field estimation is generally applied for complex electrode 
configurations which cannot be easily solved analytically 
using methods such as the Gausse’s Law or the above mentioned 
image charges method. 
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The basic principle of CSH is very easy to formulate. 

Using the superposition principle, the potential functions of 

the fields of individual charges of any type (point, line or 

ring charges, for instance) can be found by a summation of the 

potentials (scalars) resulting from the individual charges. 

Let Q . be a number n of individual charges, and <p . be the 
J ^ 

potential at any point within the space (independent of the 
coordinate system used). The superposition principle results 
in 


1 


where p. . are the potential coefficients, which are known for 
1 J 

many types of individual charges by particular solutions of 


Laplace’s or Poisson’s equations mentioned earlier 


Uhereas the potential coefficients are known, only 

additional boundary conditions enable us to relate S . with 0. 

1 J 

quantitatively. If the individual charges are placed outside 

the space in which the field is to be computed (or inside of a 
closed metal electrode, whose surface is an equipot ent ial , the 
magnitude of these charges are related to the distributed 
surface charges which are physically bonded by the electric 
flux leaving or entering the surface of any electrode or 
conductor surrounding these charges. If n charges are 

assumed, we require also at least n known potentials to solve 
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eqn.(i) for the priori unknown charge 
easily be done by identifying the 
potentials on the surface of the 
points”), which are adequately placed 
configuration. If this potential is <p ^ 
eqn . ( i ) as 


magnitudes. This can 

potentials with n 

conductors ( ” contour 

at a given electrode 

= (p , we may rewrite 
c 


n 



P - .Q . 
1 J J 



Cii) 


This equation leads to a system of n linear equations for the 
n unknown charges 


■ ^11 

^12 

••• ^in ■ 





^21 

*^22 

P2n 

•* 

^2 



. ""nl 

Pn2 

• • • ^nn _ 




^ J 


IpliQ) = 


After this system has been solved, it is necessary to 
check whether the set of calculated charges fits the actual 
boundary conditions. It must be emphasized that only n 
discrete contour points of the real electrode system have been 
used to solve eqn. (ii) and thus the potentials at any other 
contour points considered in this calculation might still be 
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different from Therefore, eqn. (i) must be additionally 

used to compute the potentials at a number of "check points" 
located on the electrode boundary (with known potential). The 
difference between these potentials and the given boundary 
potential is then a measure of the accuracy and applicability 
of the simulation. The development and introduction of special 
objective functions in thus an important procedure within the 
optimization of the CSM . 

As soon as an adequate charge system has been adopted the 
potentials and the field strength within the space can be 
computed. Uhereas the potentials are found by superpos i t ion , 
i.e. by eqn. (i) or the corr esponding set of linear equations 
[compare with eqn. (ii)], the field stresses are calculated by 
superposit ion of magnitudes and directional components. For a 
cartesian coordinate system for resistance, the x-coordinate 
would then be for a number of n charges. 

n . . n 

"x ' 2 ' 1 '^''3 

j=l j=l 

where ^ ij are "field intensity coefficients” in the 
x-dir ect ion . 
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1 • 3 Steps i nvol ved 1 ri CSM Pr ogr anuni ng 

This CSM programming normally involves the following 

s t eps : 

Step 1 : Assumptions about the type and number of charges to 

be used i.e. point, line and/or ring charges and 
proper location of these charges within the electrode 
surface . 

Step 2 : Choice of representative points i.e. "contour points" 
on the surface of the electrode. 

Step 3 : Calculation of pot ent ial ~ coe f f i c i ent s 

i j 

Step 4 : Inversion of pot ent ial~matrix and multiplication with 
the voltage vector of the representative point to 
calculate the charge vector. 

Step 5 : Choice of "check point" on the electrode surface for 
the calculation of voltage at these points to find 
out the deviation from the actual voltage. 

Step 6 : Location optimization based on certain objective 
criteria . 

Step 7 : Choice of one of the methods for calculation of 
equi potential lines. 
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Discussion about above steps : 


Step 1 : In making assumptions regarding the type of charges 
to be used, one has to see the types of electrode surfaces ‘to 
be reproduced or represented. Generally ring charges are 
suitable for reproducing spherical surfaces and line charges 
for cylindrical surfaces while point charges are appropriate 
for various types of surfaces as shown in Fig. 5. The 
complexity of computation, in general, increases with the 
complexity of the simulation charges used. This is because 
the potential coefficients become more difficult to compute 
numer i cal ly . 


step 2 : The ” 

contour 

points” should 

be so 

chosen such 

that 

they fully take into 

account the 

cri t i 

cal points on 

the 

surfaces such 

as curves and corners 

etc . 

Larger number 

of 

contour points 

should 

be taken near 

such critical points. 



Step 3 ; 

Calculation 

of 

potential co 

knowl edge 

of analytical 

expr 

ession for vol 

discrete 

charges such 

as 

point , line 

Pot ent ial 

coef f ici ents 

for c 

ertain discret 

given in 

the following 

[1]. 



efficients involves 
tages due to various 
and ring charges . 
e charges have been 
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For Point Charges CFlg> 6> 


P. . = 
1 J 


in (r,z) coordinate systein 


: / r . 

T 1 


2 2 
4n£_/ r , + (z^-Zj) 


where C0,Zj) is the location of the point charge and (r^,z^) 
is the 'contour point'. 

This expressions for field stress components become. 


n +0 . 

2 ^ 


j = l 


r 2 ^ , >2,3/2 

4n£ [r^ + ] 


'z L 

j = l 


n Q . (z . - z . ) 

Y J i J 


/I r 2 , .2-3/2 

[r^ + (z^-Zj) ] 


For st-ralght line charge CFig. 7!> 


= 4ncCz.2 - 


In 








2 = A? 




= A? 


- (=.i - 


+ C Z - , - Z . ) ‘ 

1 Yi jl 1 




+ CZj2 - 
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the electric field components become 


” Q: f^j2'^i 

Y J — i. — 

2 4n£(Zj2"^jl^ [ 


jr"i , - 1^4^ 

r.?-'2 


^ i^'2 


and 



j = l 
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where E(k) is the complete elliptic inte£ral of the first 


kind . 

The field stress components become 


and 


n 


2 

j=i 


4n£ nr . 

1 


[r^-rj + K(k^) - E(k^) 


1 ' 1 


+(z.+zp^] K(k^)-/^2 ^^^2^ 


a,/?? 

2^ 2 



j = l 


where K(k) is the complete elliptic integral of second kind. 

Step 4 : "Check points” should be such as to correctly reflect 
the errors. The error is likely to be largest somewhere near 
middle of two "contour points”. This assumption becomes more 
and more valid if the contour points are sufficiently close. 
If this is not so, then one can choose more than one point 
between two contour points for the calculation of error. 

Step 5 : Error calculated at the check points can be used for 
the optimization of location of discrete charges within the 
electrode surface. One can also use some other objective 
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function for error minimisation. Greater is the accuracy 
required, greater is the requirement of error minimizat ion . 
One such objective function is the cummulative square error at 
the "check points”. 

i.e. optimization of U = X CV - * 

j 

For HV electrodes location-optimization. may be carried out 
by simply hit and trial method. But optimization becomes 
necessary in EHV electrodes where slightest deviation means 
meaningful errors. 

One can also assign different veightages to the errors at 
various points on the electrode, i.e. more veightage to 
crucial points and less to unimportant points by minization 
of 




where are the veightage coefficients 


In this work objective function used is the accumulated 
squared error of the electrostatic potential at the electrode 
surface. This objective function is suitable if the field 
gradient distribution between the electrode configuration is 
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the aim of the computation. On the 


other 


hand 


this 


optimization criterion may lead to greater errors if the field 
gradient at the electrode surface is the main aim of the 
calculation (Fig. 9 ). 

In general the optimization variables, which are primarily 
the position of the charges and their values, are subject to 
the following equality and inequality constraints. 

f(x.) = A; f(x.) < B; fCx.) < C; D < fCx^ < E 


where f (x^ ) can 

be any one 

of the variables x., 

1 

or , 

a linear 

or non-linear 

expr ess i on 

involving 

a number 

of 

var iabl es 

(e-e - , 

+ 5x3 = F). 

A, B, C, 

D and F 

ar e 

constants 


related to the physical system. 

The final consideration is the choice of the optimization 
technique or algorithm. The technique must be capable of 
handling highly nonlinear objective functions, equality and 
inequality constraints, and constraints that are described by 
linear or non-linear functions of the variables. Furthermore, 
it should be possible to change the constraints or the 
objective functions without modifying the optmization 
algorithm . 
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The availability of the first and second derivatives of 
the objective function U, determines whether or not gradient 
techniques that require these derivatives are suitable for 
use. A number of optimization subroutines are available in the 
FORTRAN and UATFIV scientific subroutine manuals. Rosenbrock’s 
method is one of the earliest and most reliable technique, but 
has a relatively slow rate of convergence. One of the 
f ast-converging techniques is Davidson's method as modified by 
Fletcher and Powell, combined with the Created Response 
Surface Technique of Carroll [2]. 

Other optimization algorithm generally used is given by 
Fletcher which is currently considered as one of the most 
powerful techniques for unconstrained optimization. This 
algorithm has the advantage of rapid convergence by a skillful 
use of the gradient g (where g = VU with U being the objective 
function). It has been an additional advantage due to the 
fact that as the gradient g is computed, the electric field 
intensity 2 (= ~V<^) on the desired boundary is implicitly 
computed. Therefore it can be achieved without any extra 
computation. This is a very desirable property especially for 
the applications in high voltage engineering [2], 

This method of optimized simulated charges usually gives 
surprisingly good accuracy in capacitance calculations. The 
accuracy in potential and electric fields is relatively not as 


24 



good. The accuracy in field calculations deteriorates when 
corners and edges are encountered [2]. 

The computation time depends mainly on the rate of 
convergence of the optimization method used to minimize the 
objective function. Fast converging techniques should be used 
if computer time is an important parameter. Other factors that 
could influence the computation time are the initial values of 
the optimization parameters and the effectiveness of the 
objective function [2]. 

The error minimization can either be done in an open loop 
by hit and trial method, or in a closed loop whereby the 
computer itself decides the location of the simulation charges 
for obtaining minimum error. 

Step 6 : Calculation of equl potential surfaces 


Several methods of plotting equipot ent i al surfaces have 
been considered. 

Method I : By finding out potentials at several points in the 
space around the electrode and thenjoining 
thepoints of the same potentials. But this would 
require very large number of calculations as the 
number of points would be very large, it would also 
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involve problems for plotting. Hence this method is 


considered impractical . 


Method II; Other method could be to derive an equation for 
equipot ent ial surfaces and plot these equations 


using computer 

. But in 

this 

case due 

to large 

number of charges involved 

the equation 

of 

equipot ent ial 

surfaces 

is not 

explicit in 

one 

of 

the variables 

(r or z). 

Hence 

it cannot be 

plotted 

using digital 

computers 

usi 

ng methods 

such 

as 

Newt on-Raphson 

etc. which are 

appl i cabl e 

only 

in 


case of explicit equations. 

Method III: In this method we exploit the fact that the 
electrode surface itself is an equipot ent ial 
surface. Also, we keep z constant i.e. we first 
find points on one line parallel to the r~axis and 
repeat the process along several such lines. 

Let us consider an electrode surface as shown in Fig. 10. 

Consider points a,b,c,d on the elctrode surface. Their 

potential = V . Therefore, potential at infinitsimal 

r 

distance dr from the electrode surface = V . . . 

r + dr 
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Ue know from Taylor’s series expansion that. 



If we take V ^ , 
r + dr 


V “ xV i.e. smaller by xV , then 
r r r 



Thus knowing s, V and (i^V/<^r]) , dr can be found out, 

r X , <2!!> 

i.e. the distance in which the voltage would drop by some 
given factor x. 


The new r value considered is 
value remains the same. 


r . . + (dr ) 

initial r 


Now we apply this to several points on the electrode 
surface and join the new points obtained as shown. This 
process can again be repeated at a", b’, c’, d’ and further 
new points can be found out. Thus the equipot ent ial surfaces 
can be plotted. 
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Method IV: It is known that potential along an equipot ent ial 


surface is constant. 
Also V = £(r,z) 


dV = aat/Sr') dr -f CSf/Sz) 


Along an equi pot ent i al dV = 0 


.dz _ -(df /Sr') 
" +(Sf/Sz) 


or 


dz 


(4^f /^r)/(df /i!>z) dr 


E /E dr Since f -St/Sr = E 

^ r , z \ ^ 

and -St/Sz ^ ^ j 


where E^ and E^ are electric field component at (r,z). By 

taking dr = constant value (small), we can find dz by 

calculating Electric field components at Cr,z). The new 

point is given by (r+dr, z^dz). Thus starting from a 

fixed point (r ,z ) the entire equipot ent ial surface can 
o o 

be found passing through that point. 

However, this method has one limitation. It cannot be 
applied to those problems where the equipot ent ial surface is 


28 



expected to be vertical or of negative slope at some place for 
example in rod-plane and sphere-sphere case; because (dr) 
being fixed, any small step dr from such an equipot ent i al 
would push the calculations out of that particular 
equipot ent i al surface. This method has been used only for the 
' Rogowski- el ectrodes ' case and the rod-rod case. 

Method V t The method described below is devoid of any 
limitations of earlier methods. 

An element dL along the equipot ent ial surface can be 
written as 


(dL)^ = (dr)^ + (dz)^ 


or 


(dL)' 

(dz)' 




1 + (dr/dz)' 


or 


/ 


dL/dz = ± Jl + (dr/dz) 


or 


dz 


± dL 


/ 


1 + (dr/dz)' 


Since (dz/dr) = -(E^/E^) 


(i) 


Therefor e 


± dL 


dz = 


/ 


1 + (E /E )■ 
z r 


(ii) 
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Thus, by taking (dL) to have a very small constant value, the 
equipot ent ial surfaces can be calculated by using equations 
(i) and (ii) . At the tip of the electrode, is equal to 
zero. This may give overflow error in case of equation (ii) 
when using digital computer. Thus, care should be taken to 
incorporate additional statements in the programme to avoid 
this, dz has two possible values. Correspondingly, dr also has 
two values. Thus correct value of dz should be determined by 
storing the previous point in the programme. In case of 
monotonical ly rising equipot ent ials , dz would always be +ve. 
This method has been used in the spher e-spher e case and 
mul t idi electric case. 
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CHAPTER 2 


SAMPLE PROBLEMS FOR CSM APPLICATION TO CONFIGURATIONS 
INVOLVING SINGLE DIELECTRIC 

1. Problem 1 : A conducting cyl inder / conduct or running 

parallel to the ground. The image charge theory gives the 
location and value of charge to be used. This has been 
converted into a program with certain additional features. 

Prog- 1 : This program asks the person running the program 
to give the number of charges N as input. It calculates 
the location and value of the charges given by the 
image- theory . It also calculates the error at the tip of 
the sphere. 

2. Problem 2 : Two spherical electrodes of opposite polarity 

have been taken to show the rapid fall in error with 
increase in the number of image-charges as shown in Fig. 
16 . 

Prog. 2 : The program asks the radius and voltage of the 

cylindrical conductor. It also asks for the distance of 

centre of the sphere fromthe zero potential line. It then 

calculates the position and value of image charge using 

image theory- It also calculates the ratio F = (E/E ) at 

max 

large number of points from conductor tip to the ground. 
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3. Problem 3 : Two spherical electrodes with three charges 

have been taken and CSM applied. Further, equipot ent i al 
surfaces have been plotted to check the accuracy of CSM as 
shown in Fig . 12. 

As seen earlier, atleast three charges are required 
for minimizing the error. Though the location points 
could be chosen at the image points, but in order to show 
that the error does not increase substantially even if 
charges are located at some other points; the charges were 
placed equidistant along the axis i.e. at 0.15, 0.2 and 

0.25 m. (for a sphere of 0.1 m. radius). 

Two contour points have been taken at the tips of the 
diameter along the axis of sphere and one at right angle 
to this diameter on the sphere surface i.e. at (0.0, 0.1), 

(0.0, 0.3), (0.1, 0.2). This has been done so that the 

entire sphere can be given a good representation and to 
minimize the error. 

Test points have been chosen on the entire sphere 
surface, 15^ apart from each other and error has been 
calculated on these points with voltage of the sphere 
assumed to be 1 V. 
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Prog. 3 : This program asks for the charge location 
inside the spherical electrodes and the contour location 
of points. It also asks for the radius of the sphere and 
the distance of its centre from zero potential surface. 
The program then calculates 

(a) Simulation Charge values. 

(b) errors on the electrode surface at points 15^ 
apart covering the electrode. 

(c) the cummulative square error. 

(d) Equi pot ent i al surfaces. 

4 . Two spherical tip-shaped rods of opposite polarity are 
considered. CSh has been applied using line and point 
charges. Accuracy is then checked by plotting the 
equpotential surfaces in the region of interest as shown 
in Fig- 14. 

Spherical part of the rod can be best represented by 
considering a point charge placed at the centre while the 
cylindrical part can be represented by line charges most 
appropriately. Now one can choose line charges of equal 
length, progressively decreasing length or progressively 
increasing length. To decide this, error reduction process 
was resorted. It was accompanied with successive variation 
of parameters such as, length of first line charge and 
multiplying factor for the length of line charges. In this 
way only one parameter was changed at a time while the 
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others were kept constant . 


The result 


which fiave the 


minimum error for a radius of unit length was 
length of first line charge = 0,1 

mulyiplying factor A1 = 0.6 


This is not to suggest that these are the optimum values 
but they are only one of the several possible sets of 
values for low error. Since the multiplying factor is less 
than 1, therefore, the line charges are of progressively 
decreasing length. Here 9 line charges were considered. 

Prog- 4 : This program asks for the radius of the 

spherical part of the electrodes and also the distance of 
it's centre from the zero-potential surface. It also asks 
for the length of first line segment (1) and the 
multiplying factor (Al) for successive 8 other line 
segments. The program itself then calculates the contour 
points, some of them on the spherical part and others on 
the cylindrical part. It then calculates the value of 
point and the line charges and their location. It also 
calculates the equipot ent ial surfaces. 

5. Problem S ; Rogowski -Prof i 1 ed electrode has been taken 
and equipot ent ial electrode surfaces for different gap 
spacings plotted. Ring-charges have been considered for 
this program as shown in Fig. 17. 
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CHAPTER 3 


THE MULTIDIELECTRIC CASE 

For a field space containing only one dielectric material, 
the application of the CSH to three-dimensional problems does 
not present fundamental difficulties. Even unsymmet r i cal 
electrode configurations can be treated by means of discrete 
charges [ 1 ] - 

In contrast to the simple solutions within the FDM or FEM 
for treating mul t id i el ect r i c cases, the CSM when used for 
field calculations in systems composed of two or more 
dielectrics increases the cost. This may be understood by 
considering the fundamental mathematical solutions and the 
physical mechanisms involved. The CSM is directly based upon 
physical charges and in every dielectric material polarization 
processes take place. Whereas in a homogeneous material placed 
between electrodes the absolute value of its permittivity does 
not contribute to the field strength (or potentials), but only 
the flux density D, the field distribution at the boundaries 
of different materials is heavily distorted due to the dipole 
charges at the boundaries which may not have counterparts at 
the adjacent medium. The law of dielectric refraction results 
from the physical effect and is associated with an infinitely 
thin layer of bonded charges located in the two media. The 


35 



free surface charges phys ically present due to electrical 
condition of the interface also contribute to field 

distortions, but the common dielectric refraction is not 
related to such additional charges. 

This realignment of dipoles within different dielectric 
materials must therefore be considered within the CSM. An 
exact solution with CSH must be based upon the physical dipole 
surface charge density. But continuous surfaces can also be 
simulated by discrete charges by replacing the surface charge 
density to metal electrode. This method, originally presented 
by P, Ueiss, will be presented briefly through a simple 

examp 1 e [ 1 ] . 

Figure 12 displays a cross-section of a part of an 

insulation system, in which a metal electrode with fixed 

potential, , meets two adjoining dielectric materials I 

c 

and II. The actual shapes of the two-dimensional surfaces of 
the three different boundaries (electrode-dielectric I, 
electrode-dielectric II, dielectric I-dielectric II) determine 
the optimal types of discrete charges simulating the problem. 
Thus, the localized charges 1-7 will represent point charges 
as well as intersections with line or ring charges. From 
earl'ier considerat i ons it is obvious that a part of the 
charges (nos. 1-3 denoted as n^) have been placed inside the 
elctrode, i.e. behind the metal surface. However, the same is 
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true for the charges placed on both sides of 


the 


dielectric interface (nos. A-7), beccause the influence of 

the dipolar charges within dielectric I upon the field in 

dielectric II can be simulated by the discrete charges 

nos, 4 and S within dielectric I and vice versa. A limited 

number of contour points placed at a = constant boundary 

is necessary, which is equal to the number of simulated 

charges within an electrode, and thus a number (n_ = 3) 

h, 

contour points (nos. 1~3) are adequate. For the dielectric 

interface, however, it will be sufficient for this example 

to place only two contour points corresponding to the two 

pairs of simulation charges (nos, 4 and 6, nos. 5 and 7), 

as each contour point belongs to dielectric I as well as 

to dielectric II. Equal number of charges, designated by 

n- , on both sides of the dielectric interface are thus 
Jd 

convenient and they should be placed at positions 
uniformly distributed between the mutual contour points 
and adjacent charges respectively. For our example, n is 
equal to 2 only. 

Now it is possible to set up a system of equations for our 
unknown charges based upon well-known boundary conditions. 
These boundary conditions can be subdivided into three parts 
as follows : 
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(1) The e 1 ec t rod e- d i el ec t r i c interface is a boundary with 

known potential, absolute magnitude of the 

surface charge density at this electrode is dependent upon the 

relative permittivity of the dielectric materials due to 

the polarization mechanisms in both dielectric materials. 

Since D = sE = s s- E, where s is the permittivity of vacuum, 

r o o 

the absolute magnitudes of our simulation charges would depend 
upon these material charact erist ic parameters. It is necessary 
to take these physical effects into account, which are 
included within the potenial coefficents. For any homogeneous 
dielectric material, the electric field may be computed 
independent of relative permittivity ^ , and the potential 
coefficients are in general always computed by assuming 
The absolute magnitudes of the discrete charges used within 
our system are based upon a superposition of potentials. And 
thus we can use the known potential at the electrode interface 
to derive two sets of equations based upon n^ contour points 
taking only dielectric I into account , for which the charges 
within dielectric II can be neglected : 



Cl~3) (4-5) 


Using eqns . (iii) and (iv) subject to two new boundary 
conditions, the electric field within dielectric II could be 
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computed. All Qj charges within eqn . ( i ) , which are not yet 
known, define the potentials within this material. 


For the computation of the field distribution within 
dielectric I, the same considerat ions are applicable. But now 
we neglect the charges within dielectric I, which resultsin an 
equal set of three or n^ equations, 


n. 




I".'., • z 


Q.P. . = 4 > 

J 1 J c 


(ii) 


j = l 
( 1 - 3 ) 


j =ng + ng+ 1 


(6-7) 


(2) The potential at the dielectric interface is unknown. Ue 
know, that due to the continuity of the potential at either 
side of the interface, the potentials must be equal at each 
contour point. As the charges within the electrode (nos. 1-3) 
will not disturb the continuity condition, the potentials due 
to the charges within the dielectric materials must satisfy 
the condition. 


n. 


j = l 
(4-5) 




• I 


Q .P. . = 

J J J C 




(6-7) 


(iii) 


39 



This equation refers to a number of (=2) contour points, 

o 

giving an equal number of new equations. In these equations 
charges are involved, which have not yet been used within 
eqn . (i) or eqn . (ii) respectively. It should be noticed that 

this potential continuity condition implies that the field 
stress components tangential to the interface are equal. 


3. Finally, the third boundary condition refers to the 
continuity of the normal component of the electric flux 
density crossing the dielectric interface or the discontinuity 
of the normal components of the field intensity. To include 
this condition, the ”field intensity coefficient” f. . must be 

1 j 

considered, which is the contribution of the charge j to that 

component of the field vector, which is normal to the 

dielectric boundary at a contour point i. Then for any normal 

component (E ). = Q.f. this condition may be written as 
ni jij 




Iv. ^ I 


Q.f.. 

J ij 


j=l 
Cl-3) 




(6-7) 


1 1 


n. 




Y' Q .f . . + V Q.f.. 


j=i 
(1-3) 


j=nE+l 


(4-5) 


(iv) 
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where 

and j j are 

the 

permittivities of the 

two 

di el ectr i cs . 

This equation 

ref ers 

again to a number of 

"b 


contour points, and thus a total number of (n^+Zn^) linear 
equations are given for the calculations of the same number of 
unknown charges. This procedure demonstrates the difficulties 
involved with the implementation of dielectric boundaries, as 
a significant number of additional charges increase the 
computational efforts. 
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CHAPTER 4 


PROGRAMMING FOR THE MULTIDIELECTRIC CASE 

Thia program ia quite different from the proceeding 
programa becauae of the preaence of aurface charges on the 
dielectric. Eaaential detaila of the program are given below : 

Charge locations ; The program involves three sets of 
charge locations. First, those simulation charges which are 
inside the electrode; second, those simulation charges which 
are in the air near the air dielectric interface; third, those 

which are inside the dielectric near the air dielectric 

\ 

interface . 

Cb3 Contour points : There are tliree .sets of cotitour polnt.a. 

First, those on the electrode air interface; second, those 
on the electrode dielectric interface; third, those on the ait- 
di el ectr i c Interface. 

Cc!) Potential matrix : The potential matrix has coefficient 
arising out of the following reasons : 

(1) due to potential on the air electrode Interface. 

(li) due to potential on the electrode dielectric 

interface. 

(lii) due to equating of potentials on the air dielectric 
interface . 

(iv) due to equating of vertical components of the 
electric field on the air dielectric interface with 
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the relative permittivity taken into account. 


Cd> Matrix inversion and multiplication with a vector of 1 " s 
and O’^s at appropriate positioris. 

Ce2> Error calculations : Calculates the error at desired 

points on the electrode surface. 

Cf3 Equi potential surface calculation : Calculates the 

equipot ent ial surface passing through a desired initial 
point in air or dielectric. 

CgD Tangential field calculation : Calculates the tangential 
field at the desired points on the interface between air 
and dielectric. 

The program asks for the various charge location, contour 
points on the electrode, the coefficients of the two 
boundaries, the relative permittivity of the dielectric with 
respect to air, the points on the electrode surface where 
error calculations are desired, the points on the interface 
where tangential field calculations are desired, the backup 
and the initial point of the equipot ent ial surface. 

The above programme has been run for a sample problem 
namely a spherical electrode of radius = 10 cm., with its 
centre 20 cm. above a dielectric block of height 10 cm. and 
width 20 cm. as shown in Fig. 19. 

The equipot ent ial surfaces plotted by using the programme 
have been shown in Fig. 20. 
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CONCLUSION 


(a) In practice, the appropriate shaping of electrode 

configurations within insulation systems is an essential task, 
as field stresses may well be reduced and kept low by this 

method. Such electrode shape optimization techniques are 
either based upon an iterative process in which the contour 
points are shifted after each computation of the field 

stresses or are based upon a superposition of fixed simulation 

charges representing the original system and additional 
"optimization” charges by which the field distribution 
changed due to a given objective function [1]. 

(b) The application of discrete simulation charges used in 

the charge simulation technique provide at least a very 

reliable and efficient method to solve many two- and 
thr ee-dimensional problems. However, it should be recognized 
that the option of surface charges at electrodes or 

dielectric boundaries, i.e. distributed layers of charge 
sources in free space, offers definite advantages, because in 
this simultation method the role of physical charges, which 
are the origin of electric flux densities, is taken directly 
into account. 
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The laain contributions made in this thesis are as follows: 

(a) CSM has been presented in a unified manner. Evolution of 
CSM from image field theory has been successfully presented. 

(b) The thesis has presented successful programs for 
parallel conductors, spher e*~sph er e and rod-plane electrode 
configurations involving single dielectrics. 

(c) Important basics for further work on CSM have been 
developed namely the mathematical theory for plotting the 
equipot ent ial . The missing links in the theory of 
mul t i el ectr ic case with rotational symmetry have been brought 
out and the problem of programming the mul t i el ectric case has 
been successfully developed. The problem of incorporating the 
vertical components of field on any arbitrarily inclined 
interface and from any arbitrarily placed charges has also 
been solved in the program. Apart from this the program has 
also solved the problem of equipot ent ial surface plotting at 
the interface and the proper choice of sign of dz. 
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SUGGESTED WORK 


(1) The program on Mul t i di el ectr i c case which this thesis 
has presented involves only point charges inside as well 
as outsidethe electrode surface. The program can be 
further extended to include line and ring charges within 
the electrode surface also. 

(2) The program is valid for axially symmetric cases. 

Programs involving electrode with non-axial symmetry can 
also be developed. 

(3) The charge location has essentially been done by hit and 
trial method which is an open loop programming process. 
For higher accuracy it is important to adopt closed loop 
method of error minimization which will have to be 
incorporated into the program. 

(4) The charge location optimization criteria used in this 
program is the cummulative square error. Some better 
methods can be used as suggested in Chapter 2 in this 
thesis. 
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APPENDIX 


The realization of homogeneous fields within a finite 
volume of Insulating material is very difficult- Using 
parallel metal plates of limited dimensions creates the 
problem of a proper stress control at the edges of the plates. 
The field problem becomes thus three-dimensional, though a 
rotational symmetry exists if the parallel plates are circular 
discs . 

Depending upon the material to be tested, the breakdown 
strength may be very sensitive to local high fields within the 
whole electrode arrangement. Therefore, the highest stress 
should only be present in the homogeneous field region, where 
the plates are in parallel. A certain profile of the 
electrodes is necessary outside the plane region to limit the 
dimensions, but the field strength at the curved edges should 
never exceed the value E = V/d, if V is the applied voltage 
and d the distance between the parallel plates, Rogowski 
proposed electrodes for uniform fields for axially symmetrical 
systems whose profile follows the analytical function first 
introduced by Maxwell, 

z=~Cw4l-»‘e'*^) (a) 

where z and w represent the complex coordinates in the z- and 
W“planes. Substitution of the coordinates for the complex 
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values z = x+iy and w = u+iv and separation of the real and 
imaginary parts gives 

a u , 

X=: — (u + 1 e cos v) ; 

fT 

y — (v •+ e sin v). (b) 

n 

Aseuinlne two Infinite, parallel ’’plates” in the w-plane, 
the coordinates of which are given by v = ±n = const, it can 
be recognized from eqn. (b) that these plates are transformed 
into the z-plane to the left half-plane only. All other lines 
V = const with -n < V < +n can be assumed to be other 
equipotential lines, and all lines u = const with -oo < u < +oo 
can be assumed to be field lines in the w-plane, representing 
a uniform field distribution. These lines appear in the 
z-plane as shown in Fig. (a) providing the electrical field 
distribution of parallel plates terminating at x = 0. The 

concentration of the equipotential lines, v = const, within 
the z-plane may well be recognized at, or in the vicinity of, 
the edges of the plates. 

The parallel plates, v = ±n , are thus inadequate to 
fulfill the demand for field distribution whose intensity is 
limited to the field strength within the homogeneous part of 
the arrangement, i.e.for u ^ -n . It is obvious that the field 

is 

strength along equipotential lines for which -n < v < +n 
provide better conditions. For quantitative assessment the 
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field atrength within the z-plane may be computed in several 
ways , as shown : 


From the conjugate complex field strength in the z-plane 


E = E - i E 
z X y 


. dw „ . 1 

dw 


(C) 


the absolute values could be computed by |E | = /T2 | ^ 

z ' b + b 

A second possibility is given by 


E = E +iE = -grad v = 
z X y 



h 

< 

1 

[dxj 

H^JJ 


Cd) 


which needs a partial differentiation only. 


Finally, the absolute value of E may be computed by 

z 


£y 2 


Ce) 


a method which is easiest to apply to our separated analytical 
function, eqn . (b). Combining eqns . (b) and (e), we easily may 
find the field atrength as 


n 


/7~ 2u~i 

1 / 1 + e +2e cos v 


f (u;v) 
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PROGRAMMES 



PROCRAM PARCYL; 

CONST K«18.0E0S; 

VAR 2 1 i V , ZC , R I 20 , 2 , L , PO /EM , VZ , E2 , F , D2 : REAL j 
U, I .'INTEGER; 


FI ,FE:TEXT; 

BEGIN 

ASSIGN(F1, 'PAR.PRN' J; 

REWRITE<F1); 

WRITELNC 'VOLTAGE, CENTRE, RADIUS ARE'); 

READLN(V,2C,R) , 

URITELNCFt , 'VOLTAGE OF THE CYLINDER*',V:8:4); 

MRITELNCFI , 'DISTANCE OF CENTRE FROM THE PLATE* ' , 2C : 8 : 4 ) ; 
WRITELNCFI , 'RADIUS*' ,R:8:4); 

WRITELN; 

20:*SORTCSQR(ZC)-SQRCR) ) ; 

WRITELNCFI 'LOCATION OF LINE CHARGE FROM THE PLATE- ' , 20 : 8 : 4 ) 
21 :-{2C-R) ; 

L:»LN( (20+21 )/(Z0-Z1 ) ) ; 

PO:-V/(K*L); 

WRITELN(F1 , 'CHARGE PER UNIT LENGTH -',P0); 

WRITELNCFI , ' *»«»«»=»«*««««*«»*«*»««**»««*«**«**«**«*«**«*«*» 
WRITELNC 'NUMBER OF POINTS BETWEEN CYLINDER AND PLATE 
READLN(U) ; 

WRITELNCFI , 'NUMBER OF POINTS BETWEEN CYLINDER AND PLATE- ',U 
D2:*(2C-R)/U; 

EM:«K*po'kC C1/(Z0“ZC+R) )~< 1/C20+2C-R) ) ) ; 

WRITELNCFI ); 

WRITELNCFI , 'MAXIMUM FIELD INTENSITY «',EM:8:4); 

WRITELNCFI ) ; 

FOR I:*0 TO U DO 

BEGIN 

2:»I*0Z; 

WRITELN CFl , '2-',Z:6:4); 

VZ:*K%PO*CLNC20+Z)-LN(ZO~2) ) ; 

WRITELNCFI , 'VOLTAGE AT 2*',V2:8:4); 

EZ:*K*PO*C ( 1/(20~2) >-C 1>(20+Z) ) ) ; 

WRITELNCFI , 'ELECTRIC FIELD AT Z« ' , EZ : 8 : 4 ) ; 
f :»EZ/EM; 

WRITELNCFI , 'F AT Z*',F:8:8J; 

WRITELNCFI ) 

END; 

END. 



VOLTAGE OF THE CYLINDER*! 0000 . 0000 
DISTANCE OF CENTRE FROM THE PLATE* O.EOOO 
RADIUS* 0.1000 

LOCATION OF LINE CHARGE FROM THE PLATE* 0 . 1 73S 
CHARGE PER UNIT LENGTH « A - £1 84762083E-OT 


BKXKSSSSBSSVa 


NUMBER OF POINTS BETWEEN CYLINDER AND PLATE* £0 


MAXIMUM FIELD INTENSITY -7593£.57t8 


2 - 0.0000 

VOLTAGE AT 2* 0.0000 

ELECTRIC FIELD AT 2* 0.0000 

F AT 2*0.00000000 


2* 0.0050 

VOLTAGE AT 2*438.5187 
ELECTRIC FIELD AT 2*2533.1967 
F AT 2*0 . 033361 13 

2* O.OtOO 

VOLTAGE AT 2*877.7700 
ELECTRIC .FIELD AT 2*5079.1018 
F AT 2*0.06588963 

2* 0.0150 

VOLTAGE AT 2*1318.4936 
ELECTRIC FIELD AT 2*7650.6370 
r AT 2*0 . 1 0075567 

2 * 0.0200 

VOLTAGE AT 2*1761.4443 
ELECTRIC FIELD AT 2*10261.1583 
F AT 2*0 . 13513514 

2* 0.0250 

VOLTAGE AT 2*2207.3998 
ELECTRIC FIELD AT 2*12924.6931 
F AT 2*0. 17021277 

2« 0.0300 

VOLTAGE AT 2*2657.1691 
ELECTRIC FIELD AT 2*15656.2004 
r AT 2*0.20618557 

2* 0.0350 

VOLTAGE AT 2*3111 .6020 
ELECTRIC FIELD AT 2*18471.8680 
F AT 2*0.24326672 

2* 0.0400 

VOLTAGE AT 2*3571.5995 
ELECTRIC FIELD AT 2*21389.4568 
F AT 2*0.28169014 

2* 0.0450 

VOLTAGE AT 2*4038.1259 



ELECTRIC FIELD AT Z«S44E8.7094 
F AT Z«0.3£!7158£ 

2 « 0.0500 

VOLTAGE AT 2*451 £.££85 
ELECTRIC FIELD AT 2*87611.8443 
F AT 2*0.36363636 

Z « 0.0550 

VOLTAGE AT 2*4985. 0£39 
ELECTRIC FIELD AT 2*30964. 1687 
F AT 2*0.40776499 

2 * 0.0660 

VOLTAGE AT 2*5487.7775 
ELECTRIC FIELD AT 2*34514.8053 
F AT 2*0.45454545 

2 * 0.0650 

VOLTAGE AT 2*5991.8671 
ELECTRIC FIELD AT 2*38897.7084 
F AT 2*0.50436469 

2 * 0.0700 

VOLTAGE AT 2*6508.8410 
ELECTRIC FIELD AT 2*48358. 8£69 
F AT 2*0.55776898 

2 - 0 . 0750 

VOLTAGE AT 2*7040.4488 
ELECTRIC FIELD AT 2*46727.7365 
F AT 2*0.61538468 

2 * 0.0800 

VOLTAGE AT 2*7588 . 66£4 
ELECTRIC FIELD AT 2*51479.7097 
F AT 2*0.67796610 

2 * 0.0850 

VOLTAGE AT 2*6155.8385 
ELECTRIC FIELD AT 2*56678. 5387 
F AT 2*0.74643849 

2 * 0.0900 

VOLTAGE AT 2*8744.5855 
ELECTRIC FIELD AT 2*68410.3389 
F AT 2*0.88191781 

2 * 0.0950 

VOLTAGE AT 2*9358.0615 
ELECTRIC FIELD AT 2=68788.7885 
F AT 2*0.90584089 

2 * 0.1000 

VDLTACF AT 2*10000 0000 



PROGRAM SSERIES; 
yB»s crt; 

CONST K-9.0E09; 

|jl*3. 1419; 

VAR B,V,R, 00.01 ,XNI.QNI.0T.E1.XNF.QNF.C,X1 .V0.V1 . VT , EROR : REAL ; 
ch : char ; 
n,l -.INTEGER; 

FI, FS: TEXT; 

BEGIN 

ASSXGNtFl, 'SS.PRN' ); 

REMRITE«F1 ); 
clrBcr; 

'Give the distance between the centres of spheres'); 
REAOLNiB); 

writeln( 'Give the voltage of the spheres'); 

READLN(V); 

nrltelnt 'Give the value of Radius of the spheres'); 

REAOLN(R); 

writelnt ^*«x»* = s*xa»«a» = s=:e = =: = = *« = s= = a:s**« = «B ' ) ; 

IfRITELNif 1 , ' Distance between the centres of sphere* ',8:12:6); 

writelntft,' Voltage of the sphere* ',V:1£:6); 

writelnift,' Radius of the sphere* ',R:14:6); 

writeln(f1 , * «**«x»xxx«==**«**««««**»«**««***»» ' ); 

writelnC f 1 ) ; 

writelnt 'Press any key to Continue'); 
ch :*reaiikey ; 

writelnt 'Give the no of Charges'); 

READLN(N); 

WRITELNCFt , 'NUMBER OF CHARGES* ', N ) ; 

O0:-ReV/K; 

ilRITELNCFI , 'Q0«',Q0); 

V0:*V-(K*Q0/(B-R) ) ; 

WRITELNCF1 ); 

XI :«SQR(R)/B; 
yRITELNCFl , 'XI*' ,X1 ) ; 

Q1 :«Q0*R/B; 

WRlTELNCFt , '01*' ,Q1 ) ; 

VI ;«K*Qie( ( 1/tR-XI ) )-( 1/ (B~X1 -R) ) ) ; 

XNI:«X1 ; 

QN1:*Q1; 

WRITELNCFI); 


OT: *00+01 ; 

VT:«V0+V1 ; 

El :»K*t CaO/SQR(R) l-tOO/SORCB-R) ) + (Q1 /SQR(R-X1 ) )-( Q1 /SQR ( B-R-XI ) ) ) 

FOR I: *2 TO N DO 

BEGIN 

XNF:*SOR(R)/<B-XNI >; 

WRITELNtFI , 'X*' ,XNF) ; 

QNF:-QNI*R/(B-XNI); 

WRITELNtFI , '0*' ,QNF) ; 

VT *VT+K*QNF+( (1/(R-XNF) )- ( 1 /(B-XNF-R ) )); 

QT:*QT+ONF; 

El : »EI +K+ ( ( QNF/SQR ( R-XNF ) ) - ( QNF/sqr C B-R-XNF ) ) ) ; 
xni :*xnf; 
qni :*qnf ; 

WRITELNtFI ) 

END; 

WRITELNtfl , 'TOTAL CHARGE ON SPHERE* QT ) ; 

WRITELNtFt . 'ELECTRIC FIELD AT THE TIP*', El); 

C:*OT+0.5/V; 

WRITELNtFI , 'CAPACITANCE* ' , C) ; 

EROR : * C ( V- VT ) /V > 1 0 Q ; 

WRITELNtFI , 'ERROR AT THE TIP*', EROR); 

CLOSE (FI ) 

END. 



0.400000 


Oistanc# between the centres of sphere* 
Voltage of the sphere* 10000.000000 
Radius of the sphere* 0.100000 

at iK s ac s «E 3E » » ae s 3c s s 3 s: ss as s as s « s s s s SB a: s 3B s =s s 


NUHBER OF CHARGES*10 
Qo* 1 . m n n 1 1 1 e-ot 

X1« S.5000000000E-08 
QI* C, 7777T77776E-0© 

X« a.B6066S6667E-0E 
Q« 7 .4074074074E-0? 

X« E. 678571 4e8feE“0E 
Q* 1 .9841E69841E-09 

X« E.6794258373E-02 
Q* 5.316321 1 057E~1 0 

X* 2. 6794871 795E-02 
1 . 42430 1 4C45E- 1 0 


X* 2. 6794915836E-02 
Q* 3.61693957782-1 1 

X« 2.6794918998E»02 
Q« 1 . 0227458681 E-n 

X* 2.6794919225E-02 
Q* 2.7404392924E-12 

X* 2. 679491 9242E-02 
Q* 7.3429849526E-1 3 

X* 2. 6794919243E-02 
0I» I . 967546e881E-13 

TOTAL CHARGE ON SPHERE* 1 . 4900657388E-07 
ELECTRIC FIELD AT THE TIP* 1 . 4575598322E+05 
CAPACITANCE* 7 . 4S032S6939E-1 2 
ERROR AT THE TIP* 6 . 4816325903E-05 



PROCRAM MATCOFF; 

CONST N«3j m«l; 

typ* »at»arrBy Cl..n,1..n3 of real; 
cat*array E1..n/1..23 of raal; 
rat*array C1..n#1..m3 of raal; 

VAR RZ:cat; 

PQ : cat ; 

charga :rat ; 

NSTtmat; 

ZC,R,PJ^QJ, VI, VJ,DH,DL,RI,ZI , THETA , ER , EZ , DZ , DR , EROR , ERROR , A , B , C , D : REAL 
LI , II , J1 , THETA J,THETAO, VI ,VM,C, U1 .INTEGER; 

NSTl ,N8T£:REAL; 

I , J , K , L , S , u , V , T : I NTEGER ; 

Ft,FH,F3:TEXT; 

procadura inverse ( a ; liat ; VAR b. rat); 

CONST NM*50; 

N*3; M*1; 

VAR IPIV: ARRAYCt . .NK3 OF INTEGER; 

INDXR -.ARRAY Cl . .NM3 OF INTEGER; 

INDXC: ARRAY Cl . .NM3 OF INTEGER; 

BIG,DUM,PIVINV:REAL; 

I, J, V,U,LL,K,L,P,Q, ICOL, IROW: INTEGER; 

BEGIN 

C FOR P;«1 TO N DO 
BEGIN 

FOR Q:at TO N DO 
BEGIN 



RE'AD(ACP,Q3 ) j 
END; 

REAOLN 

END;! 

FOR P:*l TO N DO 
BEGIN 

FOR Q:-l TO M DO 
BEGIN 

BCF,Q3 :■! . 0; 

END; 

END; 

FOR J:*l TO N DO 
BEGIN 

IPIVCJI :«0; 

END; 

FOR I:*1 TO N DO 
BEGIN 
BIG:«0.0; 

FOR J:*1 TO N DO 
BEGIN 

IF C IPIVCJ3 <>1 ) THEN 
BEGIN 

FOR K:«1 TO N DO 
BEGIN 

IF (IPIVCK3«0) THEN 
BEGIN 

IF C ABSC A CJ,K3 ) >*BIG ) THEN 
BEGIN 

BIG:«AB5(ACJ,K3 ) J 
IROW;“J; 

I COL : »K 
END; 

END 

ELSE 

IF CIPIVEK3M ) THEN 
BEGIN 

WRITELNC 'SINGULAR POINT'? 
READLN 
END; 

END 

END 

END; 

IPIVCICOL3 :«IFIVCICOL3+t ; 

IF CIROU <> ICOL) THEN 
BEGIN 

FOR L:*l TO N DO 
BEGIN 

DUf1;»ACIROW,L3 ; 

AIIROWiLa :*AC:iCOL,L3 ; 
ACIC0LiL3 :»DUH 
END; 

FOR L:«1 TO M DO 
BEGIN 

DUM:*BCIR0U,L3 ; 

BCIROW.LI :«BCICOL,L3 ; 
BCICOL,L3 :«DUM 
END; 

END; 

INDXR Cn : = IROW; 

INDXRCI3 :«ICOL; 

IF (ACICOL, ICOL3=0.0) THEN 
BEGIN 

WRITELNC 'SING HATRIXM ; 
READLN 
END ; 


PIVINV:*(1 .O/ACICOL, IC0L3 ) ,* 

ACICOL, ICOL3 :«l . 0; 

FOR L:-1 TO H DO 
BEGIN 

ACICOL,L3 :»Ai:iCOL,L3»PIVINV 
END; 

FOR L:«1 TO « DO 
BEGIN 

BCICOL,L3 :«BCICOL,L3*PIVINV 
END; 

FOR LL:*1 TO N DO 
BEGIN 

IF (LLOICOL) THEN 
BEGIN 

DUM:«ACLL, ICOL3 ; 

A ELL, ICOL3 :«0.0; 

FOR L:*t TO N DO 
BEGIN 

ACLL,L3 :«ACLL,L3-ACICOL,L3*DUH; 

END; 

FOR L:«1 TO H DO 
BEGIN 

BCLL,L3 :-BCLL,L3-BCICOL,L3*DUri 
END; 

END; 

END; 

END; 

FOR L:«N DOWNTO t DO 
BEGIN 

IF (1NDXREL3 0INDXCEL3 ) THEN 
BEGIN 

FOR K:*l TO N DO 
BEGIN 

DUIi:«ACK ,INDXREL3 3; 

ACK, INDXRCL33 :«ACK, INDXCCL33 ; 

ACK, INDXCCL33 :=DUM 
END; 

END; 

END; 

C FOR U:«1 TO N DO 
BEGIN 

FOR V-l TO n DO 
BEGIN 

WRITEtBCU, V3 :8:3) 

END; 

WRITELN 

END;> 

END; 

C»ae«»aBaE*aB»ssa:3s**ac»5=s» = *e»END OF PROCEDURE INVERSE *«*»«*»«***«*«»*»»*»■**«*»« 3 
BEGIN 

ASSIGNCF! , 'S31 .FRN' ) ; 

REURITE(F1 ) ; 

WRITELNt 'TEST POINTS 
WRITELNCFI , 'TEST POINTS '); 

FOR 1:*1 TO N DO 
BEGIN 

FOR J:*t TO 2 DO 
BEGIN 

READ(PQCI, J3 ) ; 

WRITECFI ,PQCI , J3 :8:4) ; 

END; 

READLN; 

WRITELN(F1 ) ; 

WRITELNCFT) 

END; 

READLN; 






JE LOCATIONS ARE 
lARCE LOCATIONS ARE 'J; 

DO 

DO 


,L3 




g.jBE5BSJBs«;ssi**ss**a.*s*as=sasa*ssKa*a»i:!a:*a:asx!as*se*sBaKata:*«»*s**»asa6a!»as * } ; 

ISTANCE MATRIX IS'); 

DO 

DO 

(SQR(PQCS, 13-R2CT, 13 ) +SQR (PQ CS , 23 -RZ CT , 23 ) )) ; 
CSQRCPQCS, 13-R2CT, n )+SQR (PQ CS , S3 +R2 CT , 23 ) ) ) ; 

1-W5T2; 

S,T3 ;8:4) ; 

) ; 


***CALLINC THE PROCEDURE INVERSE ♦♦*♦♦*♦♦**♦♦♦**♦#♦**♦*♦*> 
Charge) ; 

KSKSBBSSESEasKsxsssEssszsisssaEssssssssssBESSsaBSsssssasssssscxacsssestaBKsaES ' ) * 

harge vector is'); 

I ' do 

i do 

large Cu, 13 : 12:8) 


tsaesxaBaEsaEacaissssiasssasraaexsxsxsasiBsszxsarsssrsssxssaBssescxssassssaBBSsacaEsacs ' ) ) 

C A L CUL A T I ON ***«»=»**»»*»*« ®* *« > 

CENTRE « RADIUS ARE '); 
f, ' CENTRE & RADIUS ARE '); 

R); 

I, 'CENTRE*' ,2C:8:4, ' RAD lUS* ' , R : 8 : 4 ) ; 

5; 

31 V THETAO); 

1 ) J 

. iy. 

TO (V«+1> DO 

yi-l )eTHETAO; 

'THETA«',THETAJ:8); 

THETA*', THETAJ:8); 
tTHETAJ»3. 142/180); 

COS ( THETA J*3. 142/1 80 ) ; 

'PJ»',PT); 

I , 'OJ*' , QJ) ; 


TO N DO 


EELl , I3*< < l/SORTCSQRCRZai , 13 -P J ) +SQR ( RZ CL 1 ,23-QJ) ) ) -< I /SORT ( SQR C R2 CL I , 1 3 +P J ) +SQR ( RZ CL1 ,23+QJ)) ) ) 



WRITELN(F1 , 'VOLTAGE AT POINT PJ,QJ IS',VJ:8:4); 

EROR:«{ 1-VJ) ; 

ERROR :-ERROR+SQR(EROR ) ) 

URITEln(F1 , 'ERROR AT P J , Q J* ' , EROR : 8 : 4 > ; 

WRITELNCFl , ' ' ) 

END; 

WRITELNCFl ) ; 

WRITELNCFl . 'CUHMULATIVE SQUARE ERROR* ' , ERROR : 8 : 8 ) ; 

WRITELNCFl); 
writeInCFl ); 

Writ#lnCFI , '*»=:**:***aasa*a*a:**a**»««*»»a*«»a»*«»«a»»***«*»»»««*«*«»» 
«rit#lnCF1 , ' **:«**ffiS5a:*»**a»equipot#nt ial «urf*ce calculat ion********® 
0H:«(ZC-R)/5; 

DL:«R/T.5; 

FOR G:»1 TO 5 DO 
BEGIN 

RI:«0.0005; 

2I:*(G)*DH; 

WRITELNCFl , ' R I* ' , R I : 8 : 4 , ' Z I ® , Z I : 8 : 4 ) ; 

WRITELNCFl , ' ); 

WHILE RI>0.0 DO 

BEGIN 

ER:*0.0; 

EZ:»0. 0; 

FOR U1;*1 TO N DO 
BEGIN 

A:«C 1/(SQR( RZ CUT , I 3-RI )+SQRCRZ CUl , 23 -21 ) ) ) ; 

B :«( l/CSQRCRZCUl . 1 3 -RI ) +SQR ( R2 CUl , 23-i'2I ) ) ) ; 

ER:«ER+CHARGECU1 , 13*2^{A-B); 

C:«( (ZI-R2CU1 , 23 ) / ( SQR ( R I -RZ CU 1 . 13 ) +SQR C ZI-RZ CUt , 23 ) ) ) ; 

D:-( (ZI-fRZCUl ,23 ) / ( SCJR C R I-RZ CU 1 , 13 )+SQR ( ZI+RZ CUt , 23 ) ) ) ; 
E2:«EZ+2*CHAR6ECUI , t 3 * ( C-D > ; 

END; 

WRITELNCFl , 'ER=' ,ER : 12:4, ' E2= ' ,E2: 12:4) ; 

WRITELNCFl , ' 

IF ER»0.0 THEN 

BEGIN 

D2:«0.0; 

DR : »DL 
END 
ELSE 
BEGIN 

DZ : «DL/SQRT C I +SQR C EZ/ER ) ) ; 

DR :*-CEZ/ER)»D2 
END; 

Z1:*2I+D2; 

RI :«RI+DR; 

WRITELNCFl, ' R I» ' , R I : 8 : 4 , ' 21® ' , 21 : 8 : 4 ) ; 

«rlt«lnC ' ' 

END; 

WRITELNCFl); 

WRITELN CFl) 

END; 

WRITELNCFl , ' » = asaas: = = *a*aa* = aa« = s = = xa:3:a = «:eas«=t»*s:*«=:aECxaa*«»=aE = *a*aB ' ) 

END . 




TEST POINTS 

0.0000 0.1000 

0.1000 0.2000 

0.0000 0.3000 


CHftPCE LOCATIONS ARE 
0.0000 0.1500 

0.0000 o.eooo 

0.0000 0.2500 


DISTANCE MATRIX IS 
16.0000 6.6667 3.8095 

6.1971 7.5746 6.7750 

4.4444 8.0000 18.1818 

at ac as; s K 3s X as s 3E 3x s: s ss r; s ae SR s; 3s c SE s s s B a: ss s a B % as s = s; =s s: s 'B a X 3S as as K s s ae K s x « s sc e s 

ch-arg^ vector is 
0.01074476 
0.12595306 
-0.00304584 


as X ce at Z S S ae a S a S X S S S S B S S B SS S B S B SS B B SS SE S B se sa S B B B B B SS SB as B se SK B as B X B B B SB SB X 

CENTRE 6 RADIUS ARE 
CENTRE* 0.2000 RADIUS* 0.1000 


THETA* OPJ* 0 . OOOOOOOOOOE+OOQJ* 1 . OOOOOOOOOOE-Ot 

VOLTAGE AT POINT PJ,QJ IS 1.0000 
ERROR AT PJ,QJ* -0.0000 


THETA* 15PJ* 2.5885183382E-02QJ* 1 . 03408296 OOE-0 1 

VOLTAGE AT POINT PJ,QJ IS 0.9935 
ERROR AT PJ,QJ= 0.0065 



THETA* 30PJ* 5. 00058794E4E-0EQJ* 1 . 1 34 0085437E--0 1 

VOLTAGE AT POINT PJ,QJ IS 0.9833 
ERROR AT PJ,QJ* 0.0167 


THETA* 45PJ* T . 07 1 78T8687E-0EQJ* 1 . 292965831 8E-0 1 

VOLTAGE AT POINT PJ, QJ IS 0.9802 
ERROR AT PJiQJ* 0.0198 


THETA* 60PJ* 8.6609328687E-02QJ* 1 . 500 1 1 7S954E-0 1 

VOLTAGE AT POINT PJ,QJ IS 0.9843 
ERROR AT PJ,QJ* 0.0157 


THETA* 75PJ* 9 . 6596 9741 1 3E- 020 J* 1 . 74 1 3449030E-*0 1 

VOLTAGE AT POINT PJ,QJ IS 0.9919 
ERROR AT PJ,QJ* 0.0081 


THETA* 90PJ* 9.9999997925E-02QJ* 2 . 0002036732E-0 1 

VOLTAGE AT POINT PJ,QJ IS 1.0000 
ERROR AT PJ,QJ* -0.0000 


THETA* 105PJ* 9.65864E9876E-02QJ* 2 . £59 0485599E-0 1 

VOLTAGE AT POINT PJ,QJ IS 1.0066 
ERROR AT PJ,QJ* -0.0066 


THETA* 120PJ* 8.6588958971E-02QJ* 2 . 5002351 631 E-01 

VOLTAGE AT POINT PJ,QJ IS 1.0105 
ERROR AT PJ,QJ* -0.0105 


THETA* 135PJ* 7 . 0669 0720 1 3E- 02Q J* £ . 70T3227762E- 0 1 

VOLTAGE AT POINT PJ,QJ IS 1.0110 
ERROR AT PJ,QJ* -0.0110 


THETA* 150PJ* 4.9970599425E-02QJ* 2 . 866 1 9508 1 6E-0 1 

VOLTAGE AT POINT PJ>QJ IS 1.0080 
ERROR AT PJ,QJ= -0.0080 


THETA* 165PJ* 2 . 5845834952E-02Q J* 2 . 9660224022E-0 1 

VOLTAGE AT POINT PJ,QJ IS 1.0028 
ERROR AT PJ,QJ« -0.0028 


THETA* 180PJ*-4 . 0734639662E-05Q J* 2 . 99999991 TOE-01 

VOLTAGE AT POINT PJ/QJ IS 1.0000 
ERROR AT PJ.QJ* -0.0000 


CUMMULATIVE SQUARE ERR0R=0 . 0 0 1 3697 0 


IK « a: 3B s K 3B s SE ae ae s s S' sc = sc K B ss sc s x sc se s: = s s s s Sc s as ae sc sc s s 3C s ix K sc ir s s s s K ae «E SB SB ar SE s oc 

)»s«aB»*******ssc.quipot €nt i«l surface calculat lon**»=******» 
RI« 0.0005 21* 0. 0200 


ER« 3.0666 E2* -2.7872 


RI* 

0.0095 

21* 

0 . 0299 


ER« 

4.6862 

E2= 

-2.8177 


0.0163 

21* 

0.0413 


ER« 

6.7173 

E2» 

-2.8648 

RI* 

0.0216 

21* 

0 . 0536 


ER« 

9 .2044 

E2* 

-2.9360 

RI* 

0.0256 

21* 

0 . 0663 


ER= 

12.2973 

E2= 

-3 . 0353 



RI- 

0.0288 21* 

0. 0792 


ER* 

16.2657 

E2- 

-3. 1655 

R1- 

0.0314 ZI- 

0 . 0923 


ER« 

21 .5402 

E2- 

-3.3279 

Rl* 

0.0334 ZI- 

0. 1055 


ER« 

28.7774 

E2- 

-3.5186 

Rl* 

0.0350 ZI* 

0.1187 


ER» 

38.8662 

E2- 

-3.7186 

RI« 

0.0363 Zl» 

0 . 1320 


ER» 

52.5904 

E2- 

-3.8807 

RI« 

0.0373 ZI- 

0 . 1453 


ER* 

69.9061 

EZ- 

-3.9426 

RI* 

0.0380 ZI- 

0 . 1586 


ER* 

91 . 1 068 

EZ* 

-3.8710 

RI* 

0.0386 ZI- 

0.1719 


ER* 

118.7264 

EZ- 

-3.5263 

RI- 

0.0390 21- 

0 . 1852 


ER- 

149.6647 

EZ- 

-2.4939 

RI* 

0.0392 21* 

0 . 1 986 


ER- 

166.0374 

EZ* 

-0 . 5650 

RI- 

0.0392 21- 

0.2119 


ER* 

150. 1214 

EZ* 

1 . 4559 

RI* 

0.0391 ZI* 

0 . 2252 


ER- 

1 14 . 8871 

EZ* 

2.5972 

R2- 

0.0388 ZI- 

0 . 2386 


ER* 

81 . 3016 

EZ- 

2 . 8785 

RI* 

0.0383 ZI- 

0.2519 


EE- 

56.861 I 

EZ- 

2.7225 

RI- 

0.0377 ZI- 

0.2652 


ER* 

40.9184 

EZ- 

2.4300 

Rl* 

0.0369 ZI- 

0.2785 


ER* 

30.6781 

EZ- 

2, 1405 

RI* 

0.0360 ZI* 

0.2918 


ER- 

23.7790 

EZ- 

1 . 8903 

RI* 

0.0349 ZI- 

0.3051 


ER* 

18.8864 

EZ* 

1 . 6789 

RI* 

0.0338 Il- 

0.3184 


ER* 

ls. 2907 

EZ- 

1 .5000 

RI* 

0.0325 Zi- 

0.3317 


ER* 

ta. 5800 

EZ- 

1 .3479 

RI* 

0.0310 21- 

0.3449 


ER* 

10.4941 

EZ- 

1 .2178 

Rl* 

0.0295 21- 

0 . 3582 


ER* 

8.8605 

EZ* 

1 . 1 058 



RI* 

€R* 

0.0E78 21* 

7.5610 

0.3714 

E2* 

1 . 0087 

RI« 

0.0E6! 

21* 

0.3846 


ER» 

6.5130 

E2- 

0.9241 

RI» 

0.0E4S 

21* 

0.3978 


ER* 

5.6573 

E2* 

0 . 8500 

RI« 

0 . OEee 

21* 

0.4110 


ER« 

4.9508 

EZ* 

0.7846 

RI« 

O.OSOl 

21* 

0.4242 


ER« 

4.3615 

E2* 

0.7267 

RI* 

0.0179 

21* 

0.4373 


ER* 

3.8656 

E2* 

0.6751 

RI* 

0.0157 

21* 

0.4505 


ER* 

3.4449 

E2* 

0.6289 

RI* 

0 . 0 1.33 

21* 

0 . 4636 


ER* 

3. 0852 

E2* 

0.5874 

RI* 

0.0108 

21* 

0 . 4767 


ER* 

£.7757 

EZ* 

0.5499 

RI* 

0.0082 

zi* 

0 . 4897 


ER* 

2.5076 

EZ* 

0.5160 

RI* 

0.0055 

21* 

0 . 5028 


ER* 

2.2742 

EZ* 

0 . 4852 

RI* 

. 0 . 0027 

ZI« 

0.5158 


ER* 

2.0697 

EZ* 

0.4571 

RI* 

-0 - 0002 

ZI* 

0 . 5289 


RI* 

0.0005 

ZI* 

0.0400 


ER* 

6.5816 

EZ* 

-2.8824 

RI* 

0 . 0058 

ZI* 

0 . 0522 


ER* 

9. 1830 

EZ* 

-2.9743 

RI* 

0.0100 

ZI* 

0 . 0649 


ER* 

12.4886 

EZ* 

-3. 1 009 

RI* 

0.0132 

ZI* 

0 . 0778 


ER* 

16.8634 

EZ* 

-3.2715 

RI* 

0.0157 

ZI* 

0 . 0909 


ER* 

22.9731 

E2» 

-3.4988 

RI* 

0.0177 

21* 

0.1041 


ER* 

32.0924 

E2* 

-3.7992 

RI* 

0.0193 

21* 

0. 1 173 


ER* 

46.8322 

E2= 

-4.1878 

RI* 

0 . 0205 

21* 

0 . 1306 


ER* 

72.2681 

E2= 

-4. 634T 


RI*= 0. ©E13 ZI=s 0.1 -<139 



ER* 

1 10.8825 

E2* 

“4.9221 

RI« 

0.0E19 

ZI» 

0. 1573 


ER« 

148.6864 

E2* 

-5. 0689 

RI* 

0 . 0284 

ZI« 

0 . 1706 


ER- 

204. 7403 

E2* 

“5.6044 

Rl» 

0 . 0227 

21* 

0 . 1839 


ER« 

334 .3655 

E2« 

-5.4069 

11* 

0.0230 

21* 

0. 1972 


ER« 

475.3594 

E2* 

“1.5145 

Rl* 

0 . 0230 

21- 

0.2106 


ER- 

393.6564 

E2* 

3.9245 

RI- 

0 . 0229 

21* 

0.2239 


ER* 

227.1643 

E2« 

5.2630 

RI* 

0 . 0226 

21* 

0.2372 


ER* 

125.0433 

EZ« 

4.6779 

RI* 

0.0221 

21* 

0.2506 


ER* 

70.9499 

EZ = 

3.7832 

RI- 

0.0214 

21 = 

0.2639 


ER- 

46.4706 

EZ* 

3. 0175 

RI* 

0 . 0205 

21 = 

0 . 2772 


ER- 

34.3415 

EZ* 

2.5052 

RI* 

0.0195 

21* 

0 . 2905 


ER* 

26.3196 

E2* 

2. 1 381 

RI- 

0 . 0184 

21* 

0 . 3038 


ER- 

20.6131 

EZ* 

1 . 8540 

Rl* 

0. 0172 

21* 

0,3170 


ER* 

16.4666 

EZ* 

1 . 6266 

RI* 

0.0159 

ZI* 

0 . 3303 


ER- 

13.3921 

EZ* 

1.4411 

RI* 

0.0145 

21* 

0 . 3436 


ER* 

1 1 . 0651 

EZ* 

1 .2873 

RI- 

0.0130 

ZI* 

0.3568 


ER* 

9.2691 

E2« 

1 . 1583 

RI* 

0.0113 

21* 

0.3700 


ER* 

7.8583 

EZ* 

1 . 0487 

RI* 

0.0095 

21* 

0.3833 


ER* 

6.7324 

EZ* 

0.9549 

RI* 

0.0077 

ZI* 

0.3965 


ER* 

5.8212 

EZ* 

0.8737 

RI* 

0.0057 

ZI* 

0.4097 


ER* 

5.0745 

EZ* 

0.8030 

RI* 

0.0036 

21* 

0.4228 


ER* 

4.4558 

EZ* 

0 . 7408 

RI* 

0.0014 

21 = 

0.4360 





ER« 

3.9379 EZ* 

0.6860 

RI» 

-0.0009 21* 0.4^9t 


RI« 

0.0005 21* 0.0600 


ER» 

11.1661 EZ* 

-3.0576 

RI* 

ER- 

0.0040 21* 0.07E9 

15. £196 E2* 

-3. £234 

RI* 

ER* 

0.0068 21* 0.0859 

£0.8187 E2« 

-3.4481 

RI- 

ER* 

0.0090 21* 0.0991 

£9.1868 E2* 

-3.7554 

RI* 

ER* 

0.0107 21* 0.I1E3 

43.0886 E2* 

-4. 1 844 

RI* 

ER* 

0.0119 21* 0.1256 

70.3377 E2* 

-4.7943 

RI* 

ER* 

0.01E9 21* 0.1389 

135.8978 E2* 

-5.5199 

RI* 

ER* 

0.0134 ZI* 0.15E2 
£15.9730 E2* 

-5.3392 

RI* 

ER* 

0.0137 ZI* 0.1655 
£30.0174 EZ* 

-6 . 1 994 

RI* 

ER* 

0.0141 21* 0.1788 

407.6140 E2» 

-8 . £807 

RI* 

ER* 

0.0144 ZI* 0.192S 
949.4174 EZ* 

-7 .5129 

RI* 

ER* 

O.OT45 ZI* 0.E055 
1054.5332 E2* 

5. 6085 

RI* 

ER* 

0.0144 21* 0.ST88 

446.0E81 EZ* 

8.2598 

RI* 

ER* 

0.0141 21* 0.£3££ 

193.9043 EZ* 

6.4009 

RI* 

ER* 

0.0137 21* 0.E455 

83.254S EZ* 

4 . 6225 

RI* 

ER* 

0.0TE9 21* 0.2588 

45.0874 E2« 

3 . 4728 

RI- 

ER* 

0.0119 ZI* 0.27S1 
37.7078 E2« 

£.7895 

RI* 

ER* 

0.0109 21* 0.2854 

29.5534 E2« 

£.3463 

RI* 

ER* 

0.0099 21* 0.2987 

£3.0106 E2* 

2. 0094 

RI* 

ER* 

0.008? ZI* 0.31E0 
18 . 1 91 6 EZ* 

1 . 7446 



ER“ 


3.9379 EZ= 


0.6860 


RI* -0.0009 21* 0.^^91 


RI* 

0.0005 ZI* 0.0600 



ER« 

11.1861 E2» 

-3.0576 


RI* 

ER* 

0.0040 ZI* 0.0729 
15-2196 E2» 

-3.2234 


RI* 

ER* 

0.0068 21* 0.0859 

20.8187 E2» 

-3.4481 


RI* 

ER* 

0.0090 ZI* 0.0991 
29.1868 E2* 

-3.7554 


RI* 

ER* 

0.0107 ZI* 0.1123 
43.0886 EZ* 

-4 . 1844 


RI* 

ER* 

0.0T19 ZI« 0.1256 
70.3377 EZ“ 

-4.7943 


RI* 

ER* 

0.0129 ZI* 0.1389 
135.6978 EZ* 

-5.5199 


RI* 

ER* 

0.0134 21= 0.1522 

215.9730 EZ* 

-5.3392 


RI* 

ER* 

0.0137 21* 0.1655 

230.0174 EZ* 

-6 . 1 994 


RI* 

ER* 

0.0141 ZI* 0. 1788 
407.6140 EZ* 

-8 . 2807 


RI* 

ER* 

0,0144 ZI* 0.1922 
949.4174 EZ* 

-7 .5129 


RI* 

ER* 

0.0145 ZI* 0.2055 
1054.5332 EZ* 

5. 6085 


RI* 

ER* 

0.0144 21* 0.2188 

446.0281 EZ* 

8.2596 


RI* 

ER* 

0.0141 ZI* 0.2322 
193.9043 EZ* 

6.4009 


RI* 

ER* 

0.0137 21* 0.2455 

83.2542 EZ* 

4 . 8225 


RI* 

ER* 

0.0129 21* 0.2588 

45.0874 EZ* 

3.4728 


RI* 

ER* 

0.0119 ZI* 0.2721 
37.7078 EZ* 

2.7895 


RI* 

ER* 

0.0109 ZI* 0.2854 
29.5534 E2* 

2.3463 


RI* 

ER* 

0.0099 ZI* 0.2987 
23.0106 E2* 

2. 0094 


RI* 

ER* 

0.0087 21* 0.3120 

18. 1916 EZ* 

1 . 7446 




RI* 

ER* 

0.0074 ZI= 0. 
14.6471 EZ* 

3252 

1 .5324 

RI* 

0 . 0060 

ZI* 0 . 

3385 


ER* 

11.9961 EZ« 


1 . 3594 

RI* 

0 . 0045 

21* 0. 

3517 

- 

ER* 

9.9739 EZ* 


1.2160 

RI* 

0 . OOE9 

ZI* 0. 

3650 


ER* 

d.4019 EZ* 


1 . 0957 

RI* 

0. 001E 

ZI* 0 . 

3782 


ER* 

7.1588 EZ* 


0 . 9935 

RI* 

-0 . 0006 

ZI* 0 . 

3914 


RI* 

0.0005 

N 

II 

o 

• 

0800 


ER* 

18.1044 EZ* 


-3.3450 

RI* 

0 . 0029 

ZI* 0 . 

0931 


ER* 

25.1636 EZ* 


-3.6229 

RI* 

0 . 0048 

ZI* 0 . 

1 063 


ER* 

36.4853 EZ* 


-4.0141 

RI* 

0 . 0063 

ZI* 0 . 

1 1 96 


ER* 

57.8708 EZ* 


-4 . 5945 

RI* 

0 , 0073 

ZI* 0 . 

1329 


ER* 

114.0555 EZ* 


-5.5314 

RI* 

0 . 0080 

ZI* 0 . 

1462 


ER* 

356. -<1141 EZ* 


-6 . 3551 

RI* 

0 . 005C 

21* 0 . 

1595 


ER* 

E80.7800 EZ* 


-5 . 3687 

RI* 

0 , 0085 

21* 0 . 

1728 


ER* 

344.E901 EZ* 


-8. 2726 

RI* 

0 . 0088 

ZI* 0 . 

1862 


ER* 

949.3098 EZ* 

-13. 0141 

RI* 

0 . 0090 

21* 0 . 

1 995 


ER* 

3120.2399 E2= 


-1 .7146 

RI* 

0.0090 

21* 0. 

2128 


ER* 

1026.5351 EZ* 


13. 0084 

RI* 

0.0088 

N 

It 

o 

2262 


ER* 

323.2616 EZ* 


8.5153 

RI* 

0. 0085 

ZI* 0. 

2395 


ER* 

122.2365 E2* 


6 . 0728 

RI* 

0. 0078 

21* 0. 

2528 


ER* 

0.5017 EZ* 


4.0303 

RI* 

*0 . 0054 

21= 0. 

2545 


R2* 

0.0005 

21* 0. 

1 000 




El?« 

30.421 1 

E2» 

-3.8164 

RI* 

O.OOE2 ZI- 

0.1132 


ER« 

46. 1193 

E2= 

-4.3083 

RI* 

0.0034 21= 

0. 1265 


ER- 

81 .6683 

E2« 

-5.0998 

RI* 

0-0042 ZI* 

0. 1398 


ER* 

242.9013 

EZ* 

-6.7084 

RI* 

0.0046 21= 

0.1531 


ER* 

803.6458 

E2= . 

-3.8545 

RI* 

0.0047 ZI* 

0. 1665 


ER* 

290.2484 

E2= 

-6.8313 

RI* 

0.0050 21= 

0 . 1798 


ER* 

602 . 6244 

EZ* 

-1 1 . 6852 

RI* 

0.0052 21* 

0.1931 


ER* 

3387. 1527 

E2* -23.2938 

RI* 

0.0053 21* 

0.2065 


ER* 

3590. 1042 

E2* 

23. 0556 

RI* 

0.0052 21= 

0.2198 


ER* 

596 . 789E 

E2* 

1 1 . 7460 

RI* 

0.0050 ZI* 

0 . 2331 


ER* 

E06 . 3640 

EZ* 

7.4001 

RI* 

0.0045 21= 

0 . 2465 


ER* 

-69 . 0556 

EZ* 

5 . 6461 

RI* 

0.0056 21= 

0 . 2597 


ER* 

22. 1685 

EZ* 

3.3163 

RI* 

0 .0036 ZI* 

0 . 2729 


ER* 

36 . 1367 

E2« 

2 . 7892 

RI* 

0.0026 21= 

0.2862 


ER* 

29.2333 

E2= 

2.3533 

RI* 

0.0015 21= 

0 . 2995 



22.8136 

EZ* 

2.0105 



PROGRAM RODPLANE^ 

USES CRT; 

CONST N«10; 

typ« 

MATRIXt*ARRAY n..N,1..N3 OF REAL; 

MATRIXS*ARRAY C1..N,1..13 OF REAL; 

{fllXXKKSeXS&SSSacXXISSSSSSSS; KKXXXXXBXBXSSSSXSSXSSXXIBSXSXXXBXBXXSaBXBSXXBBWKXXiX} 

PROCEDURE INVERSE<A:HATRIX1 ; VAR B'MATRIXE); 

CONST NH-50; 


N»1 0; 

H*1 ; 

VAR 

INDXR.ARRAYCI . .NH3 OF INTEGER; 
INDXC:ARRAYC1 . .Nn3 OF INTEGER; 

IPIV: ARRAYCl . .NM3 OF INTEGER; 
BIG,DUM.PIVINV:REAL; 

I,J,V,U,LL,K,L,P,Q, I COL, IROW : INTEGER; 
BEGIN COF PROCEDURE) 

FOR P:-l TO N DO 
BEGIN 

FOR Q:*1 TO M DO 
BEGIN 

BEP,Q3 . 0; 

END; 

END; 

FOR j:*1 TO N DO 
BEGIN 

IPIVCJ3 :*0; 

END ; 

FOR I:«1 TO N DO 
BEGIN 
BIG : -0.0; 

FOR J:*1 TO N DO 
BEGIN 

IF ( 1PIVCJ3 <>1 ) THEN 
BEGIN 

FOR K:*T TO N DO 
BEGIN 

IF (IPIVCK3*0) THEN 
BEGIN 

IF (ABSCAEJ,K3 J>«BIG) THEN 
BEGIN 

BIG:«AB5CACJ,K3 ); 

IROW: «J; 

I COL : »K 


END; END 
ELSE 

IF (iriVCK3>l } THEN 
BEGIN 

WRITELNt ' SINGULAR POINT'); 
READLN 
END; 

END 

END 

END; 

IPIVCICOL3 :«IPIVCICOL3+1 ; 

IF (IROW <> ICOL) THEN 
BEGIN 

FOR L:»l TO N DO 
BEGIN 

DUH:*ACIROW,L3 ; 

ACIROW,L3 .“ACICOL^La ; 

AnCOL ,L3 :=DUH 
END; 

FOR L:=1 TO M DO 
BEGIN 

'DUM:=BEZROW,L3 ; 

BCIROy,L3 :*B£ICOL>L3 ; 
BriCOLAI :=DUM 
END; 

END; 

INDXRCI3 :«IROW; 

INDXR 113 :=ICOL; 

IF (ACICOL, ICOL3 *0 . 0 ) THEN 
BEGIN 



WRITELNt ' IT IS A SINGULAR HATRIX *♦♦♦♦♦♦*') ; 

READLN 

END; 

PIVINV:*( 1 . O/AEICOL, ICOL3 ) ; 

ACICOL, ICOL3 :*1 . 0; 

FOf? L:«1 TO N DO 
BEGIN 

ACICOL^LI :«AEICOL,L3*PIVINV 
END; 

FOR L:»l TO M DO 
BEGIN 

BCICOL,L3 :*BEICOL,L3*PIVINV 
END; 

FOR LL:»T TO N DO 
BEGIN 

IF (LLOICOL) THEN 
BEGIN 

DyH:*=ACLL,ICOL3 ; 

ACLL, ICOL3 ;«0 . 0; 

FOR L:«l TO N DO 
BEGIN 

A CLL , L3 : *A CLL , L3 -A E I COL , L3 *DUri ; 

END; 

FOR L:*l TO n DO 
BEGIN 

BELL,L3 :»B ELL,L3-BEICOL,L3*DUM 
END; 

END; 

END; 

END, 

FOR L:*N DOWNTO 1 DO 
BEGIN 

IF tlNDXR CL3 <> INDXC CL3 ) THEN 
BEGIN 

FOR K:*T TO N DO 
BEGIN 

DUN; SACK , INDXRCL33 ; 

ACK, INDXRCL33 :»ACK, INDXCCL33 ; 

AEK> INDXCCL33 :»DUn 
END; 

END; 

END; 

WRITELNC ' x: = = = s = s=:®se5s = *ca:a:«sasss*sss*a**a*®*a*aB*a**««saEa!»»«»»«»»a»««s«*««««»««« » ) ; 
END; 

CHAIN PROGRAHME 


ttaEXsssssaKassssaeaEaesx} 

CONST Hs£0; CH>sN+6} Q»10;CQ=M-N3 

VAR 2N: ARRAY C1..N3 OF REAL; CEND POINTS OF THE LINE ELEMENTS! 

PQ: ARRAY Cl.. 11,1.. 23 OF REAL; (CONTOUR POINTS , CHECK POINTS! 

MAT:ARRAY El. .M,1..N3 OF REAL; (POTENTIAL COFFICIENTS! 

CAT:MATRIX1 ; 

CHARGE :MATR 1X2; 

VOL: ARRAY Ct..Q3 OF REAL; 

DRI ,el«c,«ror,RO,THETA,R,G,L,A1 ,D1 , D£, ERROR , D21 ,D2£,dl t ,d1 2,111 , M2 , DH, DL , R I , 21 , ER , E2 : REAL ; 
El ,E2,G1 ,G2,DL1 , DL2 , FI , F2, F3 , F4, DZ , DR : REAL ; 

0,P,BI,C,D,I,J,K,S,T,U,V,II: INTEGER; 
f il©1 .text; 

BEGIN 

•»*ignCf ilel , 'san . prn' ) ; 
r«writa( f i lei ) ; 

URITELNC *R,GM; CR*RADIUS OF SPHERICAL PART, G»DISTANCE FROM THE PLANE! 

READLN(R,G); 

WRITELNC fi lei Radius of spherical part distance froa plat# *',G:8:4); 

WRITELNC filet , ‘ = = = = = = = =. = = = = = = = = = ss = = s = = = sr==.=:= * ); 

WRITELNC 'L,A1 ^ ; (L*LENGTH OF FIRST LINE CHARGE, A1 -MULTIPLYING FACTOR FOR LENGTH! 



READLNCL,A1 ); 

yRITELN( fi lei , 'Length of first line charge multyplylng factor «',A1:e:A) 

URITELN< f i lei , ' * = = ' ); 

ZHtU :*R + G; 

WRITELNt filet iZNUl :8:4) ; 

ZNt23 :-R+G+L; 
yRlTELN(filel ,2NE23 :8:4) ; 

FOR 1:*3 TO N DO 
BEGIN 

; 

K:*I-E; 

2NCI3 :*A1 ZNCJ3-2NCK3 )+2NCJ3 ; 
yRITELN(f ilel ,2NtI3 :8:4) 

END; 

yRITELNC f i lei , '*«**«*===*==«***«»***=«*«***»**« ' ) ; 
yRITELNCf ilel ); 

FOR B1:»1 TO 6 DO 
BEGIN 

THETA :*(B1-1 1419/180; 

PQCB1 , 13 :=R*S IN (THETA); 

PQEBt ,23 :*G+R*C 1 “COS ( THETA ) ) ; 

yRITECf ilel ,PQCB1 , 13 : 8 : 4 , PQ CB1 , 23 : 8 : 4) ; 

yRITELNCf ilel ) 

END; 

FOR C:*7 TO (N+6) DO 
BEGIN 

PQCC, 13 :«R; 

PQ CC, 23 : » ( ( ZN CC“63 +R + G 1/2 ) ; 

yRITECf Ilel , PQCC, 13 :8:4,PQCC,23 : 8 : 4 ) ; 

yRITELNCf ilel ) 

END; 

FOR D:eCN*f7) TO h DO 
BEGIN 

THETA :*C3 . 1 41 9/2) *C ( D“N“7 )/( M“N“6 ) ) ; 
wrltelncfilel , '•theta«',theta:8:4)i 

PQCD, 13 :«R*SIN(THETA) ; 

PQ CD, 23 :»G+R*( 1 -COS ( THETA > ) ; 
yRITE(f 1 lei , PQCD, 1 3 : 8 : 4, PQCD,£3 : 8 : 4) ; 
yRITELNCf ilel ) 

END; 

yRlTELNC f 1 lei , ' ***»**a*a*=e«aa*a«aa»* ' ); 

FOR U'.*1 TO M DO 
BEGIN 
8 : «U ; 

«1 :*C1/SQRT{SQRCPQC5,13 ) +SQR C PQ CS , 23 -ZN Cl 3 ))); 

H2:»C 1/SQRT(SC5R(PQCS, 13 > +SQR ( PQ CS , 23 +ZN E 1 3 ) )); 

«ATCU,13 :«(m-f12); 

FOR V:«2 TO N DO 
BEGIN 

D1 •■»1/(ZNCV3-ZNCV-13 ) ; 
if pq Cs , 1 3 =0 . 0 then 
begin 

dll ;«(2nEv3-pqCs,23 ) ; 
d12:*(2nCv-13-pqEs,£3 ); 
d21:«(2nCv-13+pqCs,23); 
d£2 : *( zn Cv3 +pqls, £3 ) ; 

®at Cu,v3 ;*d1 »in( (dl 1/d1£)*(d£1/d22) ) ; 

end 

else 

begin 

dll :*(ZNCV3-PQCS,£3+SQRT(SQR<ZNCV3 -PQCS,23 )+SQR(PQCS, 13 1 ) ) ; 
d12;«(2NCV-13-PQCS,23+SQRT(SQR(ZNCV-13-PQCS,e3 )+SQRCPQCS, 13 > ) ) ; 
d21 :*C2NCV-13 +PG[S,23+SQRT(SQRC2NCV-13-<-PQES,23 )+SQR(PQES, 13 ) ) ) ; 
d2E;a{2N!:V3+PQCS,£3 +SQRTCSQR(2NCV3 tPQES,23 )+SQRCPQ CS, 13 ) ) ) ; 

HATCU, V3 : «=d1 ’felnC ( dll /d1 2 )*( d21 /d22 ) 5 ; 

END ; 


writ «( f Hel , mat Eu, v3 : 6 : 4 ) 

®nd; 

wr i t « 1 n < f i 1 ® 1 ) 

»T>d; 

FOR TO N DO 

BEGIN 

FOR V:*l TO N DO 
BEGIN 

CATCU, V3 :*MATCU, V3 
END; 

END; 

INVERSE ( CAT , CHARGE J ; 

WritelnCfilel, * ** = 3 ==:=: = a = s: = =:=: = = = a*: = « = » = sr* = »* = =i = **a*s=««»:s«a««*:*«m:*s*B*«««*w«:«««««*« 
writ«ln( fi i®1 Charge vector is'); 
writ eln( f i le 1 , ' »*======:=*==»***«* ' ) ; 

wrltelnC f i lei ) ; 

For u:«1 to n do 
begin 

writeln<f ilet , CHARGE CU, 13:8:8); 
end; 

writ el nt f i lei , ' =:=s = = = = = = = = =:===ss = a: = = = = = ==a:s« = a:**:s*:=:®sa* = a:3: = a*ssKSEKS»***««: ' ) ; 

CERROR CALCULATION 

aaaaaasssssssasaaass} 

wrltelntfilel , 'ERROR CALCULATION'); 

ERROR *.«=0 . 0; 

FOR I :*CN+1 ) TO H DO 
BEGIN 
n :«I-N; 

VOL Cl 13 :«0. 0; 

FOR J:*1 TO N DO 
BEGIN 

VOL Cl 13 . “VOL Cn a+riATCI , J3 *CHARGECJ, 13 
END ; 

WRITELNCrilel , 'voltage*' ,VOLCn3 :S:8); 

ERROR : «ERR0R + SQR ( 1 -VOL C 1 1 3 ) ; 
eror : *1 -vol Ci l 3 ; 

wrlteln(file1 , 'error*', eror : 8 : 8 ) 

END; 

wrlteln( f i lei , ' cummulat i ve square error* ' , error : 8 : 8 ) ; 

{CALCULATION OF EQUIPOTENTI AL SURFACES) 

wrltelnCfllel , 'CALCULATION OF EQUIPOTENTIAL SURFACES'); 

DH:*G/5; 

DRI:»R/5; 

FOR 0:«1 TO 5 DO 

BEGIN 

RI :*0. 0; 

ZI:»0*DH; 

WRITELN(f ilel , 'RI*' ,RI, ' 21*' ,ZI) ; 

WRITELNC 'Rl*',RI,' ZI«',ZI); 

WHILE 2I<(G+R) DO 
BEGIN 

El :*EXP( 1 .5«frLN(SQRCRI)+SQR(ZI-2NC13 ) )); 

E£:»EXP( 1 .5»LN(SQR<RI )+SQR ( 2 I+ZN Ct 3 ) ) ) ; 

ER :*CHARGEC1 , 1 3 *R I ( ( 1 . 0/El )-( 1 . 0/E2) ) ; 

EZ: “CHARGE Cl , 1 3 * C 2I-ZN Ct 3 )♦( (1 . 0/El )-( 1 . 0/E2) ) ; 

FOR P:* 2 TO N DO 
BEGIN 

Gt :«SQRT(SQR(RI>+SQR<ZNCP3-2I) ) ; 

G2 : -SORT t SQR ( R I ) +SQR C 2N CP- 1 3 +Z I ) ) ; 

DLt :*SQRTCSQRCRI )+S0R(2NCP-13-2I ) ) ; 

DL2:*SQRTCSQR(RI ) +SQR ( 2N CP3 +2 I ) ) ; 

IF RI*0. 0 THEN 

BEGIN 

ER:*0. 0 

END 

ELSE 


BE(^IN 

FI :»(ZIS(£P3“ZI)/«RI*Gt ); 

F£:«(ZN£F-13-2I >/(RI>i'DL1 >; 

F3:«<2NCP-n+2I )/(RI*G£) ; 

F4:«(ZNCP3+ZI )/ (RI*DLE); 

ER:*ER-( (CHARGEEP. 13 3 / C ZN EP3 -ZN EP- 1 3 ) 3 * < F 1 ~FE+F3“F4 3 
END; 

EZ:«E2“( { CHARGE CP , 1 3 ) / ( ZN CP3 -ZNCP-n ) )♦< (l/Gl )“( 1/DL1 3~( 1/G£} + ( 1/DLE) ) 
END;COF ER & EZ CALCULATION} 

WRITELN( 'ER=' ,ER:8:4, ' EZ«^ ,E2:8:4) ; 

WRITELN( *ER*' ,ER :8 :4, ' EZ**' ,E2:8:4) ; 

WRITELNj 

CIF ER«0 . 0 THEN 

BEGIN 

DZ:*0.0; 

DR:«DRI 

END 

ELSE 

BEGIN} 

DZ:»-<ER/E2)*DRI ; 

CEND;} 

21 :*2I+D2; 

RI:*RI+DRI; 

WRITELN(f i l€l ,'RI«',R1,' 2I«',2I3; 

WRITELN( 'RI=SRI,' ZI»',ZI); 


writeln ( ' ' ) ; 

writeln( ' ' ) 


END; CDF WHILE LOOP) 

WRITELN C ' ssa = = = = = = *1*85 = as ass xa:Kssa:35 = 5Sss = *asassm***«K»*»asssB:»*««aK* «*«***«* * J ; 
WRITELNCf ilel ) 

END; COF 0 LOOP} 

wrltelnc '==*=»*=======*=*»==»■*«=»***==*=*«*«***»**=«»«******===»»' > ; 

CEND OF EQUIPOTENTIAL SURFACE CALCULAATION 

aBXSeasssBssxxsssiXSSssxBsasxssXBXsssssxx&axsxxKa]) 

wrlt#lnC'er\d of execution'); 
close? fi let ) 

END . 




program rogc ; 

USES CRT; 

CONST N»7; H>=1 ; 

TYPE 

HATRIXt-ARRAY C1..N,1..N3 OF REAL; 

HATRIXE*ARRAY C1..N,1..«a OF REAL; 

HATRIX3«ARRAY C1..N,1..E3 OF REAL; 

VAR XY: ARRAYC1 . .7/1 . .73 OF REAL; 

PQiARRAYH . .8, 1 . .23 OF REAL; 

HAT.MATRIXl; 

CHARCE:«ATRIXE; 

B:WATRIX2; 

RS;«ATRIX3; 

RAT:HATRIX1; 

V0L:«ATRIX2i 

ELEC, ©221 , ©222,ER00R,ER0R,X, DRI ,DX,P|Q, Y, A1 , A2, A1 1 , A1 2, K1 1 ,K1 2,E1 1 , El 2, RI > 21 , ER , EZ, AA1 , AA2 
ani,an2,E21 ,E22,kni ,k112,©ni ,©n 2 ,© 1 o,© 2 o,KKT ,KK£,EE1 ,EE2,K1K,K2K,WR1 ,WR2,WZ1 ,WZ2,D2,E1 
f, I, il, «1, t1, J, S, G, H,T, pi, q1,o,RU,u,v, mm: INTEGER; 

FI, F2: TEXT; 

PROCEDURE INVERSE( A-.MATRIXI ; VAR B:MATRIX2); 

CONST NH*S0; 

N«7; H«1; 

VAR IPIV: ARRAYC1 . .N«3 OF INTEGER; 

IN0XR : ARRAY El . .NH3 OF INTEGER; 

INDXCrARRAYCI . .NM3 OF INTEGER; 

BIG, DUM,PIVINV. REAL; 

I, J/V,U,LL,K,L,P,Q, ICOL, IROW: INTEGER; 

BEGIN 

€ FOR P:*1 TO N DO 
BEGIN 

FOR Q:«1 TO N DO 
BEGIN 

READCfilel ,AEP,Q1 ); 

END; 

READLN(f il©1 > 

END; 3 

FOR P:“1 TO N DO 
BEGIN 

FOR Q:=1 TO M DO 
BEGIN 


BCP,Q3:*1.0; 

END; 

EUBi 

FOR J:*t TO N DO 
BEGIN 

IPIVCJ3 :«0; 

END; 

FOR I:*t TO N DO 
BEGIN 
BIG:»0.0; 

FOR J:*l TO N DO 
BEGIN 

IF <1PIVCJ3<>1 ) THEN 
BEGIN 

FOR K:«l TO N DO 
BEGIN 

IF (IPIVEK3*0) THEN 
BEGIN 

IF ( ABS<^iCJ,K3 ) >*BI6) THEN 
BEGIN 

BIC:«ABS(ACJ,K3 ); 

IROW:*Ji 
ICOL;*K 
END; END 
ELSE 

IF <IP1VCK3>1 ) THEN 
BEGIN 

WRITELN( '♦♦♦*♦**** SINGULAR POINT ♦*♦*♦*♦*♦***'); 
writeln; 

READLN 

END; 

END 

END 

END; 

IPIVCIC0L3 :«IPIVCICOL3+1 ; 

IF (IROW <> ICOL) THEN 
BEGIN 

FOR L:«1 TO N DO 
BEGIN 

DU«:«ACIROW,L3 ; 

AEIROW,L3 :=ACICOLiL3 ; 

ACICOL.L3 :«DUH 
END; 

FOR L:«1 TO « DO 
BEGIN 

DUN. «BCIROU,L3 ; 

BCIROW/L3 :*BElCOLiL3 ; 

BC1COL.L3 :«DUM 
END; 

END; 

INDXREI3 i^IROW; 

INDXRCI3 :«ICOL; 

IF (ACICOL,ICOL3*0.0) THEN 
BEGIN 

WRITELN( IT IS A SINGULAR MATRIX 

READLN 
END; 

PIVINV:*=(1 .O/ACICOL, ICOL3 ); 

ACICOL, ICOL3 :*t .0; 

FOR L:*t TO N DO 
BEGIN 

ACICOL,L3 :*ACICOL, La »*'PIVINV 
END; 

FOR L:=1 TO M DO 
BEGIN 

BCICOL.L3 :*Bi:iCOL,L3*PIVINV 


END; 

FOR LL:«t TO N DO 
BEGIN 

IF (LLOICOL) THEN 
BEGIN 

DUM:«ACLL, ICOL2 ; 

ACLt,IC0t3 :*0.0; 

FOR L:*1 TO N 00 
BEGIN 

ACLL^Ll :«AI:LL,L3“ACIC0L,L3*DUM; 
END; 

FOR L:*1 TO H DO 
BEGIN 

BCLL,L3 :»BCLL,L3-"BCICOL,L3*DUM 
END; 

END; 

END; 

END; 

FOR L:*N DOWNTO 1 DO 
BEGIN 

IF ( INDXREL3 OINDXCCL3 ) THEN 
BEGIN 

FOR K: TO N DO 
BEGIN 

DUM.*ACK ,INDXREL33; 

ACK, INBXRCL33 r^ACK, INDXCEL33 ; 
ACK, 1NDXCCL33 :aDUM 
END; 

END; 

END; 

END; 


CEND OF PROCEDURE INVERSE > 

£a;tBiSEaBae«!:iecea:mBarssaBK'ataE:aE:eBSE:&taesEsea:r:j) 

BEGIN C main 3 
ASSIGN (FI , 'RAJ.PRN^ ) ; 

REWRITE(Ft); 

CLRSCR ; 

{^SiszcKSEssssesKssssKssBssssisEsscxsssxssxssssKxaEBaBaBKasBaesxBe} 

WRITELNCF1 , 'CONTOUR POINTS, RING RADII « RING LOCATIONS ARE') 

{[xXXXSXXSSXSSBSSSSXSXXXSSSESBSSSaSXSSBSXXSSXXXSSXXSSSSSSB} 

X:»0.5 ; 

DX:-0.50; 

PQE1 ,13 :«0.0; 

PQCt ,£3 :*0.5; 
yRITE(F1 ,PQC! ,13 :8:4); 

MRITECF1 ,PQC1 , £3 :8:4); 
yRlTE<F1,' '); 

XYCI ,13 :*0.S; 

XYC1 ,23 :e1 .250; 
yftlTECFI ,XYC1 ,13:8.4); 
yRITE(F1,' 

yRITE(F1 ,XYC1 ,23 8:4); 
yRITELN(F1 ) ; 

PQC2, 13 :*0.5; 

PQC2,23 :«0.52; 
yRITECFt,PQC2,13 :8:4); 
yRITECFI,' 

yRITECF! ,PQi:2,23 :8:4); 
yRITE(F1,' '); 
yRITE(Fl,' '); 

P:«G.5; 

FOR 1 r »2 TO 7 DO 
BEGIN 


X:«X+(1 .5*DX/i ) ; 

P:«P+DX; 

Y:«l .407+0 . 00£66*EXP(3.23*(X)) ; 

Q;«Y--1 .0; 

XY£I,I3:*X; 

XYCI,S3 :»Y; 

PQtJ, n :*P; 

PQCJ,e3:«Q; 

WRITE(Ft ,XYCI , 13 : 8 4 ) ; 
yRITECFI , ' ' 

WRITE(Ft,XYEI,£3 : 8 •• 4 ) ; 

URITELMCFI ) ; 

WRITE(FT ,PQEJ, 13 :8:4); 

WRITECFI,' '); 

WRITE (FI >PQCJ,e3 8 : 4 ) i 
WRITECFt,' *); 

END ; 

WRITELN(F1 ); 

WRITELNC FI , ' * = * = = = = := = = a: = » = s==s* = = **« = »=«*a»*K * ) ; 

WRITELNCFl , 'POTENTIAL COFFECIENTS *>; 

WRITELNCFI , ' «a = s = »«=s = = »:=s = = = **««» = *e*«*****»« ' ); 

FOR S:*l TO 7 DO 
BEGIN 

FOR T;*l TO 7 DO 
BEGIN 

AT :*SQRT(SQR(PQES, 13+XYCT, 13 )+SQR ( PQ CS, 23 -*XY ET, S3 ) ); 
AE:«SQRT(SQR(PQtS, 13+XYET, 13 ) +SQR ( PQ IS , 23 +XY ET , S3 ) ); 

Kt :«S*SQRT(PQCS, 1 3 »XYCT, 1 3 )/A1 ; 

K2:«S*SQRT( PQES> 1 3 *XY CT , 13 )/A2i 

IF K1«0.0 THEN 

b#gln 

E1:«3.14£/£ 

end 

ELSE 

BEGIN 

El 1 :*1+0 . HS + SQR(K1 ) + C 9 /64 ) ♦EXP ( 4*LN C K 1 ) ) ; 

ETC ■•*(25/256) ♦EXP ( 6 + LN( K1 ) ) + ( 1 225/ 1 6384 ) ♦EXP ( C*LN ( K 1 ) ) ; 

El :*C3. 142/S>*tE1 1+ E1H); 

END; 

IF KS»0.0 THEN 
BEGIN 

ES: *3. 142/2 

END 

ELSE 

BEGIN 

ESI : «1 +0.25»SQRCK2)+ (9/64) ♦EXP C4*LNCK£) >+(25/256 )*EXP(6*LNCK2) > 
ESS : « ( 1 £25/ 1 6384 ) ♦EXP ( 8*LN C KS ) ) ; 

E2:*(3 . 142/2 l^'CESI+EES) ; 

END; 

MATCS,T3 :*C2/3.142)»C(E1/A1 )-(E2/A2>); 

WRITE( FI , HATES, T3 :8:4) ; 

WRITECFI , ' ' >; 

END; 

WRITELNCFl); 

END; 

WRITELNCFl , ' ss««a = = ss»* = « = s» = ** = *«*s = = a*:J=«»» ^ ); 

CCALLING THE PROCEDURE INVERSE 

INVERSE (MAT, CHARGE ); 

WRITELNCFl , 'CHARGE VECTOR IS'); 

FOR I:*1 TO N DO 
BEGIN 

WRITELNCFl , CHARGE Cl ,13 -S :4) 

END; 

WRITELNCFl ) j 


DRI :* 0 - 1 ; 
MM:*1 ; 


{xae'SKSsasxaBsssss 


sstftssBa 


CLRSCR ; 

WRITELNC 'ERROR CALCULATION'); 

WRITELNC -i.++++++++++++++4-++++++++++-t*' ) ; 

raCt,13 .*-0.25; r*C2,13:«0.6 ; rsC3>n:«0.7; rs C4, 1 3 : «0 . 8 ; 

r»C! ,23 :«0.5 ; 

raCS, 1 3 :*0 . 9 ; rs C6 , l 3 : *1 . 01 ; r* CT , 1 3 : « 1 . 1 ; 

for 1 1 : *2 to 7 do 

begin 

r«Cii,23 : sO . 487 + 0 . 002e6*exp ( 3 . 23^rB Ci i , 13 ) 

•nd ; 

URITELNCFI , 'TEST POINTS ARE'); 

FOR 1:*1 TO N DO 
BEGIN 

FOR J:*1 TO 2 DO 
BEGIN 

yRITE(F1 ,RSCI, J3 :8:4); 


URITECFI,' ') 

END ; 

WRITELNCFI ) 

END; 

WRITELNCFI , ' ^ ' ) ; 


WRITELNCFI , 'POTENTIAL COFFICIENTS FOF TEST POINTS ARE'); 

for SI :*1 to 7 do 

begin 

f or Tt *. *1 TO 7 do 
begin 

All t :«SQRT(SQR(RS CSt , t3+XYCT1 , 13 )+SQR(RSCS1 ,23-XYCTI , 23 ) ) ; 
At12:«SQRT(SQRCRSCS1 ,13tXYCTl , 13 )+SQR C RS CS1 , 23 +XY CT1 ,23 ) ); 

K1 1 1 :a2*SQRT(RSLSt , n*XYCTl , 13 )/A 1 1 1 ; 

K1 t2:*2>»'SQRTCRSCS1 , 13*XYCT1 , 13 )/A1 12; 

IF Kll 1»0 . 0 THEN 
BEGIN 

ElO :*=3 . 142/2 
END 
ELSE 
BEGIN 

El 1 I :»1+0 . 25*5QR(K1 1 1 ) + ( 9/64 ) *EXP ( 4*LN < K I 11)); 

El 12:«C25/£S6)«=EXP{6*LNCK1 1 1 ) ) + ( 1 225/1 6384 ) ^EXP C 8*LN ( K1 11)); 

Elo:«C3. 142/£)*CE1 1 It E1T2); 

END; 

IF K1 12=0.0 THEN 
BEGIN 

E2o:*3. 142/2 

END 

ELSE 

BEGIN 

£221 :»1+0.25*SQR(K1 1 2) + ( 9/64 ) *EXP ( 4*LN (K1 12) )+(25/256)*EXP( 6*LN(Kn2) ) ; 
E222:*( 1225/1 6384 )*EXP(8»LNCK1 12) ) ; 

E2o : » ( 3 . 1 42/2 ) * ( E22 1 +£222 ) ; 

END; 

RATCSl ,T13 : =C2/3 . 1 42 ) ♦ ( ( E 1 o/A 1 1 1 )“(E2o/A1 12) ) ; 

WRITECFt ,RATCSt ,Tt3 .'8:4) ; 

END ; 

WRITELNCFI ); 

END; 

WRITELNCFI) ; 

CLRSCR; 

CHULTIPLICATION OF RAT WITH CHARGE AND ERROR CALCULATION 

EROOR;=0.0; 

FOR P1:=1 TO 7 DO 



3 :» 0 . 0 ; 

* 1 TO 7 DO 


3 :«VOLCP1 ,13+MATCPl ,Q13*CHARGECQ1 ,13 

(FI, ^VOLTAGE AT THE TEST POI NT* S VOL CP 1 , 1 3 : 4 : 1 E ) ; 

I-VOLCP1 , 13 ) ; 

EROOR+SQR( 1-VOLEP1 , 13 ); 

(FI, 'ERROR AT TH TEST POINT* ' , EROR • 4 : 1 £ ) ; . 

^ * * ) 

(FI); 

(FI , ' CUnnULATIVE SQUARE ERROR » ' , EROOR : 4 : 6 3 ; 

(FI , ' equipotent ial surface calculation*); 

;(F1 , ' ss*;as*s*«==s;i=;s = 5= = = = = a=*as = s = *»*s= = sc»** = * = = * * ) ; 

*1 TO 5 DO ' 

0; 

*0.1; 

.N( 'RI«' ,RI -.8:3, ' ZI*MI :8:3) ; 

.N(Ft , 'RI«' ,RI :8 :3, * 2 1 * * , Z I : 8 : 3 ) ; 

(2I<0.96) DO 

:*! TO 7 DO 

AA1 :»5QRT(SQR(RI+XY EH, 13 ) +5aR ( Z I -X Y EH , £3 ) ) ; 

AA£:»SQRT(SQRCRI+XY CH, 13 ) +SQR ( 2 I+XY CH , £3 ) ) ; 

PI :*SQRT(SQR{RI-XYCH, 13 ) +SQR ( 2 I -XY EH , 23 ) ); 

B£; *SQRTCSQRCRI~XY CH, 13 ) +SQR C 2 I +XY CH , 23 ) ) ; 

KKl :«2*SQRT(RI*XY CH , 1 3 J/AAl ; 

KK£:*2«'SQRT(RI«XY CH , 1 3 )/AA2; 

IF KK1»0. 0 THEN 
BEGIN 

EEl :*3. 142/E 

END 

ELSE 

C3. 142/E) *{ 1+( !/4)*SQRCKK1 ) + ( 9/64 ) ♦EXP ( 4*LN C KK 1 ) ) + ( E5/S56 ) ♦EXP ( 6^LN ( KK 1 3 1 + ( 1 £25/ 1 6384 ) ♦EXP ( P^LN C KKl ) J 

KE*0.0 THEN 

3,14£/£ 


(3. 142/2)^C l + ( 1/4)^SQRCKK2) + ( 9/64) ♦EXP (4^LN(KKE) } + ( 25/256 ? ♦EXP ( 6 ♦LNCKKE) ) + ( 1 225/ 1 6384 7 ^EXP ( 8^LN C KK2 ) 

1*0.0 THEN 

»(3.14£/£) 


(3. 1 42/23 ♦( 1 -( 1 /4)^SQR( KKl 3 “ ( 3/64 ) ♦EXP ( 4*LN C KK 1 3 3 - ( 5/256 3 ♦EXP ( 6^LN ( KKl 7 3-( 1 75/ 1 6384 3 ♦EXP ( S^LN ( KKl 3 7 7 j 

2*0.0 THEN 

*(3.142/23 


= C3 . 142/2)* C 1 ~( 1 /4 )*SQRCKKE)~(3/64)*EXP( 4^LN(KK2 3 3 ~ ( 5/256 3 ♦EXP ( 6^LN ( KK2 3 3“( 175/1 6384 3 ♦EXP ( e*LN ( KK2 3 3 } 


»(8QR(XYtH, 13 )*SQR< RI )+SQR(2I '-XYCH,£3 ) }*K1Kj 
*CSQR(XYCH, 1 3 1-SQR( RI ) + SQR(2I +XY CH> £3 7 )*K£K; 

*0.0 THEN 


N 

ER- ( CHARGE CH, 13/(3. 14£*RI >)*( C CWR1 -SQRCBT )*EE1 )/C AAl ♦SQR (B1 ) ) )-( { WR2-SQR ( B£ ) ♦£££ ) /( AAE^SQR CB£ ) > ) ) 

*( (2I“XYCH,£3 >#K1K)/(AA1*SQR(B1 ) ) ; 

«((2I+XTCH,23 )*K£K?/( (AA£*S(3R(B2) ) ); 

EZ-f ( CHARGE £H , 1 3 * C £/3 . 1 42 ) ♦ ( WZ 1 + W2£ ) > 

ELNC' ER» %ER:8:3,' E2* ',E2:8:3); 

‘ELN(' ER« '^ER:8:3,' E2* ',E2:8:3); 

::*SQRTCSQR CER>+SQR(E2) ) ; 

fELNC 'ELECTRIC FIELD * ' , ELEC : 8 • 4 ) J 

fELNC 'ELECTRIC FIELD « ' , ELEC : 8 * 4 ) ; 


*-(ER/E2)*DRl ; 

*2H-D2; 

■RI+DRI; 

rELN( ' ' >; 

TELN(' RI= '.RI:8:3,' 21* ',21:8:3); 

TELN; 

) j 

TELNCFI,' RI- ',RI:8:3,' ZI* ',21:8:3); 

TELN ; 

;C0F WHILE LOOP3 

TELN ( 'as*«s:casK« = ss = asi=: = = s=s=:saass = * = a = = »s:a:**=:*a:*3:a**a=s*a*«a:K***a*«a«« ' ) ; 
TELN (Ft , ' * = »««« ®=s = a = * = = »« = = « = = 2: = = ' ) ; 


i; COF G LOOP! 

ISECFl ) ; 

TELNC'CND OF EXECUTION') 


CONTOUR POINTS, RINO RADII It RINO LOCATIONS ARE 


0.0000 

0.5000 

0.5000 1 .2500 



0.5000 

0.S200 

0.8750 

1 .5319 



1 .0000 

0.5319 

1 .1250 1 

.5877 



1 .5000 

0.5877 

1.3125 1 

.6715 



£.0000 

0.6715 

1.4625 1 

.7865 



£.5000 

0.7865 

1.5875 1 

. 9355 



3.0000 

0.5355 

1.6946 2 

.1210 



3.5000 

1.1210 





isacseaBasaKiB 

* as as *■:«=:«=: a: 

ae as a: aa as s as *s as as 3B as s 




iTENTIAL 

COFFECIENTS 




cat SB IB 

sBt sac as as IB ass as 

as as aa as as as as SB as as as as as 




0.5600 

0.2871 

0.2174 0.1743 

0 .1450 

0 . 1233 

0.1 060 

0.4541 

0.2817 

0.2195 0.1784 

0 . 1 494 

0 . 1274 

0. 1 097 

0.3089 

0.2221 

0.1853 0.1577 

0.1364 

0.1192 

0.1 046 

0.1897 

0.1648 

0.1469 0.1317 

0.1190 

0. 1080 

0.0979 

0.1227 

0 .1209 

0.1134 0.1063 

0 . 1000 

0 . 0941 

0.0884 

0 . 0855 

0.0913 

0.0888 0.0861 

0.0836 

0 . 0813 

0.0787 

0.0637 

0.0715 

0.0713 0.0709 

0.0706 

0 . 0704 

0.0700 


CHARGE VECTOR IS 
- 54.3635 
-lefS. 1773 
15857.6265 
-67115.5745 
55338.7670 
-61866 .7254 
1521 1 .5073 

TEST POINTS 
0.2500 
0 . 6000 
0.7000 
0 . 8000 
0.5000 
1 .0100 
1 . 1000 


ARE 

0.5000 
0.5055 
0.5125 
0 . 5222 
0 . 5357 
0.5565 
0.5755 


POTENTIAL COFPICIENTS POP TEST POINTS ARE 


0.5357 

0.2841 

0.2167 

0. 1743 

0 . 1450 

0 - 1233 

0 . 1 060 

0 . 4420 

0 . 2635 

0 . 2083 

0. 1707 

0 . 1436 

0 . 1226 

0 . 1 055 

0.4055 

0.2561 

0 . 2048 

0. 1 652 

0, 1431 

0 . 1225 

0 . 1 063 

0.3776 

0,2476 

0 . 2007 

0. 1 674 

0 , 1426 

0 . 1230 

0 . 1 068 

0.3480 

0.2352 

0 . 1 566 

0 . 1 656 

0 . 1421 

0 .1234 

0 . 1 076 

0.3152 

0.2307 

0. 1526 

0. 1641 

0. 1421 

0 . 1243 

0 . 1 051 

0.E551 

0.2250 

0. 1502 

0. 1 636 

0. 1428 

0 . 1258 

0.1110 


VOLTAGE AT THE TEST POINT«0 . 95955957951 

ERROR AT TH TEST POINT*0 . 00000002045 


VOLTAGE AT THE TEST POINT»1 . 00000000370 
ERROR AT TH TEST P01NT*-0 . 000000003T3 


VOLTAGE AT THE TEST P01NT»0 . 99955597765 
ERROR AT TH TEST PO I NT*0 . 0 0 0 0 0 0 02235 


VOLTAGE AT THE TEST POINT*0 . 95955599255 
ERROR AT TH TEST POINT*0 . OOOOOOOOT45 


VOLTAGE AT THE TEST POINT*! . 000000011 20 
ERROR AT TH TEST POINT®-0 . 00000001 1 1 8 


VOLTAGE AT THE TEST PO INT* 1 . 0 0 00 00 0 0 0 0 0 
ERROR AT TH TEST POINT=0 . 00000000000 


VOLTAGE AT THE TEST POINTaO . 5999999981 4 


ERROR AT TH TEST PO I NT*0 . 000 00 0 0 0 1 86 


CUMHULATIVE SQUARE ERROR «0. 00000000 
9qylpotent lal surface calculation 


RI» 

0. 000 

ZI- 

0.100 

RI« 

0.100 

21* 

0.100 

RI» 

O.EOO 

21* 

0.115 

RI» 

0 . 300 

21* 

0 . 147 

RI* 

0 . AOO 

21* 

0.201 

RI« 

0 . 500 

21* 

0 . 285 

RI« 

0 . 600 

21 = 

0,407 

RI» 

0 . 700 

21 = 

0.561 

Rl- 

0.800 

21* 

0.716 

RI« 

0 . 900 

21 = 

0 . 860 

RI* 

1 . 000 

21* 

1.017 

RI* 

0.000 

21 = 

0.200 

RI« 

0.100 

21 = 

0.20 0 

RI» 

0.200 

21* 

0.214 

RI* 

0.300 

21 = 

0.244 

RI* 

0.400 

21 = 

0.293 

RI- 

0.500 

21* 

0.367 

RI« 

0.600 

21 = 

0.470 

RI* 

0,700 

21* 

0.595 

RI* 

0 . 800 

21 = 

0.728 

RI* 

0 . 900 

21 = 

0.862 

Rl« 

1 . 000 

21 = 

1 . 01 7 

Rl- 

0 . 000 

21* 

0 , 300 

RI- 

0 . T 00 

21* 

0.300 

RI* 

0 . 200 

N 

If 

0.312 

RI- 

0 . 300 

21* 

0 . 338 

RI« 

0 . 400 

21* 

0.381 

RI* 

0 . 500 

21* 

0 .444 

RI- 

0 . 600 

21* 

0.527 

RI* 

0 . TOO 

21 = 

0 . 629 

RI* 

0 , 800 

21 = 

0.742 

RI* 

0 . 900 

21* 

0.867 

Rl» 

1 . 000 

21* 

1.018 

RI* 

0. 000 

21* 

0.400 

RI* 

0.100 

zi« 

0.400 

R1- 

0.200 

21 = 

0.410 

RI* 

0 . 300 

21 = 

0.431 

RI* 

0 . 400 

21* 

0.465 

RI* 

0.500 

21 = 

0.515 

RI* 

0.600 

21 = 

0.581 

RI* 

0.700 

21 = 

0.662 

RI* 

0.800 

21 = 

0 . 758 

RI* 

0.900 

21* 

0 . 872 

Rl« 

1 . 000 

21* 

1.019 

RI* 

0.000 

21* 

0.500 

RI* 

0.100 

21 = 

0.500 

RI* 

O.EOO 

21 = 

0.507 

RI* 

0.300 

21 = 

0.521 

RI* 

0.400 

ZI = 

0.545 

RI* 

0.500 

21 = 

0.581 

RI* 

0.600 

21 = 

0.630 

RI* 

0.700 

21 = 

0.693 

RI* 

0.800 

21 = 

0.774 

Rl« 

O./fOO 

21* 

0.879 

RI* 

1 . 000 

21* 

1 .020 


program general? 

1 a be 1 EOG ? SOC > 30 G y 400 ? oOG ? 

COraST CG==9.QeG9? 

NE=^3? 

NB=7? 

EB=10! 

EBB=17? 

N=17? 

TYPE MAT1=ARRAYi: 1 . .EBBy 1 . .EG OF REAL: 

MAT£=ARRAYC1 . .EB, 1 . .EG OF REAL? 

MAT3=ASRAYC1 . .EEE'y 1 - .EEEG OF REAL: 

MAT4=ARRAY[:1 . . NE , 1 . .EBEG GF REAL? 

HAT5==ARR A YC 1 .. NEG GF REAL? 

MAT6=ARRAYE 1 . . EBB . 1 » 1 G GF REAL ? 

HAT7=ARRAYl1 . .EBB, 1 . .EG OF REAL? 

VAR 

C1yrvlOM,ANSUGyUyE1yEEyD11,D1E,DE1yDEEyEUlyED£yM1 1,rr.2yr-1E1,M£2,I-1ylEyP 1 , P2 , P3 y Q 1 , Q2 : INTEGER ? 
G1 y NM y F G y AN33 r LQ y H 1 1 y GO y H 1 2 y n2 1 y H22 y Ar430 y G 1 y G2 - H 1 y HE , AtTB y X 1 , V 1 y K 1 y L 1 y M “! y N 1 .. 0 1 , R 1 y I y J y K H INTEG 
RZSHAT1? 

PQsMATE: 
mat : NAT 3 y 
THIS NAl 0 ? 

THE ! MAT 4 ? 

THETA :MAT4? 

CHARGE; MAT 6? 
sc; MATE? 

CSS MATE ? 

HHSMAT7 ? 

TAf'4 ; ARRAY £ 1 . .EEEG GF REAL ? 

ERRS MAT:? 

ERORsARRAYl 1 . .EBEG OF REAL? 

VOL s ARRAYC 1 . .EEEG OF REAL; 

CUMyTANQI yRISyZISy 2Iy A0yB0yRIE,2IE,C0,A2yB2,C;r‘yERyEZ,DRyDZ,DLyEPyRIyEF£:REAL? 

D1 yDEyRIl yRIE, ZI1 y ZIEyRIAy ZIASREAL: 

FlyFJsTEXT? 


FUNCTI0F4 F < R Z ; MAT 1 ? PQ s MATE ? F 1 , FE s I.NTEGER ) s REAL ; 

BEGIN 

F s = ( 1 /S3RT ( SQR ( R Z lF 1 , 1 G-PQCFE y 1 G ) +SGR R Z CF 1 , EG~PQCFE , EG ) ) ) 
END? 

FUNCTION FER ( R 2 s MAT 1 ? CHARGE s MAT 6 ? R I , Z I s REAL ? V s INTEGER ) s r ea 1 ? 
VAR A y y 0 ; RE-Au- ? 

BEGIN 

Af=RI-RZi:Vy 1G; 

B ! = Z I —R 2 0 V y E G ? 

C s =EXF ( 1 . G-^^-LN ( SQR ( A ) +SQR ( B ) ) ) ? 

FER s =CHARGEE V , 1 G* ( A/C ) ? 

END ? 

FUNCTION FEZ ( R2 s MAT 1 ? CHARGE s MATS ? Ri y 2 1 s REAL ? Ws INTEGER ) : r eal ? 
VAR A y E y C s R E AL ? 

‘ BEGIN 

A:=RI--RZCWy 1G? 

E;=2I“RZCWy£G? 

( 


C s =^EXF ( 1 . Si^fLN ( SQR ( A ) 

FEZ ! ==CHARGEZN y 1 G-j^- ( E/ 

Emi 

B£GIs-<-= 

FV s ^CHARlEZ GL , \ 1 / 

END ? 

F'3 ( R Z ; MAT 1 ; PQ s MATE ? SC s MAT 

VAR LsREAi-? 

BEGIN 


1-0 an. ( cl ! ; i ? 

C) 

CHARGE s MAT t 

sqrt^sqf^r; 


Ej-~ERRSG1yEG) ) ) 


FUNCTION 


2 . 


L:=EXF ( 1 .5*Lr'4(SQR (RZCU, ID-PGlV, 1 3 ) h-SQR ( R2 CU , E J-PQCV , EZ ) ) ) ? 

F3s='-( (RZCU, 1D-PQLV, 13:‘*SCLW,2j-r(RZC:U,£3-PQCV,23)->^3CCW, 1 j)/L 
END t 

FUNCTION FH(R2:riAT 1 ; HH : MAT? ? CS : NATE ? PO , LO : INTEGER ) :REAL; 

V AR LI I REAL ; 

BEGIN 

LI s =EXP < 1 . 5*LN ( SQR ( R ZCPD r 1 3—HHELO ,13) -rSQR ( R 2 CPO , £ j~HHCLO , £j ) ) ) ? 

FH s = ( ( R Z C r- 0 r 1 3”HHCLO , 13)-?i-C3ELG, 13+(RZEFG,E3“~HHELur£3)"^CSELC,£'3 )/Ll 
END? 

PROCEDURE INVERSES A:MAT3? UAR B:MAT6): 

CONST NM=^00? 

N=17? 

M= 1 ? ■ • 

NE=3 ; 

NB=7? 

EB=1C? 

VAR INDXR:ARRAYC1. .NM3 OF INTEGER ? 

INDXC : ARRAYu 1 . .NM3 OF INTEGER ? 

IPlVsARRAYC 1 . .NM3 OF INTEGER? 

BIG , P I V IN V , EiUN • real ? 

1 , J , U , 'v , LL 7 l. , K , F y Q t iCGL 7 IRON s INI EGEr\ ? 

BEG X T-J C Gr” F R GC EDUR E I 
FOR P I ~ 1 TO NE DO 
BEGIN 

Fun G ! = ‘i TO M DO 
EkEGIN 

BGF „G; 

END : 

END ? 

FOR p!=(N£+ 1 ) TO N DO 
EiEOIN 

FDF TO F DO 

BEGIN 

BCF,Q3s=0.0; 

END ? 

END ? 

FOR J:=1 TO N DO 
BEGIN 

IPIVGJ3:=0? 

EF4D 7 

FOR l:=1 TO N DO 

BEGIN 

BIG?=C.C? 

FOR J : =- 1 TO N DO 
BEGIr4 

IF (IPIVCJ301) THEN 
BEGIN 

FOR K:=1 TO N DO 
EiEGIN 

IF ( I P I V G K 3 = G } THE r4 
BEGIN 

IF ( ABS ( AG J 7 K 3 ) >=£■ IG ) THEN 
BEGIN 

BIG s =AE-S ( AG J 7 K 3 ) ? 

IRON ! ~ J J 
1 CGL ? -~r[ 



NfL„ i LLN s ’ Sii4bUw,.An POIN'; 


3 


REhDLN 
END j 
END- 
END 


IPiVClCGLD :=IFIVCICOLD-i- 1 ? 

IF ( I R D W •’ > I C Gl ) THEN 
EEuIN 

FOR Li— 1 TO N DO 
BEGIN 

DUH;=ADIROW,LD; 
ACIROUI^LDJ^-ADICOL^LD ? 
ACICOL.LD:=DUNj 
END r 

FOR TC N DO 

BEGIN 

DUN ! " B E 1 R ON ? L j ; 

BE IRON, Lj :=^BEICOL,LEt 

BEICGL,LD:=DUN 

END r 

END? 

I ND XREID !i=IR ON ; 

I ND X R E I D ! I C C w ? 

IF 'i AE I'L UL , ICGlD=G . 0 ) THEN 
BEGIN 

WRIT£Lr-4C ’h-INGULAR MATRIX* ) ? 

RE ADEN 
END; 

F I V I N V I = ( 1 . 0 /’ A E I C OL r I C GL D } ; 

A E 1 C GL r 1C OL D J = '■] . 0 ; 
r GR L J — 1 TG N DO 
BEGIN 

AE ICGL , L j : =AE ICOL , LG*P IVINV 
END ; 

FOR L!=1 TG M DO 
BEGIN 

BEICGL , LD 5 =BEICOL , LD*P IVINV 
END ? 

FOR LLi=1 TG N DO 
BEGIN 

IF (LL<>ICGL) THEN 
BEGI N 

DUM I =AELL r I COLD ; 

A ELL y I COLD i — G . G ; 

FOR L:=1 TG N DO 
BEGI N 

AELL , L3 s =AELL . LD-AE ICOL , LD^DUM ? 
END T 

FOR LJ=1 TG n DG 
BEGIN 

BELL . LD - =BELL ? L j~BE I COL , LD-^DUM 
Er-4D ? 

END? 

END ? 

ENEL' 

F Gr» L J ~N DGNtvTG 1 DG 
BEGIN 

IF ; I r4 E* X R L L r. 1 r4 D >. -C E L J . 1 H E r--i 
BEGIN 

FDR = 1 TG ri DO 
BEGIN 

BUM 2 — AEK , INDXRELD j ; 


ACK , INDXRCLD j : ^ACK , INDXCCL jD ; 

AEK,INDXCCL33s=DUM 
END .! END ? 

END? 

WRXTELN( ’ = = = = ====: = = 5=-=:===: = = :i.===s:=;======s:=::=-=;===iss;=:=s;=. ’ ) 5 

END?-CQF PROCEDURE} 

BEGINCCF NAIN PROGRAMME} 

ASSIGN ■; FI , ’MAT. DAT ’)? 

REWRITE ( FI ) ? 

ASSIGN (FJ , ’DAT.PRN’ ) ? 

RESET (FJ ) ? 

WRITELN(F1, ’GIVE THE LOCATION OF SIMULATION CHARGES INSIDE THE ELECTRODE’) 
WRITELN< ’GIVE THE LOCATION OF SIMULATION CHARGES ir^SIDE THE ELECTRODE’)? 
FOR E1 TO NE Du 
BEGIN 

FOR E£:=1 TO E DO 
BEGIN 

READ < FJ . R Z C£ 1 y EEC ) ? 

WRITE(FIyRZ[:EUEEJJS!:4y ’ ’ ) ? 

END; 

READLN(FJ ) ; 

WR ITELN (FI) ? 

END? 

WR ITELN ( F I , * GI VE THE LOCATION OF SIMULATION CHARGES INSIDE THE AIR’)? 
WRITELNC ’GIVE THE LOCATION OF SIMUuATIGN CHARGES INSIDE THE AIR* ) ? 

FOR D1i:=1 TO NE DO 
BEGIN 

FOR D IE 2 = 1 TO E DO 
BEGIN 


READ ( F' J . R Z C ( D 1 1 +NE ) , D 1 E 3 ) ; 

WRITE(FI,RZC(D M+NE) ,D1E3s£24y ’ ’ ) ; 

END: 

WR ITELN (FI) ? 

READLN ( F J ) ? 

END? 

WRITELNvFIy ’GIVE THE LOCATION OF SIMULATION CHARGES INSIDE THE DIELECTRIC’ 
WRlTELNl ’GIVE THE LOCATION OF SIMULATION CHARGES INSIDE THE DIELECTRIC ’ ) ? 
FOR DE1:=1 TO NE DO 


BEGIN 

FOR D£E:=1 TO S DO 
BEGIN 

READ ( FJ , R Z C ( DE 1 “i-EE ) . DEE 3 ) ? 
WRITECFIyRZKDEI+EB) yD223sSs4y ’ ’ ) ? 

END? 

READLNCFJ)? 

WRITELNvFI ) : 


END ? 

WRITELNC FI) ; 


WRIlELNCFIy ’GIVE THE NUMBER OF CONTOUR POINTS 
WRITELNC ’GIVE THE NUMBER OF CONTOUR POINTS GN 
READLN ( FJ y ED 1 ) ? 


ON THE AIR--E INTERFACE 
THE AIR~"£ I^4TERFACE’' ) : 


WR ITELN ( F I y ED 1 ) ? 

EDS I =NE"£D 1 y 

IF (EDI >0.0) THEN 

BEGIN 

WRITELNCFIy ’GIVE THE LOCATION OF ABOVE CONTOUR POINTS’)? 
WRITELNC ’GIVE THE LOCATION OF” ABOVE CONTOUR POINTS’ ) ? 

FOR M 1 1 ! =■' 1 TO EDI DO 
EtEGIN 

FQK n'lDU — "i TD E DU 
BEGIN 

READ i F wi y P Q 3 M 1 1 y n '1 E 3 ) ? 


- 



WRITE(FI.FQCM1 1,M123:Ss4, ’ * ) ? 

END? 

READLN(FJ ) ? 

WRITEL^4( FI) 5 
END ? 

END? 

IF (ED2>G.O) THEN 
BEGIN 

WRITELN(FI , ’GIVE THE LOCATIDtsi OF CONTOUR POINTS ON THE THE 
WRITELN ( ’GIVE THE LOCATIOr-4 OF CONTOUR POINTS ON THE THE D~E 
FOR M£1:=1 TO EDE DO 


BEGIN 

FOR N2E’. = 1 TO 2 DO 
BEGIN 


READCPQi: (ME 1-rED1 ) , MSEC ) ? 

WRITECFI ,PQ[: (ME 1+EB 1 ) .MEED :6!4, ’ ’ ) ? 

END? 


READLN? 
WRITELN(FI ) ? 


END ? 

END? 

WRITELN' FI, ’CONTOUR POINTS ON THE AIR-D INTERFACE ARE’)? 
WRITELN( ’CGt-iTuUR POINTS ON THE AIR—D INTERFACE ARE’ ) ? 

FOR I1s=1 TO NB Du 
E’iEGlN 

FOR IE:=1 TO E DO 
BEGIN 

PQEKNE+I 1 ) , IEj: = (RZi: * NE-fH ; , IE J-j-F C L i EB--1 1 ) , lEC) /E? 
WRITE. Fi,Fo!L(N£ + I1),I£CsSs4, ’ ’ )? 


END ? 


WRITELN^FI) ? 

END? 

WRITELN FI ) ? 

WRITELNC ’GIVE THE RELATIVE PERMEABILITY OF D’)? 
WRITELN ( FI ,* GIVE THE R:ELATIVE r ERMEABIlITN' OF' D’ ) ? 
READLN(EPE) ? 

WRITELN (FI, EPESSS4) ? 

WRITELN ( FI ) ? 

I cal culat ing Ihe theta niatri;-. 


FOR P15=1 TO NE DO 
BEGIN 

IF ( R 2 C ( NE-ep 1 ) , 1 3“R 2 Z ( EE-<-P 1 ) , 1 1! ) =G . 0 THEN 
TANQ1:=-.GE1C 
ELSE 

T ANQ 1 ! = ( R 2 LNE-i-P 1 , 21!— R 2 CEB+P 1 , Eli ) / ( R 2CNE+F 1 , 1 3— R 2 CEB'+P 1 , 1 
TH 1 L P '» 3 ? ==• An C T Ar4 1 ANG 1 ) ? 

END? 


FOR P1!=1 TO NB DC 
BEG 1 14 

SCLFI , 13:=ABS(SIN(TH1[:P13) ) ? 
BC CP 1 , E3 2 =AE‘S ( COS ( TH 1 CP '1 3 ) ) 
END? 


X c o 1 c u 1 a. 1 1 r* 9 t h & p o t e ri t i a. nia t r i y-. 

IF ( ED 1 >0 . G ) THEN 
BEGIN 

FOR X 1 ? TL ED 1 DO 
'BEGIN 

FDR, V ) 2 ~ 1 TO r4E DO 


D-E INTER 
INTER r AC 


m 


BEGIN 

MATE X 1 , Y 1 D : = ( CO ) ( RZ , PQ , Y 1 , X -1 ) ; 

END? 

FOR Y 1 s = ( NE-»" 1 ) TO EB DO 
BEGIN 

MATE X 1 , Y '112=0.0? 

END ? 

FOR Y 1 2 = ( EE-t- '1 ) 1 G EBE DO 
BEGIN 

MAT E X 1 , N' '1 j ! = V C G > ( R 2 ? F G , Y 1 , X ‘1 ) ; 

END? 

ENDyluF EE5 t ROWS CALCULAT IDNj' 

END ? IGF ED 1 >0 . C } 


FOR K 1 i — ( ED 1 1 ) TO NE DO 

BEGIN 

FOR L'iS = 1 TO NE DO 
BEGIN 

MATCH 1 y L -1 G : = ( CO ) -J^F ( R 2 y FQ , L '1 y K 1 } ? 
END? 

FOR L1 : = (NE-r 1 TO EE DO 
BEGIN 

MATEK 1 y L 1 G : = ( CO ) ^F ( R2 y PQ y L'1 y K -1) ? 
END ? 

FOR L1:=(EB+1) TO EBE DO 
BEGIN 

MATEK 1 ,L 1G:=0.G? 

END ? 

END? 


F'QR M 1 s = ( NE> 1 ) TO EB DO 
BEGIN 

FOR N'i: = 'i TO NE DO 
BEGIN 

MATEM1 yN1G:=0.0? 

END? 

FOR N1: = (NE+'1) TO EE DO 
BEGIN 

MAT C M '1 y N 1 G 2 =- ( C 0 ) *F < R Z y F Q y N '1 y M 1 ) ? 
END? 

FOR N12 = <EE+'1) TO EBB DG 
BEGIN 

MATEM 1 y N 1 G 2 = ( CO ) -M-F ( R2 y FQ y r4 1 , M 1 ) ? 
END? 

END? 


FOR 01 s^CEB+l ? TO EBE DG 
BEGIN 

F OR R 1 2 = 1 1 0 NE DC 
BEGIN 

MAT E G 1 y R 1 G 2 =• 'i, EF £•“ 1 ) C C G ;■ '-^-F 3 % R Z y F G , SC y R 1 . G 1 — NE y 0 1 —EB ;• 
END? 

FOR R1: = (NE'M:- TO EB DG 
BEGIN 

MATE0 1 y R 1 G 2 =EF2-C0^FS ( F 2 y FG y SC y R 1 , 0 1 -NB y 0 1-EB ) 

END? 

FOR R 1 2 =: ( EE">" ) TO EBBi DO 
BEGIN 

MAI L I » 1 . Fi 1 G 2 =— i. CD ) £.M R Z , F G . SC . R 1 .01 —NE . O 1 — EB^ ) 

END? 

END ; 


FOR l;=1 TG EBB DO 
BEGIN 

FOR Js=1 ID EB DO 
BEGir4 

WRITE ( NATH I , J j : 6 , * “ ) ; 

END? 

WRITELN; 

E^4D ? 

WRITELN ( = = = = = = = = = = 

FOR l:=1 TG EBB DO 

BEGIN 

FOR Js=(EE41) TO EBB DO 
BEGIN 

WR I TE ( NAT 1 1 j J j « S y * ’ ) 

END? 

WRITELN 


END j 

4- + -f-r--r 

1 n V e r s e ( ni* t , c ha r g e ) j- 

*r 4- 4' "r + -i- -r -i- -r T -j- -f -f -i- -?<■ -r~ 4 4' 4' 4 4 "i 4 -i- - 

WRITELNCFly 'CHARGE VECTOR ISV 

FOR C 1 s “• "i TG N DO 

BEGIN 

WRITEL^'I FI , CHARGECC 1 y s S ) ? 
WRITELN ( CHARGECC 1 , 1 C s 6 ) ? 

END ? 

C e r r 0 r c a I c a 1 a t i o n 


, _j- 4 


w r i t a 1 n F I t * u o y o u w a. f; t I- o c a 1 c u 1 a t a t l*i a e r r' or', if y « a 'then write 1 also w r" i t a G ’ / ? 
wr'itfelriC "do yoo warjt to t alculata the error, if yea theri write elae write O'; ; 

R E A D L N ( A N S 0 ) ; 

WRITELN ( FI yANSO) ; 

IF ANSC= I THEN 
BEGIN 

WR ITELN ( FI r * GI VE THE NUNBER OF SUCH TEST POINTS ON THE AIR-~E INTERF’ACE ' ) ? 

WRITELN( 'GIVE THE NUMBER OF SUCH TEST POINTS ON THE AIR~E INTERFACE') 5 
READLN(H1 ) ? 

IF H 1 y’C . C THEM 
BEG I N 

WR I TELN < FI , ' GI VE THE LOCATION OF SUCH POINTS*); 

WRITELN ; 'GIVE THE LOCATION OF SUCH POINTS'); 

FOR = 1 TG Hi DG 

BEGIN 

FOR H1Ss=1 TO £ DO 
BEGIN 

READ ( ERR CH 1 1 y H IE J ) ; 

WRITE ( F I y ER R C H 1 1 . H 1 E C ^ S i 4 ? ' * ) r 

WRITE \ ERRCH 1 "i , Pi 1£j s £ s 4 , ’ ' ) ? 

END ; 

READLN • 

WRITELN ; 

WRITELNvFI) 

END ; 

END ? 

WR ITElN ( FI y 'GIVE IHL NuME'ER GF i EBT POINTS ON Tr-:E D— E INTERFACE' ) ; 

WR 1 . wu..:'-, 'v ' v’l- i HE r-i CL i:- E R Ui" i Lb i rUli-H i d GN 1 riS D~L 1NTE^--.FACE *' ; V 
READLN \ hE ■ i ■ 

j» r r'l d w . ’w' 1 T'i c, ! ‘1 

BElIn 

f Ui'- nu 2 =•■ ■! i U riL Du 


BEGIN 

FOR H22s=1 TO 2 DO 
BEGIN 

readcerrche-uhhe:]) ? 

WRITECERREHEI ,H£EI;ai4, * '' ) ? 

NR 1T£ (Fly ERR C H2 1 ? HE2 1 s £ s 4 , * ’ ) ; 

END ? 

READLN? 

WRITELN(FI) ? 

WRITELN 

END? 

END? 

FOR Gl: = 1 TD'(H14-HE) DO 
BEGIN 

vOLCGIO :=C.O? 

END? 

IF H1>C.O THEr4 
BEGIN 

FOR TO Hi DO 

BEGIN 

FOR GE!=1 TO NE DO 
BEGIN 

VOLCG 10! =vOLlIG 1 0+CO*FV ( R 2 y ERR ? CHARGE , GE , G 1 ) 
END? 

FOR GE: = (EEt'1} TO EBB DO 
BEGIN 

VOLCG 10! =VOLCG 1 0-J- CO-^^rFV ( R 2 y ERR , CHARGE , G£ y G 1 ) 
END? 

END; 

END; 

IF HE >0.0 THEN 
BEGIN 

FOR G 1 ! = ( H 1 + 1 ) TO ( H '1 +HE )■ DO 
BEGIN 

FOR G2:=1 TO NE DO 
BEGIN 

VOL C G 1 0 : -VOL C G 1 0 + C 0* F V ( R 2 , ER R , C H AR GE , GE , G 1) 
END; 

FOR GEl=(NE+1) TO EE DO 
BEGIN 

VOLCG 10! =VOLCG 1 j-i"CG#FV ( R 2 y ERR y CHARGE , GE y G 1 ) 
END; 

END; 

END ; 

WRITELN C 'THE ERROR VEVTOR IS' ) ; 

WRITELNCFIy 'THE ERROR VEVTOR IS'); 

CUM!=0.C; 

FOR G1:=1 TO (H1+HS) DO 
BEGIN 

ERORCG10:=VGLCG10-l; 

WR 1 TELN (Fly ER OR C G 1 0 s £ ! 4 ) ; 

WR ITELN ( EROR CG 1 0 ! S ! 4 ) ; 

CUM ! ^CUM-^SQR ( ERORCG 1 0 ) 

END ; 

WRITELN ( FI y ' CUMMULATIVE SQUARE ERROR IS'); 
WRITELN ( ' CUMMUu_ATI VE SQUARE. ERROR IS* ) ; 
WRITELN(FIyCUM!£!4: ; 

WRITELN(CUM:£!4) ; 



WR ITELN ( ’*D0 YOU WAt^T TO CALCULATE THE EQUXFOTENTIAL i:>URFACL« IF YES THEN* ) j 
WRITELNC ’WRITE 1 ELSE WRITE C*): 

WRITELNvFI , *DO YOU WAf-'T TG CALCULATE THE EQUIPOTENTI AL SURrACE. IF YES THEjN*)? 
WRITELNCFI,* WRITE l ELSE WRITE 0*;? 

READLNC ANSI ) ; 

WR ITELr-4 (FI yANS 1 ) t 
WRITELN(ANS 1 ) ; 

IF AI4S1=^"1 THEN 
BEGir4 

WRITELN (FI ,* cal cu la tin 9 the equi potent la I surfa ce*)? 

•C calculate rig this ecjui potential surface 

EOOi WR XTELt4 ( FI y * GI VE THE Ir4ITIAL POITmT OF THE DESIRED EQUIPGTEr4TiAL SURFACE’)^ 
WR ITELTl ( * G1 VE T'HE ITmITIAL POINT OF THE DESIRED EQUIPGTENT I AL SURFACE* ); 
READLN(RISy 21S) ; 

WR I TE Lhi ( F I , R I S i S J 4 , 2 1 B s S i 4- ) “ 

WRITELNvFIy ’GIVE THE EXTENTION POINT OF THE EGUIPOTENT lAL SURFACE*); 
WRITELN( *GIVE THE EXTENTION POINT OF THE EQUIPOTENTI AL SURFACE*); 
READLN(Rly2i} ; 

WR ITELN (FI,RIsS = 4,2I:oS4:»; 

WRITELN(RlsSs4,ZliS:4) ? 


WR ITELN ( F I y ' GI v’E THE. C0NSTAf4TS OF THE FIRS>T BOUNDARY' LINE s a by-i“C=0 SUCH THAT THE* 
WR I TELi'4 (Fly * b i AR 7 1 NG P G I NT F ALLS 1 N THE R EG I ON s a :< by + c G 0 * ) ; 

WRITELN( ’GIVE THE CGr43TANTS OF THE FIRST BOUNDARY LINE s a;;+bv-J- SUCH THAT THE*)? 
WRIT£Lr4( *STARTir4G POINT FALLS IN THE REGION : ax-i-byY-c-.:>0 *) ? 

WR ITE ( * A“ * ) ; 

READLh4 ( AO ) ; 

WRITELN(A0;e:4) ? 

Wit ITE ( * B” * ) ; 

READLri ( BG ) ; 

WR ITELN ( BO s S 2 4 ) ; 

WRITE ( * C = * ) ; 

READLN ( CO ) ; 

WRITELN<C0s8;4) 

WRITE(FIy *A=* ) ; 

WRITELN(FIy A0;S:4) ? 

WRITE (Fly *B=* ) ? 

WRITELN(FIyB0sS:4) ? 

WRITECFIy *C=* ) ? 

WRITELN(FIyC0:8:4) ? 

WRITELNC ’GIVE THE CONSTANTS OF THE SECOND BOUNDARY LINE : ax+ b>H- c=0 SUCH THAT THE*); 
WRITELNC ’STARTING POINT FALLS IN THE REGION : ax + by+cOO ’) ; 

WRITE(’A=’)? 

READLN (A2) ? 

WRITELN(AEsS:4) ; 

WRITE( ’B=-* ) ; 

READLN (B2) ; 

WR ITELN ( BE ; 8 ; 4 ) ; 

WR ITE ( * C= ’ ) ; 

READLN (02) ; 


WRITELN(CE:Ss4) ? 

WR1TELN( ’GIVE THE RELATIVE PERMEABILITY OF THE INITIAL POINT*); 

READLN i EP ) ; 

WR ITELr4 ( " GI VE. TrlE! LENGTH OF INFINITSIMAL ELEMENT ALONG THE EGUIFDTENTI AL SURFACE*); 

READur4 (DL) ; 

WR ITELN ( DL - S « 4 ) ; 

WR ITELt'H ( F'l y ’ GIVE THE CGr4STANT3 OF THE SECOrJD BOUiY^DARV LINE « ax by •+-c = 0 SUCH THAT THE"); 
WR ITELN ( F'l y * STArdlNG P 0ir4T FALLS IN THE REGION s ax"i-by-i-c<—C *) j 


WR I TELN (FI, BH i S s 4 ) j 
WRITE (FI, i 
wRITELN(Fl ,CE:S;4 ) ; 

WR ITELt-4 ( FI , ’ GI vE THE RELATIVE PERMEABILITY OF THE INITIAL 
WRITELrM(FI,EPsS:4) ? 

WRITELr-4(Fl , "GIVE THE LENGTH OF INFINITBIMAL ELEMENT ALONG 
WRITELN(FI,DLS8:4;- ? 

500SIF ( (AO-J^RI / + (BO^-2I)-CO) <G.O THEN 
BEGI^4 lOF FIRST IF LGGF} 

IF ( ( A2*R I ) + ( ES^' 2 1 ;> +0 E ) < G . 0 THEN 
EEGir4 GOF SECOND IF i_OaF} 

ER 2 =0 . 0 ? 

E2s=0-C: 

FOR U:=1 TO NE DO 
BEGIN CGF FOR LOOP I 
ER : =ER^- ( CO ) #FER ( R Z , CHARGE , RI , Z I , U ; ? 

E2 5 =E 2+ C C G ) ‘^'FE Z ( R Z , C H AR G£, , R I , Z I , L’ ) ; 

END^ lOF FOR LOOP} 

IF EP=EFS THEN 

BEGIN lOF SECOND IF LOOP!' 

FOR Us = (NE-r1) TO EB DO 
BEGIN lOF SECOND FOR LOOP I 
ER ! ==ER + ( CO / -i^FER ( R Z , CHARGE , R I , Z I , U ) j 
E2 s=-E2+(CO)*rE2 (RZ , CHAR GE , R I , Z I , U ) j 
END? COP SECOND FOR LOOP} 

END 

ELSE 

BEGIN 

FOR U ! =■ ( EB'+- 1 ) TO EBB DO 
BEGIN 

ER ! =ER ( C G ) '>FER ( R 2 , C HAR GE , R I , Z 1 , U ) ? 

E2 : =E2+ ( CO ) *FEZ (RZ, CHARGE, RI ,21, U) ? 

END? 

END? {OF SECOND IF LOOP] 

WRITELr'4( ’ER=’ ,ERsa, ’ EZ=%EZ5S)? 

IF ER=0.0 THEN 
BEGIN 
DZ ? =0 . 0 ? 

DR:=DL; 

RI1:=RI+DR: 

RIE s =R I “-DR ? 

ZI:=ZI? 

RIEs=^E^PI”RIS: 

ZIE:=E*ZI--ZIS? 

D1 £=ABS(RI 1-RIE) ? 

D£ 5 “ABS ( RIE'—RIE ) 

Er4L 

E’.LSE 

BEGIN 

D2 ? = DL/SQRT ( l+SGR ( EZ /ER ) ) ? 

DR:=^~(EZ/ER)*DZ? 

RI1?=RI-J-DR; 

RI£:-^RI~DR? 

4 1 E - i. I “■ L Z ? 
r:. 1 E ? — ■ E F. 1 — R I L ; 

DE^'^E-ir \ F 4 4 ""F 4_L ■ ZIL'—Io-E; 


POINT" ) ; 

THE EQUIFOTENT I AL SURFAC 


2IS:=2I; 

IF DKDE THEN 

BEGIN 

Rl:=RI1? 

Zl:=2I1 

E^4D 

ELSE 

BEGIN 

RIs=RlE; 

2I;=ZI2 

END; 

IF <A0-i?-RI) + (E0^2I ;+CO>C.OOCOO THEN 
BEGIN 

XPt s == ( ( R 2 1~R IB ) -t-CO-s- ( R IS—R 1 ) ) / ( AO* ( R I— R IS ) +E;0* ( 2 1 — Z I S ) ) ; 

2IA- = C AO* ■: ZI3*RI~ZI*RIS ) +00* ( ZIS-ZI ) ) / ( BO* ( Z 1-2 IS ) +A0* ( Rl-RIS ) ) ; 
WRITELNC ’RIA=’ ,RIA:8:4, * Z1A==’ ,2IA:S:4) ; 

WRITELNCFl, 'R1A=’ ,RIA:Ss4, '' 2IA=’ .ZIAiSsA) ; 

END? 

IF ( AE*R I ) •4- ( BE* 2 I ) h-CE >0 . OOOOG THEN 
BEGIN 

RIA! = (BE*(R1B*Z1~R1*2IS)+CE*(RIS-RI/ • / ( AE* ( RI-RIS ) -i-BE* ( ZI -ZIB ) ) ? 
ZIAE = ( AE* ( ZIS*RI-ZI*RIS)-»'CE*i ZIS-ZI) ) / ( BE* ( 2 1~Z IS ) +AE* ( RI-RIS ) ) ? 
WRITELN( ’RIA==’ ,R1 h:Ss 4, * Z1A=‘ , ZIA:&s4) ? 

WRITELN<FI, ,RIA:Ss4: ' 2IA=:“ ,2IAs3:4) ? 

END? 

NR I TELN ( ’ R I =■ * , I s £ s 4 , ' 21=' . 2 1 2 £ : 4 ;• ? 

WRITELN; 

WRITELNCFI, ' RI= ' , RI : S : 4 , ' Z 1= ' , 2 I s S : 4 ) ? 

GOTO 500; 

END 
ELSE 
BEGIN 
GOTO 400? 

END?-t:OF OF SECOND IF LOOPI 

END 

ELSE 

BEGIN 

GOTO 400? 

END?{OF FIRST IF LOG?} 

400: WRITE-; 'DC YGU WANT Tu CALCULATE MORE EQUIPOTENTIAL SURFACE'): 
WRITELN ('IF YES THEN WRITE -1- ELSE WRITE-0~')? 

WRITE;FI,'DG you WANT TO CALCULATE MORE EQUIPOTENTIAL SURFACE')? 
WRITELN; FI, ‘IF YES THEN WRITE -1- ELSE WRITE-0~')? 

READLN(ANS) ? 

WRITELN ( FI T ANS ) ? 

IF’ ArMS= "I THEN 
GG i G LOO 
ELSE 

GGTLi 300 ? 

END 

ELSE 

BEGIN 


GOTO 300 
END ? 



WRITELNCFI ,Ar-4S.3) ; 

WRITEL^4;:AN33) ^ - 

IF ANS3=1 THEN 
BEGIN 

WRITELNKFI , ’GIVE THE NUMBER OF SUCH PDINTS’ ) ; 

WRITELNC 'GIVE THE NUMBER OF SUCH POINTS’)? 

RE ADEN ■( F J ? r4DM } ? 

WR ITElN (Fly f'4GM / ? 

WRIT£LN(NQM) ? 

WRITELNvFIy 'GIVE THE LOCATION AND SLOPE OF NORMAL AT SUCH FD 
WRITELN( 'GIVE THE LOCATION AND SLOPE OF ^4DRMAL AT SUCH POINT 
FOR riM2 = 1 TO NON DO 
BEGIN 

FOR 00: = 1 TXi 3 DO 
BEGIN 

READCFJ yNHCMMyOGj ) ? 

WRITE<riyHHCMMyG0G:S:4y ' ' ; ? 

END? 

WRITELNKFI; ? 

READLN ( FJ ) ; 

END? 


FOR MM:=1 TO NON DO 
BEGir'4 

CSOMM ,13: - ABS ( CDS ( AF'CTAN ( HHEMM , 33 ) ) > ? 

CSOMli , E3 : ~AB3 ( SIN ( ARCTAN ( HPICMM , 33 ) ) ) 

END? 

WRITELr4i 'TANGENTIAL FIELD VECTOR IS')? 

WRITcLi4 UFI » * TAi'4GEr4TI AL FIELF VECTOR IS*)? 

FOR LO:==-'i TO NGM DO 
BEGIN 

TAN E LOO : --C . 0 ? 

FOR P0:=1 TO NE DO 
BEGIN 

TANCLOO : =TANCL03+ COi-s-CHARGECLC , 1 j-^^FH ( RZ , HH , CS , PO , LO ) ? 
Er-4D ? 

FOR P0:=(NE+1) TO EB DO 
BEGIN 

TANCLOO :=TANi:LOO+CO-ifCHARGECLG, 1 0*FH ( RZ , HH , CS , PO , LO > 
END? 

WR I TELN ( F I , T ANC LOO: 12:6) ? 

WRITE CHHCLO, 1 0 : S : 3 , TANCLOO : 1£ : 6 ) ? 

WRITELN? 

Er-4D ? 

END 
ELSE 
BEGIN 
GOTO 60C 
END? . 


600:WRIT£LN(FIy 'END OF EXECUTION ')? 
C 1 os e ( f i ) ? 
c lose ( f j ) ? 


M CO 


XNPUT file = DAT. PRN 


tCCATloKiS mi:^DE 


“THE ELBCT^^ODE 


. 105 O . 09 - 
.,105 0.01 
.09 0 . 1 1 
3 0.1 1 
39 0.1 1 
105 0.09 
105 0.01-- 
.095 0.09- 
.095 0.01 
.09 0.09 
0 0.09 
09 0.09 
095 0.09 
095 0.01- 


L I_OOA1IONS* 


the air. 


cviAftiit uocA-rioNi IN Trte bietecTRic 

J 


0.1 1£1C 


^ p O-n eleotrcae- AiV Ui 2 ^-{aai fo^r^ts) 

-^Nut^be-r ifj- Po>r^±s o-r^ ^ v i 

Loca-tion Of iontou-r Points on ik-«- electro da. -A\V A>J5^-fp.ciL . 

.'5.. Humlper- 4 Po'f'ti ^ Air- Bi-efecir'c vrWe. •Uvn^er'tl'oJ (UnvvjDonen-tS 4 


.01 

0- 1 

'1 E 1 0 

.02 

0.1 

1E10 

.03 

C. 1 

1E1C 

,04 

0. 1 

12 10 

.05 

0. 1 

IE 10 

. 06 

0 . 1 

1E10 

.07 

0. 1 

IE 10 

,08 

0. 1 

1E10 

'O^i 

O' I 

iE to 


'( Elecix'*^ -fieldl io de^\re.(^ . 

I 

I 

y LocA-t^^v^ ^ ^ Pomts on tJ^ AvV-bCelecVic unl£^j^<js ^Tan^en-bicj 


ir^ , . ^e-iir^d 




HE ELECTR OD 


GIVE THE 

LOCATION OF SIMULATION CHAR-GES INSIDE 

0.0000 

0 . 3500 

0.0000 

0 . 3GGG 

0.0000 

0 . E5oG 

GIVE THE 

LOCATION OF SIMULATION CHARGES INSIDE 

“0 , 1 050 

0 . 0900 

-0. 1050 

0,01 GO 

-0,0900 

0 . 1 1 CO 

0 . 0000 

G . 1 1 00 

0.0900 

0. 1 100 

0.1050 

0 . 0900 

0, 1050 

0 . C 1 00 

GIVE THE 

LOCATION DF SlMUi-ATIuN CHAKLES Ir46lDE 

-0.0950 

0 . 0900 

-0.0950 

‘0 . 0 1 GO' 

— 0 . 0900 

0 . 0900 

0.0000 

0 . 0 9 0 0 

C.C900 

G . 0900 

0.0950 

0.0900 

C.095G 

G.0100 

GIVE THE 

NUMBER GF CONTOUR POINTS GN THE AIRv — 

GIVE THE 

LOCATION OF ABOVE CONTOUR POINTS 

C . 0000 

0 . 4O0C' 

0, 1000 

0.3000 

0 , 0000 

0.3000 

CONTOUR POINTS GN THE AIR— O' INTERFACE ARE 

-0. 1000 

C . 090'C‘ 

-0. 1000 

0 . G 1 GO 

-0.0900 

0 . 1 GOO 

0. OOO'O 

G . 1000 

0.0900 

0 . 100 0 

0. 1000 

0 . 0900 

0. 1000 

0 . 0 10G 


cr. 


THE AIR 


THE DIELECTR 


ir4TbRF"ACE 


GIVE THE RELATIVE PERMEABILITY OF D 

6,0000 

CHARGE VECTOR IS 
“1 -EE-15 
1. IE- 11 
1 .5E-14 
9-4E-16 
6 « 4E- 1 6 
-1-5E-14 
“E . OE- 14 
-1 :,5E-14 
9-4E-16 
6-4E— 16 
S-7E- 15 
S.7E- 16 
-Z-HE-14 
-E,GE-14 
-E - EE— 14 
£.7E.- 15 


U y iA’ drU i u 


y CA.. c u i 3. ij e 


Ui- V C 

GI vE' 


t;' r r y r , It' y 3 3 t r - y T's w i T 1 t j : 

F(_Lr4TS CiM 'iril AIR— c. 11-41 ERFACE 


th: 


0.0701 O.S'TG'I 
0.1 000 C . 300C 
0.0701 0.2399 

0.0000 0.2000 

GIVE THE NUMEiLFi GF TEST FGIuTS ON THE D-’ E INTLFGr ACE 
THE ERROR V'EVTuR IS 
~0 . OOOO 

0 . ooss 
— 0 . 0000 
0 . 0GS6 
-O , 0000 

CUHMULATIVE SQUARE ERf.GR IS 
• G.OOQE 

END OF ERROR C ALCULA'i j. ON 

DG YOU SaJANT TO CmLC'JLATE THE EQl* b.rJTI AL SURFACE- IF YES THE.N 
WRITE 1 ELSE WRITE C 

1 

calculatir^g t l-i e e q li i p c) t e n I i a 1 s u r f a c e 

give the initial point of The DESIRi D EGUIPGTENTIAL SURFACE 
O.OOGO 0.02 DC 

GIVE THE E. XT£!.r4 T IGTh FDlr*;! Lir 1 Hu. EQUIFDTENTI AL SUF^FACE 

0.0000 0.0200 

GIVE" Trii... O', ATnl S GF THE F1RE>T BGUNDARI' LINF s ax-i- b y + c-=G SUCH THAT THE 
•S T A R T I r IS P 0 1 N T FALLS I tY T H E R £( •> < 0 • i : a. -r b y + c % =• 0 
A=^' 1 . 0000 

Eh- O.COGO 

-0.1000 

GIVE, THL, CiJNGlA!*ITS OF THE BECOIsD EOUNDAFOv’ L 1 NE ! ax + by-i- c=C SUCH THAT THE 

STARTING PGINI' F’ALLS ITm THE REGIOI'i a/i + by h- C \ =-0 

A-"‘ O-OCOO 

B=- 1 .0000 

C- -”0. ICOG 

GIVE THE FtE-LAT 1 vE F EFiflEABILITY OF THE INITIAL. POINT 
6 . 0000 

GIVE THE LENGTH OF INFINITSIMAL ELEMENT ALONG THE EQUIPOTENTIAL SURFACE 

0.0020 


Rl = 

0 . 0020 

21 = 

0 . 0200 

RI = 

0 . 0040 

21 = 

0 . 0200 

RI = 

0.0060 

21 = 

0 . 0200 

RI = 

o.ooso 

21 = 

0 . 020 1 

RI = 

0.01 CO 

21 = 

G .CEO 1 

RI=^ 

0.0120 

21 = 

0 • 02CE 

RI = 

0-01 40 

21 = 

C . 0203 

RI = 

0-0160 

21 = 

G . 0204 

RI = 

0.0160 

21 = 

0.0205 

Rl = 

0 . 0200 

21 = 

0.0206 

RI = 

0 . 0220 

2 1 = 

C , G20c 

RI=- 

C m 0240 

21 = 

G . 0209 

Rl=. 

0.026G 

21 = 

0 . 02 1 1 

RI^ 

0 . 02 SG 

ZI = 

0 . 02 13 

RI = 

0 , 03 OG 

2 1 = 

0 . 02 1 5 

RI = 

0 - 03 1 9 

21 = 

0 . 02 1 7 

RI = 

0 . 0339 

21 = 

0 . 02 1 9 

RI=. 

0.0359 

21 = 

C « 022 1 

RI = 

0.0379 

2 1 =■ 

0 . 0224 

RI=- 

0 . 0399 

11 = 

0 • 0227 


0 . 04 1 9 

ZI = 

0 . 0229 

R 1“ 

0 . v.' 9 3 <“• 


C . 0232 

RI^ 

G , 09 56 

ZI = 

0 . 0235 


C . G*9 “ S 

ZI = 

0 - 0239 

R 1 = 

C . C‘4'9'L 

ZI = 

0 . 0242 


0 . ‘"'Ej 1 7 

Z i 

0 . 0246 



RI = 

0.0537 

ZI = 

0.0249 

RI = 

0.0557 

21- 

0 . 0253 

RI = 

0 . 0576 

ZI = 

0 . OE5>7 

RI = 

0.0596 

21- 

0 . 026 1 

RI=- 

0.06 15 

21 — 

0 . 0265 

RI=- 

0.0635 

2! I — 

0 « 0269 

RI = 

0.0654 

21- 

0.0274 

RI = 

0.0674 

21 — 

0 . 0279 

RI- 

0.0693 

21- 

0 . 02S3 

RI = 

0.0713 

21- 

0 . 02SS 

RI = 

0 . 073E: 

21- 

0.0293 

RI = 

0 . 075 1 

21- 

0.0296 

RI = 

0.077 1 

ZI- 

Q . 0*3 04 

RI = 

0.0790 

21 = 

0 . 0309 

RI = 

0 . 0809 

21- 

0.03 1 5 


0 . OSES 

Z X — 

0 » 0320 

RI“ 

0 . 0847 

2 1 — 

0 . 0326 

RI~ 

0 . 0S66 

21 — 

G . 0332 

RI = 

0 . CSSs 

Z I- 

C . C33S 

RI==^ 

0 . 0905 

2 1 — 

0 . 0344 

RI = 

0 . 09£C 

21- 

G . 035 1 

RI = 

0 . 09 43 

21 — 

0 . 0357 

RI = 

C.0961 

21- 

0 . 0364 

RI = 

0 . 09 SO 

Z i- 

0 . G37 1 

RI = 

0.G999 

ZI- 

0 . 037S 

RIA=^ 

0. 1000 liA 

- C.03' 

R 1 = 

0 « 10 16 

2 i =■ 

0 . 0385 

DC YOU WAiNT 

TO C 

ALCULATl 

GIVE 

THE IN 

ITIAL 

POINT i 

0. 

0999 0 

. 0376 


GIVE 

THE EX' 

TENT I 

ON FuIN' 

0. 

1000 0 

. 037 S 


give 

THE CONSTAN 

TS OF T! 


the. c.o!ui,r w I teP' 


NTS OF THE FIRST BOUNDARY LINE 
STARTING POINT FALLS IN THE REGION ! ax+by+cC^O 
A= 1.00GC 
B= O.OOOC 
C= — O.SOOG 

GIVE THE CONSTANTS OF THE SECOND BOUNDARY LINE J a;:H"by<" c— ( 
STARTING POINT FALLS IN THE REGION : a;; + by*i-c-:>0 
A= 0.0000 
B=- 1 . 0000 

C= -0.3000 

GIVE THE RELATIVE PERNEABILITV OF THE INITIAL POINT 

1.0000 

GIVE THE LENGTH OF INFINITSIMAL ELEMENT ALONG THE EQUIP£ 


:ax + by‘i 


FANTS OF THE SECOND BOUNDARY LINEJ< 
r FALLS IN THE REGION : aj; + by-5-c-:>0 


;:, + by-<-c-^ 


GIVE 

TH 

E LEI' 

4GTH 

OF INF I! 

0. 

002 

0 



RI- 

0. 

1019 

21- 

0 . G3S5 

Rl- 

0. 

1 037 

21 — 

0.0392 

RI- 

0. 

1056 

21- 

0 . 0400 

HI- 

0. 

1075 

21- 

0 . 0407 

RI = 

0. 

1093 

21- 

0.0415 

RI- 

0. 

1 1 12 

21= 

0 . 0422 

RI- 

0- 

1 1 30 

21 = 

0.0430 

RI= 

0 . 

1 146 

21- 

0 . 0436 

RI- 

0 . 

1167 

21 = 

0 . 0446 

RI- 

0 . 

1 1S5 

Z I — 

0 . 045>4 

RI- 

G. 

1203 

2 1 — 

C.0462 

RI- 

0. 

122 t 

ZI- 

0.047 1 

Rl- 

0 . 

123? 

21 — 

0 , 0479 

R 1 — 

0 « 

1 2 5'- 7 

i. j.. — 

Q . 04 So 

R I — 

0. 

127c' 

ZI- 

0 . 0496 



RI = 

0. 

1S93 

21 = 

0 . 0505 

RI = 

0. 

1 3 1 1 

21 = 

0 . 05 1 4 

RI = 

0. 

1329 

21 = 

0.05E3 

RI=. 

0. 

1 347 

21 = 

0 . 0532 

RI = 

0. 

1365 

21 = 

0 . 0542 

R 1=^ 

0- 

13SE 

21 = 

G.G551 

RI = 

0. 

1400 

21 = 

0 . 056 1 

RI= 

0. 

1417 

21 = 

0.0570 

RI= 

0. 

1435 

21 = 

0 - 0560 

RI- 

0. 

145E 

21 = 

0.059C 

RI = 

0, 

1469 

21 = 

0 . 0600 

RI = 

0. 

14S7 

21 = 

0.0610 

RI = 

0. 

1504 

21 = 

0 . 062 1 

RI= 

0. 

15H1 

21 = 

U . GoS 1 

RI= 

0. 

1538 

21 = 

G > G642 

RI = 

0. 

1555 

21 = 

0.0652 

RI= 

0. 

157E 

21 = 

0.0663 

RI= 

0. 

15SS 

21 = 

0 . 0674 

RI= 

0. 

1605 

21 = 

0 . 06S5 

RI= 

0 . 

1 622 

21 = 

0.0a 9 6 

RI = 

0 . 

163S 

21 = 

C . 0707 

RI=^ 

0 . 

1655 

21 = 

0 . 0 7 1 7 

R 1 = 

0 . 

1671 

21 = 

0 . 0730 

RI= 

c. 

1 6S7 

21 = 

0.07 42 

RI = 

0. 

1703 

21 = 

0 . 0754 

RI= 

c. 

1 720 

Z 1 

0.0766 

RI=:: 

0. 

1736 

21 = 

0 . C77S 

RI= 

0. 

175 1 

ZI = 

0 . 0790 

RI= 

0. 

1767 

j£. 1 = 

u . Go Go, 


0. 

1763 

Z 1 = 

0 . OS 1 4 

RI= 

0. 

1799 

Z 1 = 

0.0627 

R I~ 

0. 

IS 14 

2i= 

G . 0S39 

RI= 

0- 

1630 

21 = 

0.GS52 

RI= 

0. 

1845 

ZI= 

0.0865 

RI= 

0- 

1S60 

21 = 

0 . 087 8 

RI = 

0, 

1S75 

21 = 

0 . OS 9 1 

RI= 

0. 

1891 

21 = 

0.0904 

RI= 

0. 

1 906 

21 = 

0.G917 

R I~ 

0. 

1 920 

2 1 = 

0 . 093 1 

RI = 

0. 

1935 

21 = 

0.0944 

RI= 

0. 

1950 

21 = 

0 . 0958 

RI=^ 

c. 

1 965 

21 = 

0 . 097’ 1 

RI== 

0. 

1 979 

21 = 

C . 0965 

Rl = 

0. 

1993 

21 = 

0.0999 

R1A= 

c 

.2000 ZIP 

1= 0 , 10- 

RI= 

0. 

2008 

Z 1 = 

C . 1 G 1 3 


DO YOU WANT TC< CALCULATE MORE ECUIPOTENTIAL 3URFACE1F YES THEN WRITE -I— ELSE WR1TE~G— 
1 

GIVE THE INITIAL POINT OF THE DESIRED EGUiPOTENTI AL SURFACE 
0 . 0 G 0 0 G . G 4 C C 

GIVE THE EXTENTION POINT OF THE EGUIPOTENTI AL SURFACE 
O.OOGC 0.0400 


GIVE THE CONSTANTS OF THE FIRST BOUNDARY LINE ;ax+by+c=G SUCH THAT THE 
STARTING POINT FALLS IN THE REGION : ai-;+by-^c<=G 
A= 1 . 0000 



GIVE THE RELATIVE FERMEABILITV OF THE INITIAL POINT 
^000 

INITSIMAL ELENEraT ALONG THE EGUIPGTENT I AL SURFACE 


6. 

0000 



GIVE 

THE LEr4GTH 

OF INF If 

0. 

oceo 



R 1 = 

0 - OOEO 

21 = 

0 . 04'00 

RI=- 

C « 004C 

2 1 = 

o 

o 

o 

o 

RI=-- 

C . 0060 

ZI = 

0 . 0400 

RI=^ 

0 « oosc 

21 = 

C . 040 1 

RI= 

0.0100 

2 i = 

0 . G4C 1 

RI=^ 

0.0120 

21 = 

0.040E 

RI = 

0.0140 

21 = 

0 . 0403 

RI = 

0.01 60 

Z 1 = 

0 - 5^404 

RI = 

0-01 SO 

21 = 

0.04G5 

RI = 

C - 0200 

ZI = 

0.0407 

RI = 

0 . 0220 

Zl = 

0 . 0406 

RI=^ 

0.G24G 

Z i = 

0.0410 

RI= 

0 . 0260 

21 = 

0 . 04 1 2 

RI=- 

0 * C2oG 

21 = 

0.04 1 4 

RI = 

G . 0299 

Z i = 

0 . 04 1 6 

RI = 

0.03 99 

21 = 

0 • C 1 6 

RI = 

0.0339 

21 = 

0 . 042C 

RI=^ 

0.0359 

21 = 

C . 0423 

RI = 

G . 0379 

21 = 

0 . 0426 

RI=- 

0.0399 

2 1 = 

0 » 0426 

RI = 

0 . 04 1 S 

2 1 =■ 

0 . G43 1 

Rl = 

0 - C43S 

21 = 

0 .0435 

Rl = 

0.045S 


0 . 0436 

RI=^ 

0 . 0476 

21 = 

0 . 044 1 

RI= 

0 - 0497 

21 = 

0 . 0445.- 

RI=^ 

0.0517 

21 = 

0 . 0449 

RI- 

0 . 0536 

i.- JL ^ 

0 . 0453 

RI = 

0.055o 

L J. = 

■0 .0457 

RI= 

0 . 0576 

21= 

0 . 0461 

Rl= 

0.0595 

2 1= 

C . 0465 

RI=^ 

C . C6 1 5 

21 = 

G.047G 

RI= 

0.0634 

ZI= 

0.0475 

RI = 

C.0653 

21 = 

C - 0460 

RI=- 

0.0673 

21= 

0.0485 

Rl= 

0-0692 

21= 

0.0490 

RI= 

0.071 1 

21= 

0 . 0495 

RI = 

0 . 07 3 1 

21 = 

0 . 050 1 

RI=^ 

0.0750 

21 = 

0 . 0506 

RI= 

0.0769 

21 = 

G.C512 

RI = 

0 « 07 ScF 

21 = 

0.0516 

RI= 

0 . 0SC7 

21 = 

C • u524 

RI= 

0. 0626 

21 = 

0.0530 

RI= 

0.0S45 

21 = 

0.0536 

Rl^ 

0 « 0664 

21 = 

0 « 05>43 

Rl = 

0 , 0663 

21 = 

0.055G 

Rl=^= 

0 . C902 

21 = 

0.0556 

Rl= 

0-0920 

21 = 

0.0563 

Rl=^ 

0 . 0929 

21 = 

0 . 057 C 

Rl = 

0 . 0956 

21 = 

0.0577 

RI= 

0 . 0976 

21 = 

0 . 0585 

RI = 

0.C995 

21 = 

0 . 0592 

RIA=-^ 

0 . 1C0C 

F Z 1 A 

= 0- . V 5 ' 

R I “= 

U . '1 0 '1 4 

21 = 

C . C6CD 

LL‘ YuU LAF'-L 

TQ C 

Hi_C UuAI i 

GIVE 

THE IF4I 

T I AL 

POINT i 

0. 

0995 0- 

0 1; 9 E 


GI'v'E 

T'HL EX'* 

"ENT I 

1,7' 4 rEIi'l 


; c:.)’'-. i iMi_ 


SURFACE 


/ 


0.1000, 0.0594 

GIVE THE CONSTHrMTS OF THE FIRST BOUNDARY LINE : + c = G SUCH THAT THE 

STARTING POINT FALLS IN THE REGION « a M + by*i-c<=0 
A= 1.0000 
B=- 0.0000 

C= -O.EOOO 

GIVE THE CONSTANTS OF THE SECOND BOUNDARY LINE s + t'y + c = 0 SUCH THAT THE 
STARTING POINT FALLS IF-i THE REGION s ax + bya- c< =0 
A= 0 . 0000 
1.0000 
C= “0.3000 

GIVE THE relative PERMEhEILITY OF THE INITIAL POINT 
1 . OOOG 

GIVE THE LENGTH OF INFINITSIHAL ElENENT ALONG THE EQUIPGTENTIAL SURFACE 
O.OOHG 


RI = 

0.10 IF 

2 1 ~ 

M 'w*C'0 \ 

RI = 

0. 1037 

21=- 

0 . 0609 

Rl = 

0 . 1056 

21=- 

0 . 06 1 7 

RI=.-: 

0. 1074 

21 — 

0 . 0 6 2' 4 

Rl=^ 

0 . 1093 

ZI=- 

0 . Oci E E 

RI=r- 

0 . 1 -i 1 1 

21 — 

0.0640 

KI== 

G , 1 139 

2 1 =■ 

0 . Gc'4S 

RI~ 

0-1 148 

2 j..= 

0 . 0o56 

Ri=-^ 

0 « 116 6 

ZI = 

0 . 0664 

Rl=^ 

0. 1 164 

2 i =• 

0 . 06 7£ 

RI = 

0 - 1 EOS' 

Z i =! 

0 . OsS 1 

Rl = 

0 ■ 1 EEO 

Z I- 

0.,0oS9 

RI=^ 

0 . 1 E 3 £> 

21 = 

0 . Do 9S 


0. 1E56 

11 = 

0 , 07 07 

RI=- 

0 . 1 ‘E 7 8 

21 = 

0 , C '7 1 L 


C„, 129E 

21 = 

0 - G7'£5 

K 1 = 

0 . 1 3 0 7 

21 = 

0 . 0 /'■ 6 4- 

RI = 

0 . 1GE7 

21 = 

0 .0744 

RI = 

0. 1345 

21 = 

0 . 0754 

RI=^ 

0 « '1 3 6 E 

21 = 

0 . 07 63 

RI = 

0 • 1379 

21 = 

C . 0 /' 7 3 

RI>- 

G - 1397 

Zl = 

0 . C7SS 

RI = 

0 • ’1 414 

21 = 

0 . 07v4 

RX=^ 

0 . 143 1 

21 = 

0 . 0604 

RI = 

0 . 1 4 4’ S 

21 = 

0 . OS 1 4 

RI = 

0. 1465 

21 = 

0 » 0o35 

RI“ 

0 - ; 4oE 

21 = 

0 . 0636 

RI = 

0 » '1 49S 

Z 1 = 

G . 0 S 4 7 

RI^ 

0. 15 15 

21 = 

0 . osss 

RI=- 

0- l5sE 

ZI = 

C . 0369 

RI = 

0. 154S 

2 1 = 

0 . 0S3 1 

RI=^ 

0. 1564 

21 = 

G . Ob 96 


U . 1 5 6 1 

21 = 

0 . 0904 

RI=> 

0 . 1 5 9 7 

21 = 

0 K ’w‘9 1 5 

RI---= 

0 . 1 1 3 

2 1. =•■ 

0 . 0937 

R I ™ 

0. 1639 

21 = 

0 . 0939 

R I ~ 

0 - 1645 

21 = 

0 . 0953 

RI = 

0 . '1 66 1 

Z 1 = 

C . 0V64 

RI = 

0. 1676 

21 = 

0 , 097 6 

RI=' 

0 . 1693 

21 = 

0.0939 

RI-"~ 

U . '1 / U7 

21 = 

0 . C ■- 

R I - 

. 1 . L_ 

21 = 

u . 1 w.' 'i 4’ 

Rl-=^- 

0 . 1 7 3 9 

2 1 = 

G - 1 037' 

RI=^- 

C - ■■ 7 1 3 

2i = 

C .. ■ 0-C 

Rl=^ 

U . ", /’'c-C' 

21 = 

C „ 101 3 

Rl=- 

c . *:-s3 

21 = 

0 . 1 3 c ” 

RI = 

0. 1796 

Z I = 

C'. . 1 C c 0 


RI= 

0- 1813 

21=- 

0 . 

1093 

RI~ 

0, 1827 

Zl^ 

0. 

1 1 07 

RI = 

0. 1S4E 

21 = 

c . 

1121 

RI~ 

C. 1856 

21 = 

G , 

1 1 35 

RI = 

0- 1871 

Z I = 

G . 

1 1 49 

Ri = 

0.1665 

Z I = 

G . 

1 1 63 

Rl= 

0. 1699 

21 = 

G . 

1 177 

RI = 

G . 1 9 1 3 

21 = 

0. 

1 1 9 1 

RI = 

0 . 1 987 

21 = 

0 . 

1 206 

RI~ 

0.194 1 

21 = 

0 . 

1 EEC 

RI = 

C. 1954 

21 = 

G . 

■| -p c. 

RI=- 

0.1 968 

21 = 

0 . 

1250 

Rl= 

0.1 96 1 

21 = 

C' . 

1264 

R 1=- 

0. 1995 

2 1. =' 

0. 

1 2 7’ 9 

RIA= 

0 . EOG' 

U i X 4 

1 = 

0, 12 

RI = 

0.20GS 

Z I = 

G • 

1 2 V 4 



DO YOU wANT 
•1 


GIVE THE 
0 . 0000 
GIVE THE 
# 0.0000 
GIVE 'THE 


TO CALCULATE MORE EQUIFOTENTIAL BURFACEIF YES THEM MR 
_ . IMT OF THE DESIRED EGUIPGTEr-4TIAL SURFACE 


ITE “1“ 


0 « OqUI-' 
EXTEI'mT iC 
C . 0600 
CONST AN! 


INI ' 


OF 




rSMTIAL 


OF 


THE FIRST BOUNDARY LINE 
IN THE REGION J a;;+by+c--. =C 


SURFACE 
; a;-i+b v-s- C 


BUCK THAT THE 


STARTING POINT FALLS 
A= 1 . GOOD 
B- 0,0000 
C = — G , 1000 

GIVE THE CONSTANTS OF THE SECOND BOUNDARY LINEsax+by- 
STARTIN 6 POINT FALLS IN THE REulQN ^ axTby “t"CL=0 
A = 0-0000 
B=- 1 , GOOO 

C = - 0 . 10 G 0 

GIVE THE RELATIVE PERMEABILITY OF THE INITIAL POINT 


SUCr 


THE 


6. 

0006 



GIVE 

THE LET 

4GTH 

OF INFir 

0. 

GOEC 



RI = 

O.OGSO 

Z 1 -■ 

0.C6CG 

RI=:^. 

0.0040 

Z 1 “ 

0 - 0600 

RI = 

^ 0.0060 

ZI=^ 

C . 0600 

RI=. 

0.0080 

ZI = 

0 . 060 1 

RI = 

0-C10C 

21 = 

0 . C60 1 

Ri= 

0.0120 

Z 1 

0 » C‘602 

RI=- 

0.0140 

21 = 

C . 0603 

RI=^ 

G « C i c»U 

21 = 

G « 0604 


0.01 SC 

ZI = 

U . UoCf 5 

RI=- 

0 . OEIOO 

21 = 

0 . 0607 

RI = 

0 . 0220 

21 = 

0 . 0608 

RI= 

0.0240 

21 = 

0.0610 

RI = 

0 - 0260 

Z I = 

0 . 06 1 2 


0 ■ 0280 

Z 1 = 

0.0614 

RI = 

0.0299 

21 = 

0.0616 

RI = 

0.0319 

21 = 

0 . 06 1 S 

RI= 

0.0339 

21 = 

0 . 062 1 

RI=- 

0.0359 

Z I = 

0 . 0624 

R I- 

0.0379 

21 = 

G n 0626 

RI==' 

0.0396 

21 = 

0 . 0629 

RI = 

0.0418 

21 = 

0 . C‘6ZB3 

RI=^ 

0.0438 

21 = 

G . 0a3c 

RI = 

0.0458 

21 = 

0 . 0640 

RI=- 

G . 047^7 

21 = 

C . 0643 

RI= 

0.0497 

21 = 

0.-G647 

RI=:. 

0.0516 

21 = 

0 . C)65 '1 

RI = 

0 . 0536 

21 = 

0 . 0656 

R I~ 

0.0555 

Z i. “■ 

0 . 0o60 

RI = 

0-0575 

ZI = 

C . 0665 

RI=: 

0 . 0594 

21 = 

C . Ohio 

RI = 

0-0614 

ZX = 

0.0675 

RI = 

0 . 0633 

21 = 

0 . 0680 

RI- 

0.0652 

21 = 

0.0665 

RI=- 

0 . 067 1 

21 = 

0 .069 'i 

RX^ 

0 « 069 1 

21 = 

0 , 0697 

R2=- 

0,07 10 

Z I 

'■J . U / U C' 

R 1 = 

G , 0729 

Z 1 = 

0.0709 

RI=- 

G . 07 4£ 

2 i "■ 

0 V 0 / 1 5 

Rl^.^ 

0 - 07 67 

2 Z. = 

0 . 0726.' 

,R I =■ 

G . 076 f 

Z i = 

0 , 0728 

RX~^ 

0 « 0£C4 

2 1 

C . 07' '2 5 


G » UcLct 

2 1 =' 

/ 81 ;:. 


ALONG THE EQUiPUl 


r-i iwC- 


ELSE W 


RI = 

0.0S4E 

21 = 

0-0749 

RI = 

O.OS60 

21=- 

G.G756 

RI = 

0-0679 

Z I = 

0 . 0764 

RI== 

0.0896 

Z I = 

0 . 077 1 

RI“-= 

C . 09 1 6 j 

21 = 

0.0779 

RI = 

0.093D 

21 = 

G . 07 66 

RI = 

0.0953 

21 = 

0 . 0794 

RI = 

0.0971 

ZI = 

0 . OSOE 

RI = 

0-0990 

21 = 

0 . 0609 

RIA== 

0. 100C 

' ZIA^ 

= G . 08 1 4 

RI = 

0.1006 

21 = 

0 - 06 1 7 

DO YOU WANT 

TO CALCULATE MORE 

1 

GIVE 

THE INI 

, ! lAL 

POINT Or THE 


i'ES THEN 
RFACE 

0.0990 0.0S09 

GIVE THE EXTENTION POINT OF THE EQUIPGTENTI AL SURFACE 

0.1000 G.0S-.4 

GIVE THE CONSTANTS. OF THE FIRST BOUNDARY LINE i a>;-+-by+ c=0 SUCH 
STARTING POINT FALLS IN THE REGION s + c.<~0 

A= 1.0000 

B=- 0.0000 

C= “O.EOOO 

GIVE. THE CONSTANTS OF THE SECOND BOUNDAR Y LINE '• aK + by+ c=Q SI 
STARTING F'0It'4T FALLS IN THE REGION s ax + byT-c<’.=-0 
A= 0.0000 
B= 1.0000 

Qsz —0-3000 


-CRITE -1- ELSE 


THAT 


THAT 


THE 


GIVE 

THE RELATIVE PERriE. 

1 . 

0000 



GIVE 

THE LEr 

vGTH 

Gr lilrl 

0. 

0020 



RI = 

0. 10EC 

21 = 

0 . Oo '1 7 

RI = 

0.1039 

21 = 

0 . 0Sc£ 

RX = 

0. 1056 

21 = 

0 . 0627 

Rl = 

C. 1077 

21 = 

0 . 0634 

RI = 

0. 1096 

21 = 

G . 064 1 

RI = 

0.11 15 

21 = 

0 . C34S 

RI = 

0. 1133 

21 = 

C . 0656 

Rl = 

0. 1151 

21 = 

0 . 0664 

Rl = 

0. 1 170 

21 = 

C « 0£7£ 

Rl = 

0.1166 

Z 1 = 

C.CS61 

RI= 

0. 12G6 

21 = 

C . Ooo9 

RI = 

G . 1HE4 

ZI = 

0 . OSv S 

Rl = 

0. 1241 

21 = 

'U m U T ’w / 

RI = 

0.1 £5 7 

21 = 

G • 09 1 7 

RI = 

0. 1277 

21 = 

0 - 09£6 

RX = 

0. 1294 

21 = 

0 . 0936 

RI= 

0. 131 1 

21 = 

C . 0946 

RI = 

0 . '1 3 cl 9 

21 = 

C . 0956 

Rl = 

0, 1346 

21 = 

C.0967 

RI = 

C . 1 363 

21 = 

0 - 0977 

RI = 

0. 1360 

21 = 

0.09SS 

R 1 = 

0 . 1397 

21 = 

G . 0999 

Rl = 

0. 14 13 

21 = 

0 . '1 0 1 c 

RI = 

0. 1430 

ZI = 

0 . 1 c£ 

RI = 

0. 1446 

21 = 

0. lOEZ 

R I = 

G « 1465 

2 1 

C’ 104 4 

RI = 

« ' 1 

2. 2 

0 . 1055 

R I = 

0 . 1‘ T 5 

2 2-' 

. : O'U .■ 


0 . '1 5 1 1 


C « '1 2 7 9 

Rl = 

0 - 15E7 

2 2 = 

0 '1 09 

Fv 1 

C. 1545 

2 2'= 

0 . '1 1 03 


;:)lNAiL EL 


NENT 


HL_Ui.:U ; DC. C.. UiUi.ru I C-lx ! xWi. 


= ACE 


wR I TE""0” 


RI=; 

0. 1559 

ZI = 

0.1115 

RI= 

0. 1574 


0.11 28 

RI= 

0. 1590 

21 = 

0 . 1 1 4 1 

RI= 

0.1605 

z x= 

C . 1 153 

RI = 

G . 162 1 

Z I- 

0 » 1 1 66 

RI=- 

0 . 1636 

ZI = 

0. 1 179 

RI= 

0.1651 

21 = 

0.1 192 

RI=: 

0. 1666 

21 = 

0, 1206 

RI = 

0. 1681 

21 = 

0 . 1 £ 1 9 

RI~ 

0.1695 

Z 1 = 

0,1 233 

RI = 

C. 1710 

Z I = 

G . 1 246 

Rl= 

0 . 1724 

21 = 

0 « 1 26 G 

RI = 

0. 1739 

Z 1 = 

0 . 1 274 

RI = 

0. 1753 

21 = 

0.1288 

RI = 

0. 1767 

21 = 

0 . 1302 

RI = 

0. 1781 

Z 1 = 

0.1317 

RI = 

0. 1795 

Z 1 = 

0 . 133 1 

RI=^- 

0. 1809 

ZI = 

0.1 346 

RI= 

0. 1S22 

Z I = 

0.1360 

Rl=:. 

0 . 1836 

Z 1 = 

C. 1375 

RI = 

0 - 1849 

Z 1 = 

G . 1 3 9 0 

RI== 

0 . 1863 

Z 2 = 

0.1405 

RIs= 

0 . 1876 

21 = 

0 . '1 H 0 

RI= 

0. 1889 

21 = 

G . 1 435 

RI- 

0. 1902 

21 = 

0. 1450 

RI= 

0. 1914 

21 = 

C . 1 466 

RI= 

0 . 1927 

Z I = 

0 . 1 48 1 

RI= 

0- 1939 

21 = 

0.1497 

Rl = 

0. 1952 

21 = 

0.1513 

RI=^ 

0 . 1964 

Z 1 = 

G . 1 528 

RIs 

0 . 1 97 6 

ZI = 

G. 1544 

RI=. 

0.1988 

ZI = 

0. 1560 

RI = 

0 . 2000 

ZI = 

0.1577 


RIA= 0.2000 2IA= C.1577 
RI~ 0.2012 ZI= 0.1593 

DO YOU WANT TO CALCULATE MORE EQUIFOTENTIAL SURFACEIF YES THE^4 WRITE -1- ELSE WRITE-0' 
1 

GIVE THE INITIAL POINT OF THE DESIRED EQUIPOTENTIAL SURFACE 
0-0000 0.030C 

GIVE THE EXTENTION POINT OF THE EQUIPOTENTIAL SURFACE 
0.0000 O.OSOC 

GIVE THE CONSTANTS OF THE FIRST BOUNDARY LINE :d.x+by+c=C SUCH THAT THE 

STARTING POINT FALLS IN THE REGION : + c<=0 

A== 1.0000 

B==- 0-0000 

C= “0. 1000 

GIVE THE CONSTANTS OF THE SECOND BOUNDARY LINE : ax+by+c=0 SUCH THAT THE 
STARTING POINT FALLS IN THE REGIOiv s a by + c<—G 
A=- 0 ■ OOOC 

B= 1 . 0000 
C= -0-1000 

GIVE THE RELATIVE PERMEABILITY OF THE INITIAL POINT 


6.0000 

GIVE THE LENGTH OF INF INITSIMAl. 

0 « 0020 

RI- G-004C Z, 1“ G.06D0 

R I " Lf' « W i-f' O fi-f X ~ . 1, " O L.'* 

R X **'■■ O - 0 W C -L '«J' W* 

R X ™- 0 . C 00 2 1 0 . OSO ^ 

RX=-^. G.u‘.EC C;.GO.O:'1 


ELEMEr4T ALOrYG THE EQUIPOTENTIAL SURFACE 


RI = 

0.0160 

ZI = 

0 . 0603 

RI = 

0.0180 

2I~ 

0 . 0803 

RI= 

0 . OECC 

21 = 

0 . 0804 

Ri= 

0 « 0E20 

21 = 

0 . 0806 

RI = 

0.0240 

21 = 

0.0807 

RI=^ 

C.026D 

21 = 

0.0609 

RI= 

0.02S0 

21 = 

0 . 06 1 '1 

R 1 = 

0.C30C 

ZI = 

O.OS13 

RI = 

0.0319 

21 = 

0.0815 

RI=^ 

0.C339 

21 = 

0 . 08 1 7 

RI= 

0.0359 

21 = 

0 . 0820 

RI= 

0.0379 

21 = 

0 . 0623 

RI = 

0.0399 

2 1 = 

C . 0S26 

RI= 

0.0418 

21 = 

0 . CS3C 

RI = 

0.043S 

21 = 

0 . 0833 

RX=- 

0.0458 

21 = 

C . 0837 

RI= 

0.0477 

21 = 

C . 084 1 

RI=^ 

0.0497 

21 = 

0 . C845 

RI = 

0,0516 

21 = 

0 . 0850 

RI= 

0.0536 

21 = 

G .0855 

RI= 

0.0555 

21 = 

G . OSoO 

RI==- 

0.0574 

21 = 

0 . 0865 

RI^ 

C . 0594 

2 !•= 

0.C67G 

:R1 = 

G • 06 ■'! 3 

21 = 

0 ,087 6 

RI=^ 

0.0632 

Z 1 = 

0 , CooE 

RI = 

0 . 065 1 

21 = 

G . GS68 


0 . 0670 

il X~ 

0 . 0895 

RI=- 

0 . 0688 

2 1 =■ 

0 . 0902 


0 . 0707 

21“ 

0,0909 

RI = 

0 . 0726 

21 = 

0.09 ‘7 

RI== 

0 . 07 44 

21 = 

0 . 0925 

RI~ 

0 . 0762 

21 = 

0 . 0933 

RI = 

0,0780 

21 = 

0 . 0942 

RI = 

0.0798 

Z 1 = 

0 , 0952 

RI = 

0.0S15 

21 = 

0 . 0962 

RI= 

0 . 083 1 

21 = 

0 « 0973 

RI = 

0.0847 

Z 1 = 

0 . 0985 

RI = 

0 . 0862 

21 = 

0.0999 

RIA^= 

0,0663 ZIA 

= C . 1 0C’i 

RI=^ 

C.0S7S 

ZI = 

0 , 1C 15 

Du YOU WANT 

TG C 

ALCULATE 

GIVE 

THE in: 

ITIAL 

POINT Of 

G . 

C'S62 C 

. 'w'9‘t 7 ' 


GIVE 

THE EX TENT 1 

GN F'GINT 

, C. 

0863 0 

. 1 GOO 


GIVE 

THE CGi 

1 f-j , H 

T3 GF THi 

STAR 

TING Pu. 

IfvT FALLS Its 

A=- 

1 , OGOu 



B=^ 

0 « 0000 



c= - 

0 . 2000 



GIVE 

i HE CGI 

’4 ST AN 

TS OF THi 


3URFACEIF YES THEN WRITE -I- 


WR] 


EGUIr CTENTI^ 




OF THE EQUIFuTENTIAL 


r j.n.w i r'iJ*w‘UjL-'i-in i 
E R E G 1 G r'4 * a y-. b v ■' 


c <=0 


3Yi£||^7lr^j(3 pr 

0*0C‘0C 
£- l.ODC 

~ G . 3 c ; 

uX'v'E. "’.'HE. 

“ » GuGu 

GI\'E "THE w_Er. 

0 . 


FALLS 


E SECOND BOUNDAFiY LII4E 
THE R EG I ON s a b y + c <! — C 


t ur 


SURF ACE 
: as-i + bv+c^O 


a.;-; + bv"^ c =0 BL 


THE 


THAI 


C.%jrW j. r L.‘ i C-IV . 


Rl= 

0.09E1 

Z 1 = 

0 . 1 016 

RI = 

0.0941 

ZI=^ 

0.1019 

RI = 

0.0960 

Z 1 = 

0. 1022 

RI=- 

0 . 09BG 

Zl = 

0 . 1 026 

RI = 

0 . 0999 

Z 1 = 

0 . 1 0 3 1 


0- 1018 

ZI~ 

0 . 1 039 

RI = 

C . 1 037 

2I~ 

0. 1043 

RI = 

0. 1056 

Zl=- 

o 

q 

O 

RI=- 

0. 1075 

21“ 

0 - 1058 

RI = 

0.1093 

21=. 

C . 1 066 

RI = 

0.11 1 1 

Z 1 — 

C . 1074 

RI=^ 

0.1 129 

2 1 " 

0 « 1 0 S 3 

RJ.= 

0. 1 147 

21- 

G . 1 0'-/c‘ 

RI=- 

0.1165 

Z I — 

0.1 101 

RI = 

0.1 182 

21 = 
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RI= 

0.1817 

ZI = 

G - 2737 

RI= 

0.1820 

21 = 

0 . 2757 

RI== 

0 . '1 823 

21 = 

0 . 2777 

RI = 

0-1825 

21 = 

0 . 2797 

RI=^ 

0.1827 

21 = 

C . 2 8 17 

RI = 

0.1829 

21 = 

0 . 2836 

RI=- 

0, 1831 

2 1 = 

G . 2856 

RI = 

0. 1833 

21 = 

0 . 267q 

RI = 

0 . 1834 

21 = 

C . 2896 

RX = 

G - 1 835 

21 = 

0.2916 

RI = 

0.1836 

21 = 

0 . 2936 

RI- 

0. 1837 

ZI = 

0 . £956 

RI” 

0. 1838 

Z I = 

0.2976 

RX= 

0- 1836 

Z I =•■ 

0.2996 

RIA= 

0. 1S3i 

B 2IA 

= O.SOi 

RI = 

0. 1836 

71 = 

0 . 30 16 

DO Y 

DU WANT 

TO CALCULAT: 

WRITE 1 EL: 

1 

c a 1 c a I a t i r'! c 

SE WF 

HE C 

the 

equi pat i 

GIVE 

0. 

THE IN 
0000 0 

itial 

. 14CG 

FGINT i 

GIVE 

0. 

THE EX 
0000 0 

TENT I 
« 1400 

ON FOIN' 


M . t i a I 

jF the 


IF YES THEN 


s u r i a L t 

3ESIRED EGUIFGTENTIAL SURFACE 


GIVE THE C0r^4STANTS OF THE FIRST BQUr4DARY LINE :aK+bv + c=C SUCH THAT THE 
STARTING POINT FALLS IN THE REGION s aK+by+c<=G 
A= 1 . 0000 
0.0000 
C= -0-2000 

GIVE THE CONSTANTS OF THE SECOND BOUNDARY LINE:a;'i + by + c=0 SUCH THAT THE 
STARTING POINT FALLS IN THE REGION s a.K + by + C'--.=0 
A= G.OOOO 
B= 1.0000 
C= ”0-3000 

GIVE THE RELATIVE FERriEABI;-!! V OF THE INITIAL F’UlNT 
1-0000 

GIVE THE LENGTH OF INFINITSINAL ELEMEN ■ ALGlYGi THE EQUiPLH ENTIAL £!URFACE 
0.0030 
0.0030 
0 . GC 6C 


RI = 

RI =^ 

RX= G.0G90 


-L JL - 

ZI = 


RX= 

RI~ 

, KI ==' 


0 « 0 1E0 ZI = 
0.0 i DO ri= 

0 ■ C '“i o 0 Z I - 


0 . 1 40 C 

0 . 1 4GE 
Go 1403 


Rl=- 

FI - 




RI=r 

RI^- 
RI = 


0.053 

0.05c.( 


t-' . "t H-OL- 

G . 1 4 7 3 
C. 1433 


RI = 

C.G5S8 

Z 1 — 

C. 1503 

Fi 1 

0.06 ' c> 

Z I 

D . 1513 

RI=- 

C . 0644 

Z I “ 

0 . 1525 

RI===^ 

0 . 067 

Z — 

0 . 1 3'co 

RI=-~ 

0 . 069 S 

Z 1=^ 

0. 1549 

RI = 

0 . 0736 

21 = 

0.1562 

RI = 

0.0752 

ZI = 

0.1575 

RI = 

0.0779 

ZI = 

0. 1589 

RI=^ 

O,0SC5 

Zl = 

0. 1603 

RI = 

0 . Cq3c 

ZI = 

0.1618 

RI==^ 

0.0357 

21 = 

0 . 1 633 

RI = 

0.0863 

Z 1 = 

0 . 1649 

RI=^ 

0 « 0908 

Z 1 = 

0 . 1 665 

RI = 

0.0933 

21 = 

0 . 1682 

RI=' 

0.095c 

2 1 = 

0.1 699 

RI = 

C . 0962 

ZI = 

0.1717 

RI=^ 

0. 1006 

2 1 = 

0 . 1735 

RX = 

0 « 1CE9 

i! 

0. 1753 

RI=^ 

0- 1053 

Z 1 

0 . 1772 

RI = 

0. 1076 

ZI = 

0. 1792 




RI = 

0. 1096 

21 = 

C. 1£1 1 

RI = 

0.1120 

ZI = 

0 - 1 S3E 

RIn 

0. 1142 

21 = 

0. 165E 

RI = 

C . 1 1 63 

21 = 

0. 1S73 

RI=^ 

0. 1 1S4 

21 = 

0. 1695 

RI = 

0- 1205 

2 I = 

0 . ‘1 9 1 7 

RI-= 

0.1 220 

ZI = 

G. 1939 

RI = 

0- 12 A- 3 

2 X = 

0.1 96 1 

RI^ 

C • 1 L64 

Z 1 = 

G. 1964 

RI = 

0-1 2S3 

ZI = 

C . cc’vyc 

Rl = 

0- 1301 

21 = 

U . E'03 1 

Rl = 

0. 1319 

ZI = 

0 . 2055 

RX=^- 

0 . 1 337 

ZI = 

0 . 2080 

RI=- 

0. 1354 

ZI = 

0.2104 

RI = 

C . < vl? < • 

Zl = 

0.2 129 

RI=^ 

0 . 1 3ci-7 

21 = 

0.2155 

RI = 

0. 1402 

21 = 

0,2160 

RI = 

0 . 1418 

Z 1 = 

0 . 2206 

Rl== 

0. 1432 

ZI = 

0 . 2232 

RJ=r 

0 . 1 447 

21 = 

0.2259 

RI^ 

0.1 4 6 0 

A S. 

0 . 2265 

Rl = 

0. 147 A 

2 i = 

0.2312 

RI- 

0. 1466 

2 1 = 

0 . 2339 

RIs; 

0. 1499 

2 1 = 

C . 2367 

RI« 

0. 15 10 

T 

jL 

0-2394 

Rl = 

0. 15E2 

21 = 

0 . 2422 

RI« 

0. 1530 

i. 1. = 

C.2A5C 

RI = 

0. 154? 

2 1 = 

0 . 2A 7c 

RI = 

0. 1552 

2 1 = 

C . 2507 

RI- 

0- 15c 1 

21 = 

C . 2535 

RI« 

0- 1570 

2 I ”■ 

0 . 2564 

Rl« 

G. 157S 

21 = 

0 . 2593 

RI=; 

G. 15S5 

21 = 

0 . 2622 

RI= 

0- 1592 

21 = 

0 . 265 1 

RI= 

0. 1599 

21 = 

0.26S1 

RI = 

0.1605 

21 = 

0 . 27 10 

RI« 

G. 1610 

21 = 

0 . 2739 

RI = 

0.1615 

21 = 

C.2769 

RI = 

0 « ’1619 

21 = 

0 . 2799 

RI= 

0.1 623 

2 j.. = 

0 . ES29 

RI- 

0 . 1 626 

Z 1 = 

0 . 2856 

RI~ 

0.1 629 

ZI = 

0 . E3aS 

RI = 

0.1631 

21 = 

0.29 16 

RI*=- 

0. 1632 

Z 1 = 

G , 2946 

RI^ 

0- 1633 

21 = 

C , 2976 

RIA*^ 

0 

o 

‘1 2 1 A 

= 0 . 30; 

Rl«- 

0 - 1 634 

Z 1 = 

0 « 3006 


DO YOU WHr.iT TCi CALCULATE MORE EGUIP GTEWTIAL SURFACEIF YES 
t 

GIVE THE I!''4IT1AL FGI.i4'! OF 1 HE DESIRED EsUIF JTENTIHL SURrFVv 
0 m GOOD G - 1 oOv'. 

GIVE THE rr A TEf4'l I Gr.i F'uX.Ur Gr THE EQUlPuTENTiAL SURFACE 
G . GOOD C « '1600 

GIVE THE CCiNSTAr^TS OF THE FIRST BOUraDARY LINE :aM+by+c==^C 
STAF'TILG POINT FAi-LS IN THE REGlUiX - by -rcT AC- 



TIVE 7 r-, 7 67-7'; E GT" TrlE S£7;Or4r' BOUND' Hr- 7 LlNu ^ a rH-Ly'i" c — s 

S T A r. 7 1 rs r' j , i 7" i ' t-i ... i... - N T h E R i.. U i. •.7- ‘ i-. i 7' y ■ r ' —U 

A=- C.G'.-'IC- 


Th’Er4 WRIT 


SUCH THAT 



THE 



ELSE WRITE—0' 


C= -0.3000 


GIVE THE RELATIVE PERMEABILITY OF THE INITIAL POINT 



COCO 



GIVE 

THE' L£!h 

igth 

OF INFINITSIMAL ELEMENT ALONG THE EQUIPQTENTIAL SURFACE 

0. 

0G3G 



RI = 

0 . 0030 


0.1600 

Rl== 

0 . 0060 

7.1- 

0.1601 

RI = 

0.0090 

hsi 

1! 

0.1 60E 

RI=- 

0.0120 

zi==^' 

0.1604 

RI= 

C.0150 

21 = 

0.1606 

RI= 

0.01S0 

Z I- 

0.1609 

RI= 

0 . 0209 

2I~ 

0. 1613 

RI-- 

0 . 0239 

21== 

0.1618 

RI= 

0.0269 

21 = 

0. 16E3 

Rl=- 

0,0298 

ZI = 

0. 16ES 

Rl = 

0 . 0327 

21 = 

0. 1635 

RI--=- 

0.0357 

21 = 

0. 164E 

RI = 

0 . 0386 

21 = 

0. 1649 

RI = 

0.04 15 

ZI = 

0.1657 

RI = 

0.0443 

21 = 

0 . 1 6)66 

R 1=' 

0.0472 

2 i = 

0. 1q75 

Rlr.- 

0.0500 

21 = 

0. 1685 

RI=^ 

G . 0526 

21 = 

0 . 1 696 

RI = 

0.0556 

21 = 

0 . 1 707 


0.0584 

21 = 

0 . '1 7 1 9 

RX=r 

C . 06 1 1 

21 = 

0 . 1 73 1 

RI- 

0.0638 

ZI = 

0.1744 

RI*:^ 

0 . 0665 

21 = 

C . 1 758 

RI^ 

0 . 069 1 

Z I = 

0. 177E 

RI = 

0.07 17 

Z I = 

0 . 1786) 

RI*^ 

0.0743 

2 I = 

0 . 1 SC) 1 

RI^ 

0 - C7 69 

21 = 

0 . 1817 

RX=- 

0 . 0794 

Z 1 = 

0 . 1 833 

RI = 

0.GS19 

Z 1 = 

0. 1850 

RI = 

0 . 0843 

21 = 

0.1868 

ei=:, 

C . 0S6S 

21 = 

0. 1SS5 

RI=^ 

0.0891 

21= 

0.1904 

RI = 

0.0915 

21 = 

0 . 1 9E3 

RI=^- 

0 . 0938 

21 = 

0. 194E 

RI = 

0 . 0960 

21 = 

0. 196E 

RI=^ 

0 . C9S2 

ZI = 

0 . 1 9SE 


0. 1004 

Z 1 = 

0 . ECG3 

R 1=' 

0. 1025 

21 = 

0,E0E4 

RI = 

0 . 1046 

21 = 

0 . E04t.. 


0 . '1066 

21 = 

0 » E068 

RI=' 

0, 1086 

21 = 

0 . EC90 

RX=^ 

0. 1 105 

21 = 

C.E1 13 

RI.-Z 

C - 1 1 2:4 

ZI = 

0 . E 1 3 7 

Rl=-- 

0. 1 142 

21 = 

0 . E 1 6 C 

RI== 

0 « 1160 

21 = 

0.E1S5 

RI=- 

0 -1 •} 7 V 

Zl = 

0.EEG9 

RI = 

C . 1 194 

ZI = 

0 . EE34 

Rl = 

0 . '12 1 1 

Z 1 = 

0.EE5? 

RI = 

0. 1226 

21 = 

0 . 2ES5 

RI-- 

0 « 124 1 

zx= 

0 . E3 1 1 


. 1 L t> 6 

Z J. = 

0 . Eow 

RI==- 

U . 127 w 

2 “ 

O.ESoC 

R I “ 

0» ■■64 


. £390 

R i. "" 

0 . »CL. Z 

1 ~ 

G . E 1 

R !-■ 

0 . 1 : ■: - 

ZI = 

0. 3444 

FU, 

C . . 1 1 

5 1 = 

■6 L4‘:’£ 


0 . ■ 6. 1 1 

1 i = 

0 . E500 





RI = 

0. 1342 

21 = 

0 . 252S 

RI = 

0. 1352 

21 = 

0 . 2556 

RI^ 

0. 1362 

21 = 

0.2585 


0 . 1 37 C 

21 = 

0 . 26 1 3 

Rl = 

0 . 1 379 

21 = 

0 . 2642 

R 1 = 

0 . 1 3S6 

21 = 

0 . £67 1 

RI = 

0 . 1 393 

21 = 

0.2701 

RI=- 

0. 1400 

21 = 

C . 2730 

RI = 

0. 1405 

21 = 

0.275-9’ 

RI = 

0,1410 

ZI = 

0.2789 

RI= 

0. 14 15 

21 = 

G.2819 

RI=. 

0 - 1419 

21 = 

0 . 2S4S 

RI = 

0. 1422 

21 = 

0 . 2S7S 

RI = 

0. 1424 

ZI = 

0 . 2906 

RI^ 

0. 1426 

21 = 

0.2936 

Rl = 

0. 1428 

ZI = 

0.2966 

RI = 

0. 1428 

Z I = 

0 . 299S 

RIA=^ 

0 » 1 428 2 1 A= 

=■ 0.30 GO 

RI=^ 

0 . 1 426 

II 

0 . 3026 


DO YOU WhM i TO ChLCULhTE MORE EGUIFDTENTIAL SURFACEIF" VES THEN WRITE - I- ELSE WRITE-O- 
1 

GIVE THE INITIAL POINT Or THE DESIRED EQUIPGTENTI AL SURFACE 
0-0000 0.1S00 

GIVE THE EXTENT lOrj POINT OF THE EQUIPOTENTIAL SURFACE 
0-0000 0.1600 

GIVE THE CONSTANTS OF THE FIRST BOUNDARY LINE : aK->- &v-i-c=0 SUCH THAT THE 
STARTirL::i POINT FALLS IN THE REGION s aK + bv+c-:>0 
A^ 1 . OCOO 
B- 0.0000 
C- -0-6*000 

GIVE I HE". CGN31 P'tN I S OF THE SECOND BOUNDARY LINE s ax + by'i'C^O SUCH THAT THE 
STARTING F’CINT FALLS IN THE REGION s by+ cT—C 
A= 0 . 0000 
1.0000 
C« -0.3000 

GIVE THE RELATIVE PERHEABXLITY OF THE INITIAL POINT 
1 • OOGO 

GIVE THE LENGTH OF INFINITSIMAL ELEMENT ALONG THE EQUIPOTENTIAL SURFACE 
0,0030 


RI = 

0 , 0030 

2 1 = 

0. 1800 

RI = 

G . OOdC 

21 = 

0 . 1 80 1 

RI = 

C.OGFC 

Z I = 

0 . 1602 

Ri=: 

0 . 0 1 2 0 

2 1 = 

0-1 Ei04 

RX= 

0.01 50 

21 = 

0 . 16'07 

Rl = 

0.0180 

21 = 

G- 1E1 1 

Rl = 

0.0209 

2 1 = 

0 . 1816 

RI = 

0.0239 

21 = 

0. 182 1 

RI = 

u ■ OcoE.’ 

1 

C . 1 6-27 

RI = 

0 , 0297 

2 1 = 

0. 1633 

RI = 

0 . 032c; 

2 1 = 

U - 1 6 4 1 

RI= 

G . 0355' 

i.. i. -■ 

0 « 164 7 

RI = 

0 . 03^84 

1 = 

0. 1658 

RI = 

L « G *^1- 1 2 

2 1 

0 . 1 667 

RX = 

0 . 0 4' 4' 1 

2. 1 

0. 1877 

RI = 

0 . 0469 

Z 1 =- 

0. 1cir-c 

1! 

0 . G4'-- c 

4L -• 

G . 1 9 G G 

R 1 = 

G . C I . *■ 

I . ■■ 

C . 1 9 • £ 

Rl = 

O' . 0 . . 1 

^ J *“ 

G . l9'£.. 

R I = 

0‘ « 05"' 

II- 

G., VFGF 

R I 

0 . 

I j 

I-' . 1 9 5 3. 



R 1 “ 

C.C6S 1 

Zl=^ 

0 . 200C' 

RI- 

0 . 07 C'U 

ZI- 

C . EO 1 7 

RI = 

C.G73G 

z :.~ 

C . 2034 

Rl=- 

G » C 7 1-‘ 4 

Z 1~ 

G.2Q33 

RI- 

C . G777 

Z 1 -~ 

0 ■ 2G7 1 

R I~ 

G .. 0600 

Z I~ 

0 . cG9C 

RI- 

C . 06c_G' 

2 1 - 

0-21 1 G 

RI = 

0 . 064 O' 

Z I “ 

0.2130 

Rl= 

0 . OS 6 7 

Z .1 - 

0 . E 1 5 1 


‘ G.0S66 

71“ 

0.2173 

RI- 

0 . 0906 

Z I =' 

0.2194 

RI = 

G.09E6 

Z 1 -- 

0 . 22 17 

RI=^ 

G . C 9' *+ , ’ 

Ak j. “• 

G . 2240 

RI==- 

0 , 0966 

Z I “■ 

0 . 2263 

RI- 

C . C9S4 

Z I 

0 . £2o7 

RI-^ 

G, IGOE 


G . 23 1 1 

RI^- 

G. 1019 

Z I” 

C . 2336 


0. 1037:' 

Z i -■ 

C « 236 1 

RI^- 

0 . 

Z 1 = 

3 . 23&'" 


C. iCwjCi 

Z I “ 

Q • Z-4 1 3 


0 . 1 06 1 

Z 1 “ 

0 . 2439 

RI=^' 

C . 1 096 

Z I 

0.2460 


0.11 06 

Z I 

0 . 24 92 

RX = 

0.11 S 0 

Z I 

G « 2.!i;'2G 


C . 1 1 36 

21- 

0 . 234 ; 

Rl--^ 

G • 1 143 

Z 1 

0.237Z" 

RI^- 

0. 1 ’(1)4 

2 1 “ 

0 . 26C3 

RI='- 

0 » 1 1 i- ‘-1 

Z j 

C.EtoE 1 


0. 1 ■<73 

A — 

o.Ls.: 

R i, ~ 

C' » 1 1 G " 

Z 3, 

Z « c ' 

Rl = 

0. 1 169 

Z 1 -*■ 

0.2718 

RI=^ 

0. 1 19to 

2I=' 

C . 27'47 

RI = 

O 

h) 

o 

W 

ZI~ 

0 . £77*6 

RI=- 

0. 1c07 

Z I “ 

0 • 2806 

RI = 

0 . 1 2 1 E 

Z I — 

0 . 2636 

RI=- 

0- 1216 

Z 1 = 

0 » 2863 

RI^-' 

0 . 1 2HG 

A. J. “■ 

0.2895 

RI=- 

C. 12Z£ 

A. J. “ 

0 - 2925 


0 . 1224 

Z 1 — 

0 . 2955 

Rl = 

0. 1£E3 

Z 1 = 

0 . 2985 

RIA- 

0 . 1 £ U 21 H 

“ 0 • 3000 

RI=- 

0. l2c!o 

Z 1 — 

0.30 15 

DO YOU IvmY. i 

1 u CAi.,.CULATE 

WR I 

1 

TE 1 EL L 

/i. i'vP. 

ri'E, 0 

C. C' 1 U U 3 c-t t i f. Ij 

tJ-.r 

fe'aui po'len' 

GJ'.’E 

THE liU.I 

' T I AL 

POINT OF 

V/ . 

V Vi' \J W \J • 

C 0 ij G 


ul v‘E. 

1 HE. L .* 1 

t 61 > 

gIN r'ClIOiT i 

0 , 

(.‘G‘. V. Cl, 

ECO.: 


ul VL 

1 HE 3 i-u 

ISr. "1 A1 

iT’6 ur THE 

STAR 

1 INu PEG 

, Ml r 

ALL 3 XN Tl 

A=' 

1 , OC'CG 



C^. 

0 « G 7'., 



'3 1 ■■■ ■:. 

1 Hi-; ‘.6i' 

YZ": r'li' 

Z' C'Z T'rZZ, 


TrC EQI'IPDTENTIAL SURFACE. IF YES THEN 


r bv*i"C- 


SUCH 1 HAT TH 




MH, 1 t™ .i 


ci: ic 


CF jrjrir.l'’ SIMA,.. ELEMENT ALC'vIG THE EaUIPOTENTIAL SURFAC-E 


ilV 3 i '_'X r X/ i » X. Mt. 


URFACEIF YES TKEK WRITE - 1 - ELSE WRITE- 


rEf^TlAi- SURFACE 


n* ' L*r ’*( HL E' 3 v ’ 1 r “ * £ I hl surface. 

1 HE. ]?'•£.'» J 4 *.'hE i u. li'ic. s ij V** c cvUCH THAT THE 

IK Tr:r lu::., 10!.lui.'^-bv4c...-^=F 


.,v i..= l:H„tv^KI,'Mr\’V L. I uL I c. 4 - Ly 4 C "0 SUCH THAT THE 


Jb^ i.Uuv'w 

C=^ - 0 . 39 S 0 

GIVE THE RELATIVE FERflEAEU L I TV OF THE INITIAL POINT 

1 . 000 c 

GIVE THE length OF INF INI TSI ELEMENT ALONG THE EGUIFOTENTIAL SURFACE 


f’\ i w , U u 



( , , ‘ f ' 


f * •. i . 

h ' > , jy ^‘' ? 

f* , t . ui. 

^ I U • U 4i 

^ * » ( jCi -* 

f' i C . %»<&t • 

;*. 'W # W '^ ' *’ 

hi 

K 1 ■'-' u » 0 ■*' w / 
PI- C.u7£« 

r .o:«i 

H J - 
hi 
u I 
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R I ‘ 
h i «' 

HI 
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hi 
H 1 *»■ 

H 
K 
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Rl* 

RI» 

Rl,a 

Rl «. 

Rl «. 


, % G.i . IS 
^ i *' i , r 1 i' 7 

;:^ ( ... 147 

* ; 0.: 16^ 

7 i .' Kjm*^ lOU 
7 1* «£ "hi 

7 1 0 m £tw4 

A '"“ h » 1 S 

;]- O.Zcc-B 
::■ 0.EH41 

i. 1 * ijmC.C.Zd^ 

71'.: 0-E£6fc 

7i- O-cSSE 
71s- 0»£29t 

7 1 »" 0 - £ 3 1 C 

21* 0 . £ ac 4 

2.1* C.t33^ 

21 0.23!i4 

4 1 ‘ 0 • 4*C» / U 

21* ()*£35t> 

2 1-*- 0-240 1 

2 i 0 - £ 4 1 7 

0 « £434 
G.£4!)G 
0.2467 
0 . E4S4 
0 . 21)0 * 
0-2S19 
0.2036 
C.2SS4 
0-E57E 
0-2590 
0-2606 
Op £627 
0-2645 


U.O '46 
C' . C *ti. 

0 - C ' 2 * ^ 
C*-C/'S6 
♦ OoviO 

V » . 0£. 1 1 

y , 0624 
6-0636 21* 

ij.Gt*4.' 2 1 

u*063<' 2i*' 

O-OScS ;j*' 
I* 0.01.6*4 21» 

r o.cea-/ 

0«0c’‘»c> 2i*“ 
w - 0 9 w / 21* 

0*0916 21* 
0*0925 21* 
0-0933 21* 
0.0941 21* 
0.0946 21* 


Rl«» 0-0955 21* 
Ri«: 0.0962 21* 
Rl«- 0.096S 21* 
Ri** 0.0974 21* 
Rl*a. 0.0979 21» 
RX« C.0984 21“ 
RI*^ 0.0989 ZX* 
Rl-ik 0.0993 21* 
Rl». 0.0997 21*» 
fU- 0.1001 Zl* 
R i » 0*1 004 21* 

lU* 0. 1007 21* 
Rl* 0.1009 21* 

RX«4 0. 1012 21* 

Rl* 0. 1013 ZI* 
Rl* 0. 1014 21“ 
Rj;. 0- 1015 21* 

r '{ i ” 0*10 1 6 4- i 

K X *■ V(’ * * W I Cf 1 """ 


0.2664 
0.2683 
0-2702 
0.2721 
0.2740 
0.2760 
0.2779 
0.2798 
0.2818 
C.283S 
0.2857 
0.2877 
0.2697 
0.2917 
0.2937 
0-2957 
0 - £977 
C . £997 
C . 3 0 1 7 
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